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ABSTRACT 


The kinetics of diffusion in cubic metals is treated in terms of statistical 
mechanics from the point of view of absolute rate theory. An approximate 
description of the activated state leads to the conclusion that in the two 
degrees of freedom orthogonal to the jump direction only the lower vibrational 
energy levels are occupied. This circumstance gives rise to a negative contri- 
bution to the experimental entropy of activation for diffusion in solids which 
is negligible for the vacancy mechanism but very appreciable for the ring 
mechanism. It is suggested that a low D,, perhaps of the order of 10-4 em2/sec, 
is an experimental criterion for the 4-atom ring mechanism and that the 
self-diffusion in pure chromium and uranium (b.c.c.) may occur by this 
mechanism. Further, a low Arrhenius frequency factor is predicted for all 
thermally activated processes in solids which require simultaneous jump 
of n atoms, e.g. in certain movements of dislocation loops. 


§ 1. INTRODUCTION 
ZENER (1951, 1952), Wert and Zener (1949) and Wert (1950) have con- 
sidered the diffusion of both interstitial and substitutional atoms in cubic 
alloys in terms of absolute rate theory. ‘They conclude that the diffusion 
coefficient appropriate for use with the activity-gradient formulation 
(Hollomon e¢ al. 1948) of Fick’s first law is given by 


D=oa*vexp(AS/k)exp(—AH/kT) .. . . (1) 


where a is the lattice parameter, v is the Einstein vibrational frequency of 
lattice-site atoms, and AS and AH are the entropy and enthalpy of acti- 
vation, respectively. The temperature-independent factors are usually 
grouped together to express 

= OEVORPNNO Keser ae 2) ted ea (2) 
The numerical factor « is determined by the geometry of the atomic jumps. 
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Ignoring correlation coefficients, it is defined by 
(4)Z7? = aa? So een 


in which Z is the number of jump paths of length /. For diffusion of 
substitutional atoms by the vacancy mechanism, « is unity in both b.c.c. 


and f.c.c. lattices. é ; 
AH is difficult to calculate and is usually determined experimentally 


from the slope of an Arrhenius plot. The AS, which may be determined 
experimentally from the Arrhenius intercept, is presumed to be somewhat 
more susceptible to theoretical approaches. Zener (1951) invoked the 
thermodynamic relationship 


AS = —(dAG/OT),, eS Roe ne) 
and evaluated the derivative on the assumption that elastic shear strain 


is the only temperature-dependent contribution to the Gibbs free energy 
of activation AG,,, for motion. The resulting expression is 


ASS = —AAH, d(plig)|@l oe eee ee od 
where ,« is the measured elastic shear modulus and 
A= (AG°/AGo)m = (AHy’/AHo)m = (AH*/AH),, : : oO (6) 


in which the superscript e refers to the elastic contribution and the sub- 
script 0 to absolute zero of temperature. Thus, in view of the negative 
temperature coefficients of elastic moduli, it was suggested that entropies 
of activation for diffusion might always be positive. However, as pointed 
out by Huntington et al. (1955), the above treatment neglects dilatation 
of the activated species. This contribution to AS may be large and of 
either sign. Accordingly, the possibility of a negative entropy of acti- 
vation in diffusion is not excluded by thermodynamic arguments. 
Nevertheless, most experimental values of AS for self-diffusion in cubic 
metals are positive and range from 4 to 14e.u. (Zener 1951), and, with the 
exception of tungsten, the observed D,’s are of the order of unity. 

Zener (1950) calculated that the potential energy of activation for the 
ring mechanism is a minimum for four atoms in the case of copper (f.c.c.). 
Thus it was conceived that diffusion might occur by the rotation of 
4-membered rings. However, the estimates (Huntington and Seitz 1949, 
Zener 1950, LeClaire 1953a,b, Huntington 1953, Fumi 1955) of the 
potential energies of activation for the various mechanisms of diffusion in 
f.c.c. lattices indicate that the vacancy mechanism is favoured. Also, the 
Kirkendall marker shift, which demonstrates the presence of the vacancy 
mechanism, has been found in many f.c.c. substitutional solid solutions. 

The calculation of the potential energy of activation for ring diffusion 
in b.c.c. lattices has not been made, but it is considered (LeClaire 1953 a, b) 
that the ring mechanism might be favoured in some of the less closely 
packed lattices. Although Kirkendall shifts have been shown to ac- 
company chemical diffusion in b.c.c. solid solutions, e.g. beta brass 
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(Landergren et al. 1956), titanium—molybdenum (Shewmon and Bechtold 
1955), potassium—rubidium (Yang et al. 1958), uranium—molybdenum 
(Adda and Philibert 1958) and iron—chromium (Paxton et al. 1959 b), it 
should be borne in mind that these measurements were conducted on 
fairly concentrated alloys, and it is possible that the ring mechanism 
might be operative for self-diffusion of the solvent atoms in more dilute 
solutions. 

Zener (1951) evaluated a for 4-atom ring rotation and found that it is 
two for f.c.c. and six for b.c.c. lattices. Also he inferred that the elastic 
strain contribution to the entropy of activation should be positive and 
probably of the same order ofmagnitude as that for diffusion by the vacancy 
mechanism. ‘Thus he concluded that the apparent entropy of activation 


ASapp=* In (D,/a?v) . ° ° . . . . (7) 


for diffusion by the 4-atom ring mechanism should be positive and perhaps 
larger than that for the vacancy mechanism. 

From a somewhat more fundamental point of view, Huntington and 
co-workers (1955) considered the entropy of ring formation in terms of 
changes in vibrational frequency upon activation, neglecting the partition 
functions of the atoms in the directions for rotation. For copper (f.c.c.), 
they calculated the entropy of formation of an activated 4-atom ring to be 
—4e.u., in contrast to + 2e.u. for the activated species in diffusion by the 
vacancy mechanism. Significantly, they found the effect of Debye 
coupling to be negligible. However, this evaluation of the entropy of 
activation assumes that the energy levels of an activated atom in the two 
degrees of freedom orthogonal to the jump direction are those of an 
oscillator which is stationary with respect to the reaction coordinate. In 
the following discussion we will refer to these as the normal vibrational 
energy levels of the activated atom at its saddle-point position. One 
notes that these are not the levels of the atom in its lattice-site position, 
which are somewhat lower and more closely spaced. 


§2. An ALTERNATIVE DESCRIPTION OF THE ACTIVATED STATE 


The present authors note that the above assumption may not be ap- 
propriate for evaluation of the partition function of the activated species 
in diffusion in solids. The difficulty is that the time of existence of the 
activated atom at the saddle point is so brief that a precise description of 
the energy states is incompatible with the quantum-mechanical un- 
certainty principle (Hirschfelder and Wigner 1939). More succinctly, the 
transition state cannot be defined because the curvature of the potential 
energy surface is too great to satisfy the condition 


(h[2m)(K/my2<kT  . . . . . . . (8} 


where K is the ‘spring constant’ and m atomic mass. This difficulty 
might be obviated by considering the moving atom as a non-localized. 


2H 2 
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particle in a periodic potential. In treating the problem of a particle in 
a sinusoidal potential, Slater (1952) showed that for energies less than the 
height ¢ of the potential barrier the solution may be approximated by 
that of a localized linear oscillator without appreciableerror. — Significantly 
he found that the levels for energies slightly greater than e are not at all 
those of a freely translating particle as assumed for the degree of freedom 
in the reaction coordinate of standard absolute rate theory (Glasstone et al. 
1941). In fact, a forbidden gap is present and only when the energies are 
appreciably greater than « may the assumption of a freely translating 
particle be invoked. Unfortunately, the results of the Slater treatment 
are in a form which is not readily amenable to application in absolute rate 
theory. One notes that an exact treatment would involve a wave- 
mechanical approach, not just for one particle as carried out by Slater 
(1952), but for all the degrees of freedom of all atoms in terms of a normal- 
mode description of cooperative motion. 

From the above discussion it would appear difficult to provide a direct, 
quantitative description of the true activated species for diffusion in solids. 
Accordingly, an alternative treatment is proposed in which the con- 
centration of activated species is evaluated by calculating the equilibrium 
number of those reactant species which are in states favourable for 
reaction and assuming this to be equal to the actual number of activated 
species. In other words, a Boltzmann distribution is used to calculate 
the concentration of species in known vibrational states which are only 
incrementally removed from the true activated state at the saddle point. 
Actually, this is similar to Rice’s (1958) treatment, except that here the 
difficult normal-mode description of cooperative motion of atoms is not 
employed. 

One next considers the most appropriate formulation of the partition 
functions to be used in determining the Boltzmann distribution. The 
treatment of Vineyard (1957) or Rice (1958) in which the cooperative 
motion of the jumping atom and its neighbours is depicted by a normal- 
mode description would be more exact, particularly for low temperatures. 
However, for high temperatures where the long wavelengths are less 
important, it is considered that the Einstein formulation of the lattice 
dynamics in terms of independent oscillators should be satisfactory. 
Also, the Vineyard—Rice treatment is not readily amenable to numerical 
evaluation. Accordingly, due to its simplicity and applicability at the 
high experimental temperatures, the Einstein formulation of vibrational 
partition functions will be used in the present analysis. 


2.1. Self-diffusion by the Vacancy Mechanism 


If AG vac is the Gibbs free energy of formation of a vacancy, the prok- 
ability of its existence at a given lattice site is 


P1= exp (—AG%\ac/kT’). ates : ee (9) 
Similarly, the probability that the j neighbouring atoms are in states 
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corresponding to the saddle-point configuration is 

sp OXD NG conjkl je ak ee OY 
The probability that an adjacent atom has sufficient vibrational energy in 
the given x direction to surmount the saddle point for motion to the vacant 
site is i: 

Ps= 3 exp(—Hy/kT)/fyzexp(—e«/kT)  . . «. (11) 
where « is the height of the potential energy barrier and f, is the Einstein 
partition function for vibration in the one degree of freedom. 

The compound probability p,p.p, is not sufficient to express the likeli- 
hood that the atom is only incrementally removed from its saddle-point 
state, because the actual motion of the atom may be in a direction other 
than that leading to the vacancy. In this case, the energy levels of the 
activated species would be appreciably elevated. Thus the predominant 
reaction path is the axis between vacancy and jumping atom. The 
condition for the atom to be moving along this reaction path, i.e. in exactly 
the right direction, is that it have almost no vibrational energy in the two 
orthogonal directions. The probability that an atom has the least possible 
vibrational energy, i.e. the zero-point energy «,=/v/2, in one of its 
orthogonal degrees of freedom is 1/f,, and the probability of zero-point 
energy in both of its orthogonal degrees of freedom is 1/f,2+. However, 
depending on the mode of interaction with its neighbours near the saddle 
point, the atom may be permitted some small vibrational energy in the 
orthogonal directions and still make the jump. That is, if the orthogonal 
vibrational energy is sufficiently small, the energy levels of the activated 
species will not be appreciably elevated. This critical orthogonal 
vibrational energy is difficult to evaluate, but one presumes that it is 
smaller in hard lattices than in soft ones where a misdirected jumping 
atom may push aside its saddle-point neighbours. Accordingly, for 
the case of very hard lattices (high Debye @) only, the probability of 
motion in the reaction path is 

Pa= lf? =exp (—AG4/kT). cee rere Ch 2) 

This term does not arise in treatments of diffusion by standard absolute 
rate theory which assume normal vibrational energy levels for the acti- 
vated atom in the degrees of freedom orthogonal to the jump direction. 
Retaining the concept of vibrational energy levels, the presence of this 
factor suggests that these levels are much higher for the activated species 
than has been supposed (Huntington ef al. 1955) and hence only the 
lower ones are occupied. In other words, high orthogonal vibrational 
frequencies of the activated atom are implied. This view is consistent 
with the idea that a proper wave-mechanical approach to the problem 
would yield high and widely spaced energy levels for the activated species. 
In this paper the term activated species is not restricted to a description 
of atoms that are in free translation over a potential energy barrier, as 


a 
+ In order to simplify the notation, the zero-point energy is taken as the zero 
datum of energy in evaluating the partition functions. 
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it usually is in standard absolute rate theory. Rather it refers to all 
species which are crossing the barrier when the system is close to 
thermodynamic equilibrium. In this connection, one notes that the solid 
system in which diffusion is occurring under usual experimental conditions 
is quite close to its thermodynamic equilibrium state. 

The probability that the atom has sufficient vibrational energy to sur- 
mount the potential energy barrier and motion in the right direction (1.e. 
that it is only incrementally removed from its saddle-point state) in 
lattices of high 6) is 


DP =P1PoPsP a= XP [ — (AG vac + AG con + €)/kT ] I/f?. (13) 


Assuming that p is the probability of occurrence of activated species and 
recognizing that the transmission frequency* is, by this model, the Einstem 
vibrational frequency v, the actual atomic jump frequency becomes 


Peivp>~ 2a. gs 2, eee 


The diffusion coefficient is obtained from the standard expression for 
three-dimensional random walk as 


D=($)Z2T = aa?vexp [ — (AG vac + AG con + €)/AT A/F? (15) 


for diffusion by the vacancy mechanism in hard cubic lattices. Upon 
differentiating the logarithm of eqn. (15) with respect to the reciprocal of 
‘absolute temperature, it emerges that the experimental enthalpy of 
activation is 


AH = ABM\a¢ + AH cont € + 2k[dln fy/d(1/T)}. «a oO 
From eqns. (1), (15) and (16), the experimental entropy of activation is 
given by 

AS = AS%\ac + AS°con— 2hk[In fy +(dInfy/dinT)]. =... (17) - 

In view of the fact that only the terms containing f, arose from the 

requirement for motion in the right direction, it is clear that the probability 

for such motion contributes the following quantities to the enthalpy and 
entropy of activation: 

AH 4=2k[dlnf,/d(1/T)], st OO acl aie pace emia ee 

AS = —2k[Inf, + (dInfy/dIn 7)}. oe i SO ie 

Alternatively, eqns. (18) and (19) may be obtained directly from eqn. (12). 

In general, AH; and AS, are negative contributions, and they are 


readily estimated on the present assumption of independent oscillators. 
The Einstein vibrational partition function for one degree of freedom is 


. f= & exp[- vik] = [1 ~exp(—hv[kT)7, (20) 


aalAsas bt he transmission coefficient is assumed to be Unity. 6 oo 


IB AT Bache sg 
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and at high temperatures where kT’ >hv, fy2kT/hv=T/6,. Subject to 
this approximation, eqns. (18) and (19) become 


AHg= —2kT ee ol) 
ASH Sok ln (1Y0,)]. se. ce ee (22) 


At the usual high experimental temperatures (~~ 1000°K) for self-diffusion 
in solid metals, fy ~ 7'/6,, is often about 3. Accordingly, in a typical case 
of self-diffusion by the vacancy mechanism in cubic metals of high Debye 0, 


AH ,~ — 4000 cal/mole and AS, ~ —8e.u. 


One notes that the sign and orders of magnitude of these contributions 
are consistent with our physical model of the diffusion process; they 
represent the losses in energy and entropy occasioned by the inaction of 
two degrees of freedom. Thus the contributions from the probability of 
correct direction of vibrational motion to the enthalpy and entropy of 
activation are not very large in the case of the vacancy mechanism where 
only one atom undergoes net movement. 


2.2. Self-diffusion by the Ring Mechanism 
Analogously to the treatment outlined above, it is readily shown that 
the expression for the coefficient of diffusion by the n-atom ring mechanism 
in hard cubic lattices is 


D=aa*v exp [— (AG®con+ ne)/kT 1 / fy” eer 23) 
in which the probability of correct direction of motion for the n atoms is 
identified as 

pa=1/fy" =exp (— AGg/kT). 5 AST Go NOLS) 


Also, the analogues of eqns. (18) and (19) for the contributions of pz to 
the enthalpy and entropy of activation are 


AH j= HE 1/T)] on ata ie inci Boge 1 0455 
and 
AS,= —2nk{in fy + (dln fy/d1n T)]. Smet 20) 
Again, invoking the approximation that fy ~ T'/@, 
| AH , = —2nkT Wey ceed a Does (2:7) 
and 
ASqg= — 2nk[{1+1n(T7/@5)] RAY PP I( 28) 


Accordingly, for the interesting case of self-diffusion by the 4-atom ring 
mechanism in cubic metals of high 0) in a typical situation where 


T = 1000°K and 7/6, ~3, 
~ AH, & —16000cal/mole and AS,~ —32e.u. 
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Thus, the probability of correct direction of vibrational motion, which 
has a relatively small effect in diffusive processes involving-jump of only 
one atom, becomes very important in processes involving simultaneous 
jump of several atoms. This circumstance arises from the compounding 
of probabilities for correct direction of vibrational motion for all of the n 
atoms. As discussed in the section on self-diffusion by the vacancy 
mechanism, the requirement of correct direction of atomic motion may be 
viewed, from a more fundamental standpoint, as signifying high vibrational 
energy levels for the activated atom in the two degrees of freedom per- 
pendicular to the direction of jump. 
It is obvious that exp (AS,/k) is by no means cancelled by 


exp (— AH kT); 


in fact, exp(AS,/k) is by far the more significant term in the expression 
for the diffusion coefficient. Therefore one chooses to establish a criterion 
for diffusion by the ring mechanism in terms of D) rather than AH. From 
the discussion presented in the Introduction, it seems reasonable to 
presume that the contributions of AS ac and AS®oon to the entropy of 
activation may be neglected in an approximate comparison of the D,’s 
for diffusion by the respective mechanisms of vacancy and n-atom ring. 
Accordingly, the expected ratio of D,’s for self-diffusion in cubic lattices 
of high @p is obtained from eqns. (2), (22) and (28) as 


Do ring/Do vae = (oring/vac)(@p/eTZ’)?"—, sf wiper epee 
which has a value of about 10-> for n=4 in b.e.c. lattices. Since the 
observed values of Dy for self-diffusion in most cubic metals, which is 
presumed to occur by the vacancy mechanism, are generally of the order 
of unity, the present authors suggest that an appreciably lower value of 
measured Dy, perhaps of the order of 10-3 to 10->cm?sec~, is an experi- 
mental criterion for diffusion by the 4-atom ring mechanism. Further, 
the ring mechanism should be operative in cases where it has a sufficiently 
low AH to compensate for its low Do, and, in view of the Arrhenius plot of 
In D versus 1/7’, low temperatures should favour the ring mechanism. 


§ 3. COMPARISON OF THE PRESENT THEORY WITH OBSERVATIONS 


Paxton and Gondolf (1959) and Paxton e¢ al. (1959) recently studied 
the self-diffusion of chromium (b.c.c.) from 1225 to 1550°x by three 
different radioactive-tracer techniques and concluded that the Dy is of 
the order of 10~*cm?/sec. 7'/0,,~3 for this case}, and from eqn. (28), the 
contribution of probability of correct vibrational direction to the experi- 
mental entropy of activation for four atoms in ring diffusion is 


AS4z —32e.u. 
a ee Se eT ENT 
+ Op =485°R. 
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According to eqn. (29) and the argument of the preceding section, the 
expected D, for diffusion by the 4-atom ring mechanism in chromium is 
=10~cm?/sec. Therefore it is suggested that the self-diffusion in pure 
chromium may occur principally by this mechanismf. 

Adda and Kirianenko (1958), Bochvar et al. (1958) and Rothman (1959, 
private communication) ef al. measured the self-diffusion coefficient in 
b.c.c. uranium between 1070 and 1310°k and found D, to be about 
15x 10% cm?/sec. The predicted D, according to the present theory is 
~10~*cm?/sec. Although the quantitative agreement is poor, it should be 
remembered that the present treatment assumes a high Debye 6, while 6,, 
for b.c.c. uranium is only 170°K (Rothman 1959, private communication). 
Actually, the departure from the theoretical prediction is in the direction 
to be expected for soft lattices. Thus it is suggested that the 4-atom ring 
mechanism of self-diffusion may be operative in b.c.c. uranium also. 

Mayburg (1955) has given arguments to show that self-diffusion in 
germanium occurs by a ring mechanism. However, the observed Dy is 
7-8cm?sec? (Letaw et al. 1955). Inasmuch as 7'/0, for germanium is 
about 4 at 1100°K, the present theory suggests that, if the ring mechanism 
is indeed operative, the number of atoms in the ring is less than 4. In 
fact, considering the diamond-cubic structure of germanium, direct 
interchange (n= 2) might be expected, for which the predicted D, =~ 10-°. 


§ 4, CONCLUSIONS 


(i) There are negative contributions to the experimental enthalpy and 
entropy of activation for diffusion in solids which may be evaluated 
through an alternative description of the activated state. In this treat- 
ment, the concentration of species in a state incrementally removed from 
the saddle-point state is calculated and assumed to equal the actual con- 
centration of activated species. The requirement that jumping atoms 
have the correct direction of vibrational motion leads to a probability 
Pa=\/fv", in the case of hard cubic lattices. This term gives rise to 
negative contributions to the enthalpy and entropy of activation 
given by AH,=2nk[dln f,/d(1/T)] and AS;= — 2nk(Infy + (dinfy/din T)]. 
Invoking the assumption that the vibrational partition function in one 
degree of freedom f,~7'/6p, these formulae simplify to AH,= —2nkT 
and AS, = — 2nk[1+1n(T/@p)]. 


(ii) These expressions do not arise in standard treatments of diffusion 
by absolute rate theory, which assume the normal vibrational energy 
levels for the activated atom in the degrees of freedom orthogonal to 
the jump direction. Their presence here suggests that the orthogonal 
vibrational energy levels of the activated atom are high and hence only 
the lower ones are occupied. 


+ The interstitialey mechanism is ruled out by energetic considerations 
(LeClaire 1953 a, b). 
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(iii) These negative contributions to AH and AS are unimportant for a 
process of single atomic jump (n=1), such as occurs in diffusion by the 
vacancy mechanism. However, they are important in diffusion by 
the 4-atom ring mechanism (n=4) where AHg~ — 16000 cal/mole and 
AS, —32e.u. in cubic metals of high Debye 6. In general, the contri- 
bution of AS, to the diffusion coefficient is more significant than that of 
NEG: 

(iv) It is suggested that a low Dy, perhaps of the order of 10~*cm?/sec 
(versus unity for self-diffusion by the vacancy mechanism), is an experi- 
mental criterion for the 4-atom ring mechanism. 


(v) The present theory suggests that the self-diffusion in pure chromium 
(b.c.c.) and pure uranium (b.c.c.) may occur principally by the 4-atom 
ring mechanism. 


(vi) There are appreciable negative contributions to the experimental 
enthalpy and entropy of all thermally activated processes in solids which 
require simultaneous jump of » atoms, e.g. in certain movements of dis- 
location loops (Seeger 1956). However, in this latter case, the appropriate 
quantities for AH, and ASq in hard cubic lattices would be just half of 
those quoted in item (i) above, because vibrational motion is restricted in 
only one of the two degrees of freedom perpendicular to the jump 
direction. Nevertheless a low Arrhenius frequency factor is predicted 
_ for such processes. 
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Total Cross Section of Lead for Slow Neutrons+ 
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ABSTRACT 


A beam of filtered neutrons with a mean wavelength of 8-4 & has been used 
to measure the total cross section of lead as a function of temperature from 
290°K to 840°x. The cross section is found to vary linearly with temperature 
in both solid and liquid phases to within 1°x of the melting point. There 
is a Jump in cross section at the melting point of (9-4+0-6)%. The results 
in the solid are not in agreement with predictions based on the Debye theory. 


§ 1. IyTRODUCTION 


In recent years, much attention has been given to the use of slow neutrons 
to investigate thermal motions in solids (for references see Kothari and 
Singwi (1959)). The present experiment was performed to study the 
behaviour of these motions in lead near the melting point, as reflected in 
the variation of its total inelastic cross section. A similar experiment of a 
preliminary nature has been performed by Squires (1954). 

The total cross section was measured for lead over the temperature range 
290°K to 840°K, by means of a beam of cold neutrons from the BEPO reactor, 
Harwell. Observations in the solid phase are compared with the predictions 
of the Debye theory and those in the liquid phase similarly analysed in 
terms of ‘effective’ Debye temperatures. 


§ 2. EXPERIMENTAL 


A collimated neutron beam was passed through a filter of bismuth metal 
in the form of shot (Egelstaff and Pease 1954) in order to remove all 
neutrons of wavelengths less than 6-5A. The resultant spectrum is roughly 
that of a Maxwellian of temperature 350°K with a sharp cut-off imposed at 
this wavelength. The filtered neutrons thus cannot suffer elastic Bragg 
scattering by the lead sample. A correction of about 0:7°% was necessary 
to allow for imperfect filtering action. 

The filtered-beam method has the great advantage (over neutron choppers 
or spectrometers) of high neutron flux resulting in good statistical accuracy 
and short counting times. It has the disadvantage of large incident 
wavelength spread. The total cross section of lead, however, is only a 
slowly varying function of neutron wavelength (Hughes and Schwartz 
1958) over the range used (6-5 A to 134, longer wavelengths being present 


in negligible intensities). 
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The 99-99% pure lead sample is contained in a steel holder, an identical 
‘dummy’ holder being used for ‘out’ counts. Both holders are heated in 
an oven specially designed for precise temperature control (to + 7°K for 
many hours) of the specimen. The effective number of lead atoms in the 
beam was continually measured by observing the attenuation by the sample 
of a beam of epithermal neutrons (obtained by transmission through a 
cadmium sheet). The cross section for these neutrons is independent of the 
structure and thermal motions of the sample: it depends only on the number 
of lead atoms in the beam and on purely nuclear properties which may be 
measured in a subsidiary experiment. ; 

The total cross section of lead is plotted (curve A) against temperature 
in fig. 1. The relation is linear in both the liquid and solid states over the 
temperature range with a discontinuity of 0-239 + 0-016 barns (9-4%) on 
melting. We have investigated the linear relation in some detail near the 
melting point, as shown in fig. 2. It is found to hold to within 1°K of that 
point in both the solid and liquid states. 


§ 3. Discussion 

The total scattering cross section is the sum of a coherent and an inco- 
herent part. The coherent scattering may be further divided into two parts, 
one representing the scattering by a single nucleus and the other the inter- 
ference effect. To a first approximation, the interference term may be 
neglected. This is known as the ‘incoherent approximation’ and implies 
that the total coherent cross section has the same dependence on the atomic 
motions of the scattering system as does the total incoherent cross section. 

For cubic crystals with Bravais lattices (e.g. lead) the total incoherent 
cross section depends only on the frequency distribution function and not on 
the detailed phonon dispersion relations. It can be shown that the use of 
the Debye approximation is quite adequate in this case. Moreover, Placzek 
(1954, 1957) has shown that the total incoherent cross section can be con- 
veniently expanded in powers of the ratio of neutron to nuclear mass (mass 
expansion). Using this expression, which is very rapidly convergent, we 
have calculated the total inelastic incoherent cross section. In the 
present experiment the elastic coherent cross-section is zero for the solid. 
The results combined with the incoherent approximation give the total 
cross-section (o,) Shown as curve B (Fig. 1). 

We now consider the effect of interference on the scattering from the solid. 
Placzek and Van Hove (1955) have shown from theoretical considerations 
that the total coherent inelastic cross section depends on the precise form - 
of the phonon dispersion relations and that the use of the Debye approxi- 
mation will not, in general, give accurate values. The present experiment 
enables us to consider the adequacy of the Debye approximation in the case 
of lead. 

The cross section was evaluated using a mass expansion. This con- | 
veniently divides into two parts: the first (and numerically the dominant) of 
which is identical with the mass expansion of the incoherent scattering 
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Fig. 1 
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cross section, and the second represents a ‘ correction’ term for interference 
effects. Marshall and Stuart (1959) have evaluated this correction term 
on the basis of the Debye approximation for an appropriate range of 
experimental conditions. When the correction term is added to the 
incoherent approximation curve B, we obtain curve C, which thus repre- 
sents the total coherent cross section (og) in the Debye approximation. 
Both og and og are average values over the expected beam spectrum 
assuming the file flux is Maxwellian. 

To calculate the cross section at any temperature and incident wavelength 
we require a value of the Debye temperature, Op. Chipman (1959) has 
measured the variation of the Debye-Waller factor of lead with temperature 
using x-ray techniques. We have taken his values of @) for lead between 
300°K and 600°K, given by 

Ost (98-5-0-045 T)°R) )) 5 oe 
where 7' is the absolute temperature. The errors in these values are about 
4°K. 

The Debye-Waller factor involves a different average value over the 
frequency distribution than does the incoherent cross section, but, using 
more realistic frequency spectra due to Leighton (1958), we have estimated 
that the two ‘effective’ Debye temperatures are consistent to about $%. 

The mean wavelength of the beam, measured by means of a thin 1/v 
gold absorber, was 8-454. The lead absorption cross section was taken as 
0-170 + 0-002 barns at 220m/s (Hughes and Schwartz 1958), with a 1/v 
variation. We have taken the value of 0-05 + 0-03 barns for the incoherent 
cross section of lead as given by Egelstaff and McCallum (1958). This 
figure gives rise to the greatest source of error in interpreting the results. 

Figure 1 shows that curve C is a closer fit to the observations than the 
incoherent approximation. It is possible to obtain a fairly satisfactory 
fit using the Debye approximation if we take 6, as given by the equation 

Gay (170-0407) ee ea 
The errors on these values are about 4°K. Values of 6p given by (2) are 
much higher [(21 + 6)°K] than those given by Chipman and other sources 
(Blackman 1955). We therefore conclude that the Debye theory is not 
adequate for predicting total coherent cross sections. _ 

We assume that the part of the scattering from a single nucleus has been 
correctly calculated, so that the difference (B—A) on fig. 1 represents the 
effects of interference in the solid. If there is no change in the effects of 
interference on melting, the observed rise in cross section would imply a 
fall in 6) of (7+ 2)°K. It seems reasonable to expect these effects to be of 
less importance in the more disordered liquid state, in which case the change 
in 6) would be an indeterminate amount smaller (and may even be of 
opposite sign). In view of the expansion on melting however, it is unlikely 
that @p in the liquid will be greater than in the solid. This discussion is, 

_of course, based on a simple-minded solid-like representation of a liquid 
in which an ‘effective’ 6, describes the atomic motions. 
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It should be noted that we have assumed that there will be no quasi- 
elastic scattering of neutrons of ) > 6-54 by the liquid lead. Whether the 
‘cut-off’ is sufficiently sharp to ensure this cannot be decided from existing 
neutron diffraction data for liquid lead (Chamberlain 1950, Sharrah and 
Smith 1953). If there is some quasi-elastic scattering, this would mean 
that the increase in the inelastic cross section at the melting point is less 
than the observed increase, again implying that the fall in 0p is less than the 
value of 7°K mentioned above. 

In spite of all these uncertainties in interpretation, the outstanding 
feature of the results in the liquid still remains: viz. the surprisingly small 
change in cross section, at the melting point. 
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ABSTRACT 


The deformation stacking-fault probability «, has been determined by the 
deformation of bulk specimens of zone-refined lead at 4:2° and 77°x. 
Aluminium (99-996%), «-brass (70:30), lead-0-1 at. % silver and 
lead-20-0 at. % indium were also examined at 4:2°x. It is shown that a 
major difference exists between aluminium and lead, the latter being 
copper-like in its value of «. The influence of increased deformation in 
increasing the value’ of « is demonstrated. Addition of indium to lead 
suppresses a measurable value of « which correlated with observations 
made on twinning in this alloy. 


§ 1. LyrropuctTion 
1.1. Brief Review 


Tue stacking-fault energy, y, is an important parameter in the theory 
of deformation of f.c.c. metals ; (see for example, Seeger (1957)). 
However, it has not been possible to predict y adequately from theoretical 
considerations or to determine y directly for any pure metal. In general, 
it is inferred that y can be related to other measurable properties such as 
the twin boundary energy, creep rate, onset of cross slip or the stacking- 
fault probability, «, etc. ; it is thus possible to compare the values of 
y for different metals by comparing the values of the properties to which 
it is related. 

Data for copper have been used most frequently as a basis for com- 
parison. Fullman (1951) determined the twin boundary energy in 
copper and from this, Seeger.and Schéck (1953) estimated a value 
You = 40 ergs/em?; other data by Fullman (1950) was used to set a lower 
limit for aluminium, y,, 2 200ergs/cm?. In the assessment by Seeger and 
Schock (1953) of the role of y in various deformation processes its value 
was considered high if y>10-*b?/s,,c, and low if y<2x 107? b?/syc 
where c is the spacing of glide planes for which the elastic coefficient of 
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shear is s,,, and b is the Burgers vector of a whole dislocation. On this 
basis, copper has a low value of stacking-fault energy and aluminium a 
high value. 

A comparison of y for different metals may be obtained from their 
behaviour during deformation. Thornton and Hirsch (1958) interpreted 
creep data for several f.c.c. metals in order to determine values of y; 
Haasen (1958) obtained a value for y,; from an analysis of tensile tests 
on single crystals of Ni. However, the conclusions of both these investi- 
gations are based on the value of y,, and the behaviour of copper as 
being typical of a metal having a low stacking-fault energy. Recently 
doubt has been cast on the relative measurements which employ yo, 
as a standard by the tensile deformation analysis of copper single crystals 
made by Seeger et al. (1959). Using a theory relating the onset of 
cross slip to the value of y, they determined an absolute value of 
You = 169 ergs/cm? by measurements of the strain-rate and temperature 
dependence of 7,,,, the stress for initiation of cross slip. This value is 
much greater than the value determined from Fullman’s data. Hence, 
a discrepancy exists. ; 

The measurement by x-ray techniques of the stacking-fault probability, 
a, which is related to y, seems to favour Fullman’s value of y,, when 
stacking-fault probabilities for copper and other metals and alloys are 
compared. Following a theoretical analysis of the x-ray pattern of 
deformed f.c.c. metals, Paterson (1952) predicted that a shift of certain 
reflections should occur in the presence of stacking faults. X-ray data 
may be conveniently described by a parameter « representing the fraction 
of close-packed planes which is in a faulted sequence. (Values of « for 
various metals and alloys are listed in table 6.) However, the true 
configuration of the faults produced during deformation which give rise 
to the shift in x-ray reflections is not known. 


1.2. Materials Chosen for Investigation 


All common f.c.c. metals except lead and aluminium are monovalent. 
Since the electronic structure of the polyvalent elements is complicated 
by Brillouin zone overlaps, it may be expected that a closer similarity 
should exist between the monovalent elements when compared with 
polyvalent elements. This might be reflected in different values of y 
between the two types. However, although y-values for the noble 
metals have been estimated, a good set of values for aluminium and 
lead is not available. In the investigation of aluminium nearly all the 
techniques give only lower limits for y,, Viz.: Twins are found 
infrequently ; and, as will be shown, « is almost too small to be measured ; 
no stacking faults are seen in thin films by transmission electron 
microscopy. 

High values of yp, have been predicted (Seeger (1957), pp. 267 and 
302) although an uncertainty has been expressed (Seeger (1955 a)) 
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The deformation of single crystals by Feltham and Meakin (1957) 
indirectly supports a high value for yp,. Only the determination of 
relative twin boundary energies (Bolling and Winegard (1958 a)) supports 
a low value for y,,~20ergs/em? which is intermediate in the range 
between 10 and 2x 10-*b?/s,,c; however, doubt can be expressed 
about the correlation between twinning frequency and y. The creep 
data of Thornton and Hirsch (1958) was used to predict a value of Vo 
within large enough limits so that it might be considered either high or 
low. 

It became of interest, therefore, to examine the value of ie BY. 
other available techniques. An investigation by means of x-ray 
diffraction is described in the present paper}. Pure lead, two alloys of 
lead, aluminium and « brass were examined. The alloys contained 
0-10 at. % Ag and 20-0at. % In respectively. The 0-10 at. °% Ag alloy was 
the same as that used earlier in the twin boundary energy experiments 
(Bolling and Winegard (1958 a)). 

The Pb-20-0at. % In alloy was selected as representative of a binary 
system based on lead in which the nominal electron-to-atom ratio 
decreases within the large range of primary solid solubility of indium 
in lead. This alloy has thus a lower nominal electron-to-atom ratio 
than pure lead, and it affords a comparison of any resulting change in 
the stacking fault probability with systems where e/a increases on alloying, 
as, for example, in Cu-Zn (Warren and Warekois (1955)). A brief 
investigation of twinning frequency in this alloy was also undertaken for 
comparison between the occurrence of annealing twins, stacking-fault 
probability and stacking-fault energy. 


§ 2. EXPERIMENTS 
2.1. Sample Preparation 


Lead was zone-refined from a starting material of 99-9999°% purity. 
The two alloys were prepared from the zone-refined lead and 99-999% 
pure silver and 99-999% pure indium. The alloy with 0:10at % Ag 
was of a composition which extends beyond the limits of primary solid- 
solubility at room temperature; however, no evidence of a second phase 
was observed either metallographically or in the x-ray study. As 
expected the alloy with 20-0at.°% In exhibited no second phase since 
the primary solid-solubility is more than 60at.% In. The aluminium 
used was of 99-996% purity. 


2.2. The Cryostat 


The cryostat used was similar to that designed by Barrett and 
co-workers (Barrett (1956)). The specimen was held in a copper holder — 
which was attached to the bottom of a liquid helium reservoir and was 


Neen EEE EEE EEE 
+ The deformation of lead single crystals at low temperatures was studied 
' for the same purpose (Bolling et al. (to be published)). 
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surrounded by radiation shields at liquid nitrogen temperature. A 
nickel foil x-ray filter was used as a part of the radiation shield. The 
x-ray windows were of 0-030in. thick beryllium. A pointed stainless 
steel rod entered the cryostat through a ball-and-socket joint and sylphon 
bellows. The cryostat was mounted on a General Electric XRD-—5 
spectro-goniometer. 

A copper-constantan thermocouple and a carbon resistance thermo- 
meter were attached to the specimen. Temperatures other than liquid 
helium and nitrogen were attained by allowing the helium to evaporate 
and admitting air to the helium chambers. The thermal inertia of the 
system was such that this provided sufficient temperature control for 
this series of experiments. 


2.3. Cold Work 


The usual method of determining « is to cold-work the material by 
filing and to follow this by an x-ray examination. In the present method 
cold work was introduced into bulk specimens by scratching the surface 
with a pointed tool (as shown in the plate). This was not as satisfactory 
as could be desired; however, with proper care, the main part of the 
exposed surface could be severely worked at any temperature. 

In order to compare the deformation which could be obtained by this 
method, an annealed bulk specimen of 70:30 a-brass was deformed at 
4-2°K and the faulting frequency compared to that reported by Wagner 
(1957 a) who deformed his specimen by filing at 77°K. 


2.4. X-ray Method 


In the examination of bulk specimens, peak positions for the (111) 
and (200) reflections, prior to and after cold-working, were determined 
as the mid-point at half-peak height. However, for the (220) reflection 
in the unworked specimen a separation between Ka, and Ka, peaks was 
apparent. Therefore, as a rough measure, the mid-point at quarter-peak 
height was used to determine the peak position for all (220) lines. 

The foregoing approach allowed a precision in the determination of 
peak separation of about 0-01° for the (200)-(111) pair and somewhat 
less for the (220)-(200) pair. The determination of peak separation for 
higher-order reflections was not attempted. 

In a few cases the counting rate constant was changed between two 
peaks in order to give equal apparent line intensities. No measurable 
change in peak separation with change in counting rate was observed. 

Copper Ka radiation (A= 1-5418 4) was used and the scanning rate was 
0-20°/min (26). 

§ 3. RESULTS 
3.1. Pure Lead Deformed at 4:2°x 


The bulk sample of pure lead was first deformed at 77°K in the 
cryostat then warmed up to between 160° and 170°K for one hour. 
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Reerystallization occurred and was indicated by the full recovery of 
the line broadening and peak intensities. (Barrett, private communica- 
tion) has observed recrystallization at ~155°K in lead of comparable 
purity.) In this way, a fine grain bulk sample (‘recrystallized’) was 
obtained for subsequent deformation at 4:2°K. It can be seen from 
the data given in table 1 that the peak separation changed with cold-work. 
As expected, the (200)-(111) separation decreased and the (220)-(200) 
increased with deformation. Moreover, it is demonstrated that the change 
in separation increased with increasing cold-work. For comparison, a 
bulk sample of 70:30 «-brass was given the same treatment as the lead, 
producing a change A20o99)431)=9°17°. 


3.2. Lead Alloys Deformed at 4:2°K 


In table 2 the data for the alloy with 0-10at.°% Ag show that there 
was little difference in A20, produced by cold-work, from that of the 
pure metal in table 1. Since the alloy had a fine initial grain size, no 
recrystallization treatment was given prior to measurement. On 
warming the 2659)-20,,, peak separation recovered almost to the initial 
room-temperature value; whereas, the line intensities were only 20% 
recovered at this temperature. Recrystallization does not occur until 
higher temperatures are reached in this alloy (Bolling and Winegard 
(1958 b)). 

The alloy with indium did not exhibit measurable faulting after cold 
work (table 2). The change in separation is on the border-line between 
a true change and a random measurement error. Therefore, the effect 
of indium has been to reduce or to eliminate the separation change. 


3.3. Pure Lead Deformed at 77°K 


In table 3 a sequence of measurements is presented for one sample 
that was deformed and examined both at 77°K and 4:2°x. The 
recrystallized state noted in the table was obtained after deformation as 
before, when the sample was warmed from 77°K to 172°K at a rate 
between 1 and 2°K/min. During this heating the lead recrystallized and 
the peak-separation change was recovered. It was noted from traces 
obtained at intermediate temperatures that the line broadening recovered 
at lower temperatures (viz.: ~135°K) than the peak-separation change 
(~170°K). An estimate of recovery may be made from changes in line 
broadening; however, only a crude estimate could be obtained. This 
is shown in the table by the half-width at half-intensity of the (111) 
trace. ; 

The peak-separation change at 77°K is comparable to that at 4-2°% 
for heavy deformation as in table 1. 


3.4. Aluminium Deformed at 4:2°K 


The aluminium specimen was deformed and annealed to produce a 
uniform fine grain size. Peak separation changes were only obtained at 
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4:2°K and are shown in table 4. As can be seen, the change in peak 
Separation is just measurable after heavy cold-work. Wagner (1957b) 
has already shown that measurements on aluminium filings prepared and 
measured at 77°K reveal no change in peak separation with cold-work. 


Table 3. Zone-refined Pb 


X-ray exam. 
temperature 


Half-width (111) 
at 


Change in 
A26, with 
cold-work 


Treatment 31, 


26 


From room 
temperature 
Recrystallized 
Ist cold-work 
2nd cold-work 
3rd_cold-work 
Recrystallized 
Recrystallized 
Ist cold-work 
Recrystallized 
Recrystallized 


Table 4. 


Aluminium at 4:2°K 


Change in A26, with 


Treatment (20 ap9-2F 44) = A208, oP oe 
Annealed 6-27; —~ 
lst cold-work 6-27, — 0-00, 
2nd cold-work 6-26, —0-01; 
3rd cold-work 6-265 —0-01, 


3.5. Stacking-fault Parameter 
Following the method described by Wagner (1957a), « may be 
calculated as a. H=A20,,,_,,,,, where H is a constant depending upon the 
particular reflections and radiation used. For CuK« radiation the 
following values were determined for H : 


Pb 200-111 H=—3-7 
Pb 220-200 H=+7-8 
70:30 Brass 200-111 H=—5-0 
Al 200-111 H=-—4-7 


In table 5 the value of « for the heaviest deformation used in this 
investigation is tabulated together with some other values taken from 
the literature. The values of « for pure lead calculated from A269 999 
are not included in table 5 because the determination of peak positions 
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of the (220) are less accurate. However, it can be seen from table 1 
that the sign of the change in separation is consistent in giving a positive 
a. (For first and second cold working, « is respectively 0-003 and 0-005 
as compared to the A26,..,,, values 0-008 and 0-013.) 

In table 6 the highest values of stacking-fault parameter determined 
at any temperature are compared with twin boundary energies. These 
are determined from relative energies (see, for example, Bolling and 
Winegard (1958a)) and the best values of grain-boundary energy. 
When the latter was not available, it was chosen as one-third of the 
surface-free energy. 


Table 6. «,,,, compared to twin-boundary free energy 


Metal yr.p. ergs/em? Method 


Highest 
a 


Twinning more frequent than 
Cu-—Zn (30) 040 20 Cu, grain-boundary energy 

probably lower. 

Estimated relative energy as 
~ 0-04 

Measured relative energy 

Measured relative energy 

Measured relative energy 

Measured relative energy 


§ 4. Discussion 
4.1. Correlation between a, y and Twin-boundary Energies 


In table 6 a correspondence between twin-boundary energies and the 
stacking-fault parameter is apparent. It was also established that the 
temperature of deformation and the amount of cold-work introduced to 
the metal are important in determining the value of «. The temperature 
of measurement can also effect the value of a if any recovery occurs. 
If all these factors could be taken into account, a comparison of different 
techniques used to assess y would have to include some type of comparative 
temperature scale for the metals of different melting points. It is 
possible that results obtained at two extremes of temperature might 
be suitably compared. For example, the relative twin-boundary energies 
for copper and lead were determined after long annealing near their 
respective melting points where the structures should be in a similar 
condition. On the other hand, it would probably be desirable to compare 
the values of « at the lowest possible temperature and for equally strained 
specimens. Nevertheless, the general trend in results confirms that lead 
should have a low numerical value of stacking-fault energy, which is 
comparable to copper in the toripergture range 0-160°K where extensive 
cold-work can be obtained. 
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The data for lead alloys again show a general correlation between 
and twin-boundary energies, but emphasize the difficulty of any direct 
correspondence. As shown in table 5, « decreases from 0-013 to 0-011 
with addition of 0-10at. % Ag to lead. This may not be a true change 
because the degree of cold-work may be smaller in the alloy which is 
harder than the pure metal. Bolling and Winegard (1958a) showed 
that the relative twin-boundary energy changed from 0-:050+ 0-014 to 
0-077+0-016 with the addition of 0-10at.% Ag to lead. However, 
they attributed this change to a decrease in grain-boundary free energy 
in the alloy caused by the segregation of silver. If this conclusion were 
true, the lead-silver alloy is probably very much like pure lead with 
respect to the stacking-fault energy. 

The addition of 20-Oat.°% In to lead suppresses a measurable a as 
shown in table 5. Subsequent to the x-ray studies, samples of the 
same alloy were treated and examined for relative twin-boundary energy 
and twinning frequency in the manner used by Bolling and Winegard 
(1958a). Even after long annealing close to the melting point scatter in 
the data prevented an accurate determination of relative twin-boundary 
energies. However, on the average, the values obtained were higher 
than could be attributed solely to grain-boundary segregation of indium 
and a subsequent reduction in the free energy of grain boundaries. 
Moreover, the ratio of the number of twin boundaries to grain boundaries 
was reduced by a factor of six in the indium-—lead alloy from that in pure 
lead. These results can be interpreted as indicative of an increase in 
twin-boundary energy, and hence they correlate with the decrease of « 
in this alloy. 


4.2. Electronic Interactions and Stacking-fault Energy 


The stacking-fault energy of a metal arises from particular electronic 
interactions within the structure. However, the way in which the 
influence of the electronic structure upon y may be assessed is not clear. 
Aluminium is three-valent, and its first, and second, Brillouin zones are 
overlapped by a complex Fermi surface (Heine 1957), unlike the mono- 
valent metals, copper, silver and gold. Hence, if the Fermi electrons 
and the nature of the Fermi surface are important in determining y, 
one might expect aluminium to differ from the noble metals which do 
not exhibit overlap. For similar reasons all polyvalent f.c.c. metals 
might be expected to resemble aluminium (Seeger 1955 b). 

There is evidence that the first Brillouin zone, and higher zones of lead, 
are also overlapped by a multiply-connected Fermi surface (Gold 1958), 
Nevertheless, it is significant that the stacking-fault probability of lead 
falls with the noble metals rather than with Al. If the long-range 
electronic forces are responsible for the stacking-fault energy, the 
difference between Pb and Al may originate in the detailed difference 
between their respective Fermi distributions. This is suggested by 
reference to the Pb—In alloy. 


X-ray Study of Deformation Stacking Faults 501 


For the purpose of assessing an alloy of lead, a group IV metal, with 
indium, a group III metal, it may be justifiable to assume that the 
electron concentration decreases on adding indium to lead. The eleetronic 
nature of lead—indium alloys has not been. studied; however, inter- 
pretations exist for the lead—thallium system which shows a very close 
resemblance to lead—indium. In the lead—thallium alloys the changes of 
various electronic properties indicate that the density of states, N(Z), 
at the Fermi surface drops rapidly from the value of pure lead, reaching 
some fairly constant value at about 15% thallium (Gold 1960). This 
decrease of NV(#) has been interpreted in terms of changes in the Fermi 
distribution as the nominal electron concentration decreases. Similarity 
between the systems lead—thallium and lead-indium suggests that an 
analogous change in N(#) may also occur in lead-indium alloys, and 
hence it is tempting to speculate about the connections between N(£) 
and a. Since addition of 20at.% In to lead produces a change in a 
from a copper-like metal to an aluminium-like metal, the above possibility 
deserves further detailed study. 


4.3. Summary 


1. At low temperatures the value of « for lead is similar to the value 
of « for the monovalent f.c.c. metals and about one-third of the value 
for 70:30 CuZn brass at 4-2°K. 

2. The value of « for aluminium at 4-2°K is barely measurable, indicating 
that even at this temperature aluminium does not exhibit appreciable 
deformation faulting. | 

3. A finite stacking-fault probability in lead at low temperatures can 
be suppressed by the addition of 20 at.% In. 

4. A correlation between the presence of annealing twins and the value 
of « is exhibited in the lead and lead alloys, even though it is clear that too 
many extraneous factors are present during the production of annealing 
twins for more than a qualitative comparison. 

5. It has been demonstrated that the stacking-falt probability increases 
with the amount of deformation in the systems studied. Hence, it is 
difficult to compare a for different metals if the amount of strain is not 
assessed. 
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ABSTRACT 


Theories are developed from two different assumptions to calculate the 
effect of irradiation growth on the creep of uranium under a uniaxial applied 
stress when the temperature is sufficiently high for the internal stresses set 
up by the differential expansions of the individual crystals to be insufficient 
of themselves to cause plastic flow. It is shown that the greater the irradiation 
growth the greater is the creep rate for a given stress ; the effect on creep rate 
is also greater at lower applied stresses for a given irradiation growth. 


§ 1. LyrropuctTion 


Anderson and Bishop (1957) have developed a theory for the effect of 
irradiation growth on the creep of uranium when the internal stress 
generated is sufficient for the material to become plastic. The theory 
presented here is similar to that derived (Blackburn et al. 1960) for thermal 
cycling effects and treats the case when the stresses are still in the elastic 
range. For simplicity, primary creep, variations with stress of Young’s 
modulus and variations of the mechanical properties with time of irradia- 
tion are neglected, though these effects may be incorporated in the 
equations. 

We obtain simultaneous differential equations, from which the creep 
rate may be determined by integration, from each of the following 
assumptions: (a) a uniform deformation; (b) a uniform internal stress 
which is the same as that generated if each crystal were prevented from 
deforming in the absence of the applied stress. 


§ 2. THEORY 


To calculate the deformation of an aggregate of non-isotropic crystals 
from their individual properties, assumptions must be made, since the 
orientation and location of each crystal is not known. For simplicity 
we assume each crystal to undergo a homogeneous deformation. In 
general this prohibits the continuity of both stresses and displacements 
between the crystals. We require however on quasi-statical grounds 
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that the mean of any component of the stress over the crystal orientations 
should equal the corresponding component of applied stress. Another 
assumption is also necessary. ‘The simplest is that the stress is uniform 
(i.e. the same in all crystals). This however leads to the conclusions that 
the creep rate is the same under temperature variations as it would 
be if the temperature were maintained as its instantaneous value, in 
contradiction to the experimental evidence (Anderson and Bishop 1957, 
Blackburn et al. 1960). 

Two other possibilities are to assume either that the deformation is 
uniform or that the stress components are the sum of those of the applied 
stress w and an internal stress which is uniform with respect to the crystal 
axes and independent of w. The first is probably more realistic but the 
second leads to simpler computation and can be pursued further. We will 
develop both. 

In each case we assume that the strain rate of a crystal under homo- 
geneous deformation is related to the stress, stress rate and rate of 
irradiation growth by 


de; l+vd v 3\ (n+/2 ay, 
=F = Vig t a 77 Pere = 271) - (5) Cexp (— B/T)o,J"1. 
| (1) 


The first term on the right-hand side represents irradiation growth, the 
second is the derivative of the linear elastic strain stress relationship and . 
the third is a generalization of the form of the one-dimensional secondary 
creep law for isothermal polycrystalline uranium as found experimentally. 
T is the temperature, ¢ the time, # Young’s modulus, v Poisson’s ratio, 
_ B, C and n creep constants ¢;;, 9j;, 7;;, 04; the components of the strain, 
metric, stress and stress deviation 


(or reduced stress 7;;—1/39,;>'7;, ‘) tensors and J*= Sio,,0;;. 
k i] 


We consider a uniaxial stress w acting along the z axis of a set of 
coordinate axes in the specimen. Without loss of generality we may 
choose the 6 axis of the crystal to lie in the (x, z) plane. We denote 
the angles between the 6 and z axes and between the (b, z) and (b, c) 
planes by 47-06 and 47-1. Then 


Y11 = y (cos? 6 — cos? / sin? 4), Yo3 = y Cosy sin cos 8, 
Viz = y C08 O sin 6(1 + cos? x), Yoo= —ysin? y, (2) 
93 = y(sin? 0 — cos? ys cos® 4), Yio= —y cos sin y sin 8. 

By adding the equations for the normal strain rates we find that, 


d 1—2vd 
al DKK = “E di 2k 


For a random orientation the mean of this is zero (since the mean of 
DT ee 
k 


is the applied stress w). Hence the creep rate de/dt is the mean rate of — 
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normal strain deviation along the z axis. We note that the mean rates 
of strain deviation in the perpendicular directions are — 4 of this and the 
mean tangential strains are zero. 

We consider first the assumption that the deformation is uniform. 
Then eqns. (1) and (2) become: 

doz, EH {"- 


isa ToES di (3)ern C exp ( — B/T )o3,J°-* 


+ y(cos? % cos? @ — sin? @) 


d E 
aah Poe mis — (3)? Cexp (—B/T)o, J" 


dt 12% | 7a 
+ (cos? us sin? 8 — cos? a} : 
do43 EK 


=a 13 {- (2)e+) 12 O exp ( — B/T)o,,J5° —y(1+ cos? ) cos @ sin 6, 


—= { — (3)2?C exp (— B/T)o.,J"— — y cos # sin cos 6}, 


{— (3)@*)2C exp (— B/T)o4,J" 1+ y cos sin sin 6}. 


(3) 
We remember that o,, may be obtained as — (a4; +053). 
If we denote the properties of a crystal whose orientations, 6 and % 
are 6,, and %,) by bracketed subscripts and the frequency of its 
distribution by Si, ;), we have the additional equation: 


tig 
DS, n)733(k, h) = 5D Six, h)* 
kh kh 


Hence for random orientation it may be deduced from the first of 
eqns. (3) that 


rr DBShe as = ($)"+ D2 C exp ( — B/T)> Se, nok, h)ogax, nh) 
kh hh (4) 


Integration of the simultaneous differential eqns. (3) and (4) for a small 
number of orientations yields the stresses and the creep rate. 

We now investigate the consequences of the alternative assumption 
that the stress in a randomly orientated aggregate is the sum of the 
applied stress and a uniform internal stress. The value of the latter 
may be determined by putting 0=7/2 and % and de/dt zero in eqn. (3). 
We denote the value of o;; when w is zero by 8. 

Seer Sap» Saco and ,, are initially zero. We see that eqn. (3) is satisfied 
if they remain so. Also s,,= —8,,_ where: 


a as {y= (BPC exp (—BIP)Bal"3} + + (8) 
and | . P= DSisSii- 


9 


P.M. 2k 
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Hence we may write ace — 8,4, =1//2 where 
aE Fe yay (gym Cexp(—B/TIL"}.- (8) 
dt l+vp 


By integrating this equation, J may be determined, but it is often possible 
to approximate it by its asymptotic value: 


Table 1. Time for internal stress to approach asymptotic value at 403 e 
in a flux of 3 x 1012 neutrons cm~? sec"? 


Table 2. Asymptotic values of internal stress 


270 0-18 549000 
305 0-14 162000 


380 0-068 16600 
400 0-046 8980 
403 0:04 7950 


1750 


The times taken to approach this value at 403°c in a flux of 
3x10" neutronscm™sec™' are shown in table 1, using the values 
n= 2:52, v= 0-23, y=0-04/yr, B=25300°K, H=1-86x 10’ p.s.i., C=8-413 
hr p.s.i.-?°*. With these constants, the asymptotic values of J at the 
temperatures for which y is given by Finniston (1958) are presented in 
table 2. At the lower temperatures these stresses will not be attained 
since the material will previously have yielded plastically. 

We now use the assumption that the stress is the sum of the applied 
and internal stresses. The mean creep rate over the various orientation 
simplifies from eqn. (1) to 


d 1 . 
= = — ()"+D2 CO exp (— “BI |T) ie ie o33 J"! cos 0.d0 dibs. "tp 


—n/2 
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Then, as shown by Blackburn et al. (1960) 

O33 = 833+ 20/3, J? =I? + 2wsy, + 2w?/3, 
and the creep rate reduces to: 

de 

dt 
when w/J is small, and in general to: 


d il 2w 2 (1 — v2) — a2 
= a (Bm mPexp ( _ B/T) or {" I. (H+ cos a S| 


‘2w[ cos? b(1—2x?)—ax?] — 2w?) bre 
i} A Ped Bess bl Sa es a wee 
«{ + i + = aedy i 2 (9) 


Thus the ratio of the creep rates in the presence and absence of 
irradiation growth is: 


SONA Ceara 


_ /2w[cos? b(1—x?)—a2?] — 2w?) be 
x 4 L- S++ __ : 1 
{143 =3 T + 3 o dx dis (10) 

With n=2-52 the constant appropriate to uranium this expression and 
the linearized approximation, 

(g)mt0/2 2(n + 4) a iad 

15 Ww 
are presented in table 3 as function of w/I., The ratios of creep rates 
presented differ little from those which would obtain if the internal stress 
were caused by thermal cycling (Blackburn et al. 1960). 


= (gyre Cexp(— Biol, . (8) 


Table 3. Effect of irradiation on creep rate 


Creep rate with irradiation Teer ed 


approximation 


2 
= 


Creep rate without irradiation 


if 


nr 


5 


BE RoOIISa 
St 


Oe OTB OO 


SCOOrFNnNas 
wy) 


0-2 
0-5 
0-7 
1-0 
1-5 
2-5 
5:0 


It has been found experimentally (Rose 1957) that the secondary creep 
rate of uranium at 450°o under a stress of 2-7 t.s.i. does not differ markedly, 
with and without irradiation. This would be anticipated on this theory 
for the conditions reported but would not be true for small stress at this 
temperature. 
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§ 3. CONCLUSIONS 


Theories have been derived to deduce the effect of irradiation growth 
on the creep of uranium under a uniaxial load when the stress is insufficient 
to cause plastic flow. The theory based on the assumption of uniform 
deformation is probably slightly more reliable, but the alternative 
assumption, of the stress being the sum of the applied stress and the 
internal stress that would exist if each crystal were prevented from 
deforming in the absence of applied stress, leads to a simpler formula 
for the creep rate. For w/J less than 1/3, the proportionate increase in 
creep rate is seen to vary approximately as ([/w)"—', while J tends to a 
value which varies in proportion to the nth root of the rate of irradiation 
growth. 
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ABSTRACT 


The optical constants of liquid germanium, tin and lead have been measured 
by a reflection method for wave numbers between 4000 and 27 000 em-! 
(wavelengths 2-5 to 0:37 microns). The temperature variation of the optical 
constants was measured for tin and lead. The experimental results follow 
approximately the Drude free electron formulae if the number of free electrons 
per atom, N,, and their relaxation time, 7, are treated as adjustable 
parameters. The values of N, lie between 4:3 and 4:7, with a slight 
temperature variation. The values of the static conductivity calculated 
from N, and 7 are compared with electrically measured values. Previous 
optical measurements on evaporated films of tin and lead have indicated 
values of N, between 1-2 and 1-4 for the solid metals. 


§ 1. [InrRopucTION 


THE electronic band structure of a solid metal is determined by its crystal 
lattice structure. When the metal melts, the long-range order of the 
crystal lattice disappears but some short-range order remains. Knight 
et al. (1959) have suggested that this short-range order determines an 
electronic band structure for a liquid metal. For those metals in which the 
coordination does not increase markedly on melting, e.g. Al, Sn, Pb, the 
electronic structure of the liquid is similar to that of the solid. They quote 
experimental values of the changes of electrical conductivity and Knight. 
shift in support of their suggestion. For example, the electrical conducti- 
vity of tin decreases by a factor 2-1 on melting and the Knight shift. 
decreases by a factor 1:03. The change in the conductivity is in 
approximate agreement with the theory of Mott (1934), assuming no 
change in the number of free electrons when a metal melts. The small 
change in the Knight shift indicates a correspondingly small change in 
the electronic wave function. The behaviour of germanium and other 
~metals for which the coordination increases on melting, is quite different. 
and drastic changes in the electronic wave function must be assumed. 
The changes in the properties of germanium on melting are exceptionally 
large, with a ratio of electrical conductivity liquid/solid of about 15 at the 
melting point. Liquid germanium has a resistivity of 6-0 x 10-° ohm-cm 
at the melting point, which is close to the value for liquid tin at the same 


temperature. 
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Knight et al. also suggest that there is no simple connection between the 
optical properties of a liquid metal and its transport properties and Knight 
shift. The only values for the optical constants of liquid metals previously 
available, those of Kent (1919), are at variance with this suggestion. The 
electron densities and relaxation times deduced from Kent's values for the 
optical constants, lead to values for the electrical conductivity nearly equal 
to the electrically measuted values. Since Kent made measurements over 
a small spectrum range and at only one temperature, this agreement might 
be considered a coincidence. To find out if this is so, the optical measure- 
ments on liquid metals are being repeated over a wider spectrum range, wave 
numbers 4000 to 27000cm—, and temperature variations are being 
measured where possible. “This paper presents the results for three quadri- 
valent elements, germanium, tin and lead. 


§ 2. EXPERIMENTAL METHOD 


The optical constants were measured by the reflection method described 
in previous papers (Hodgson 1959, 1960). The metals were melted in 
silica or carbon crucibles enclosed in a tank filled with hydrogen. The 
temperature just below the centre of the crucible was measured with a 
chromel—alumel thermocouple. The experimental results are expressed in 
values of the dielectric constant, «, and the conductivity, o, in e.s.u. for 
various wave numbers, k,inem—. 


§ 3. THEORY 


A metal contains a certain density of conduction electrons, say V per cm’. 
The motion of these electrons is damped by the thermal vibrations and 
lattice imperfections. The simplest representation of this damping is by 
a single parameter, the relaxation time, 7. Interaction with the atoms of 
the metal may also change the effective mass of the conduction electrons 
from the free electron value, m, to m*. The theoretical formula for the 
static electrical conductivity, a9, due to these electrons is: 


oye Orly 4 Wels drepsheete tes 


At a non-zero frequency corresponding to a vacuum wave number k, the 
formulae for « and o due to the conduction electrons are : 


e=1—k,?/(k? + k,?), Pare Phas a 
o/C= $kyky/(k? + kp?) MEARS yer grt tt 
where 
k=eN/m*ne* and kp=1/2mer. poo Piet Oana Ue 
The effective number of electrons per atom N4, is defined by : 
No=(N/Na) x (m/m*) 10 tal yeh Wty 


where N, is the number of atoms’ per cm. At frequencies below the 
threshold for inter-band transitions, the effect of electrons in the filled 
bands can be represented by a small additional term in the expression 
for «. This term has been neglected in the discussion of results. . 
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§ 4. Discussion or RESULTS 
4.1. Germanium 


The germanium used for these measurements had a room-temperature 
resistivity of about lohm-cem corresponding to an impurity content of 
between 10% and 10!%atomspercm?. About 40g of germanium were 
melted in a silica crucible under a pressure of about 3cm Hg of hydrogen. 
The hydrogen was constantly flowing through the tank containing the 
furnace. On first melting, the liquid germanium was covered with a 
cloudy, light scattering layer. This layer gradually cleared from the 
middle of the liquid germanium surface. The specularly reflecting area 
in the middle was sufficient for reflection measurements. The optical 
constants were measured and the germanium allowed to re-solidify. It 
was subsequently re-melted and the optical constants were measured again. 
The values were appreciably different from those of the first experiment, 
and o having changed by ~20%. This process was repeated about six 
times until the change of optical properties on solidification and re-melting 
was not more than the experimental error. These changes may be 
explained by the gradual reduction of an oxide film on the surface of the 
germanium. It was not possible to make reliable measurements of the 
temperature variation of optical properties because of the small temperature 
range available, 980 to 1050°c. 

The final experimental values of « and o/c are given by the points in 
fig.1. The estimated probable errors are shown by vertical bars where they 
exceed the diameter of the plotted circles.. The theoretical curves have 
been drawn using formulae (2) and (3), choosing ky and ky to give the best 
fit with the experimental points for log k< 4-2, i.e. in the infra-red. The 
general trend of the experimental points agrees with the theoretical curves 
but there are significant differences in detail. 


4,2. Tin 

Optical measurements were made on’two samples of tin of different 
purities, 99-999 and 99-99%. The tin was melted in a silica crucible under a 
pressure of about 5em Eg of hydrogen. As with germanium, the optical 
constants changed on successive melting and solidification until stable 
values were found. The final values for both samples of tin were not 
appreciably different. ‘Tin is a convenient metal for the measurement of 
temperature variation because of the low melting point, 232°c, and the 
high boiling point, ~ 2000°c. 

The final experimental values are given by the points in fig. 2. The 
theoretical curves have again been fitted to the experimental points for 
logk<4-2. Only a few points show significant deviations from the 
curves. ; 

Some experimerital values of « and o for evaporated films of tin are 
available for comparison (Motulevich and Schubin 1957 and Hodgson 


Liquid germanium. 
c=dielectric constant, o=conductivity in e.s.u., k=wave number in cm™?. 
Experimental points for 998°c. Theoretical curves for N )=4-30, 
o/c =40°7 x 104. 
Fig. 2 


Liquid tin, 


e=dielectric constant, c=conductivity in e.s.u., k=wave number in cm-, 
Experimental points : © for 858°o, V for 446°c. Theoretical curves : 
A for Njy=4:59, o)/¢e=43:8x104; B for Ny=4-68, o/e=52-4 x 104. 
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1955). These measurements on evaporated films are not exactly repro- 
ducible and further experimental work is needed before such values can 
be accepted as standard values for solid tin. They should be accurate 
enough, however, to show some general features of the optical properties 
of the solid metal. The experimental values are shown by the points in 
fig. 3 and the theoretical curves have been fitted to the points for 
logk<3-6. The experimental points for o lie consistently above the curve 
for logk> 3-6. The general agreement of the experimental points with the 
curves for log k <3-6 indicates that the corresponding values of N, and + 


Fig. 3 


3-0 35 4:0 
Solid tin. 


e=dielectric constant, o=conductivity in e.s.u., k=wave number in cm™}, 
Experimental points : V Motulevich and Schubin (1957), O Hodgson 
(1955). Theoretical curves for Vy=1-33, o)/¢e=1-53 x 10°. 


are significant. The optical value of o,/c is smaller than the electrical 
value, but a similar discrepancy has been noticed in measurements on 
evaporated films of other metals. For example, infra-red measurements on 
a silver film (Hodgson 1955), gave Ny=1-08 and o,/c (optical) = 0-38 o/c 
(electrical). The value of N, for the silver film agrees with the value 
suggested by other experiments and by theory. The value of NV 0 deduced 
from the infra-red measurements on a tin film should therefore be significant 


for the solid metal. 


pl4 . J. N. Hodgson on. the 
4.3. Lead 


Lead'of 99-99%, purity was melted in a carbon crucible. under a pressure 
of about 3em Hg of hydrogen. | The optical properties of liquid lead did not 
change appreciably on repeated solidification and melting. The lead could 
be heated to 800°c without appreciable evaporation on to the tank windows: 
The optical constants were measured at two temperatures about 200°C 
apart. The experimental! values for the higher temperature are given by 
the points in fig. 4. The theoretical curves were fitted 'to the points for 
logk< 4:2, but all the points lie close to the curves. ‘The agreement 
between experimental] points and theoretical curves is similar at the lower 
temperature but these have been omitted from fig. 4 for the sake of clarity. 


Fig. 4 


Liquid lead. 


¢«=dielectric constant, o=conductivity in e.s.u., k=wave number in cm—, 


Experimental points for 786°c. Theoretical curves for WN. o= 4:57, 
oy/c = 24-2 x 104. 


The optical constants of an evaporated film of lead have been measured 
by Motulevich and Schubin. Their results have been analysed to give 
approximate values of NV, for solid lead. ; 


§ 5. CONCLUSIONS 


The values of Ng, 7 and o,/¢ (optical) corresponding to the curves in 
figs. 1 to 4, are set out in the table. The latitude in fitting the curves to the 
experimental values leads to an uncertainty of about + 1% in the tabulated 
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values. For all the liquid metals, the effective number of electrons per 
atom, JV, is greater than 4-0, the number of valence electrons. In other 
words, the effective mass is apparently less than the electronic mass. The 
decrease of Ny for liquid tin and lead when the temperature increases is just 
large enough to be significant. N, for liquid tin and lead is larger by a 
factor of between three and four than N, for the solid metals. This is in 
contrast to the small difference between the Knight shifts for solid and 
liquid tin. 


o/c Gy /C 
poral : Ne i (optical) (plectrical) 
Ge (1) 998 4-30 2°47 40-7 49-5 
Sn (1) 446 4-68 3-90 52-4 56-1 
Sn (1) 858 4-59 3°52 43-8 47-2 
Sn (s) 20 1-33 37-4 153 276 
Pb (1) 567 4-68 2°30 26:8 28-4 
Pb (1) 786 4:57 2:18 24-2 25:9 
Pb (s) 20 1-26 29-9 103 145 

8) x 10-71% sec x 104 x 104 


(1) liquid, (s)-solid, ¢ temperature, NV, effective number of electrons per atom, 
7 relaxation time, oy static conductivity in e.s.u. 


The optical values of o,/c for the liquid metals are all less than the electric- 
ally measured values. The difference decreases in the sequence 
germanium, tin, lead. A similar but larger discrepancy is noted for the 
evaporated films of tin and lead. These discrepancies may be due to the 
fact that the surface layer of a metal does not have the same properties as the 
bulk metal. The surface may cause extra damping of the electronic 
motions above that found in the bulk metal. The temperature variations 
of o,/e (optical) for liquid tin and lead are closely correlated with the 
changes in o,/c (electrical). This lends support to the original view of Kent 
(1919) that in these liquid metals the same valence electrons are responsible 
for both the optical properties and electrical conduction. 
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ABSTRACT 


Alpha-particles of the primary cosmic radiation have been studied in a 
stack of nuclear emulsions exposed over northern England on July 29, 1959. 
A flux of 167+ 12 «-particles/m? ster see was found. The energy spectrum 
was examined between the cut-off energy of about 250 mev per nucleon and 
an energy of 1-5 Bev per nucleon. It was not found to be significantly 
different from that observed during solar maximum. The energy spectrum 
was also examined as a function of zenith angle. An apparently significant 
linear relationship was established between neutron monitor counts recorded 
at sea level and primary «-particle flux values. <A value of 24:0+ 2-4 em 
has been determined for the mean free path of «-particles in nuclear emulsions, 
which is somewhat higher than previously reported values. 


§ 1. InTRODUCTION 


In the last few years there have been a number of experiments to 
investigate the behaviour of the a-particle component of the primary 
cosmic radiation. In particular, much effort has been directed towards 
studying the variations in the flux and energy spectrum that have 
occurred during the solar cycle. These measurements have been 
summarized in a recent review article (Waddington 1960a). Of special 
interest is the distribution of the lower energy particles, because of the 
particular sensitivity of these particles to any alterations in the conditions 
of the inter-planetary medium they traverse. In order to investigate 
these slow particles it is necessary to make exposures at latitudes where 
the geomagnetic cut-off energy is as low as possible. Due to the position 
of the magnetic pole, this means that, over Europe, exposures should 
ideally be made at geographic latitudes higher than about 60° N, compared 
with the 45° N that is adequate over northern America. For this reason 
a series of flights was made in the summer of 1959 over northern England, 
with the intention of recovery in Northern Ireland; the results obtained 
from the emulsions carried on one of these flights are considered in this 
paper. 

These flights were planned to take place at a time when the sun’s 
activity was appreciably diminished from that at solar maximum, 
1957-58, since it was hoped that it would be possible to study the process 
of recovery from the low flux values recorded during those years. An 
observation of the energy spectrum during such a period of recovery 
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would be of particular interest. It has been pointed out by McDonald, 
and Webber (1959) that the change in the spectrum from solar minimum 
to solar maximum was in good agreement with that predicted by the 
electric de-acceleration theory of Nagashima (1953),.in which each particle. 
passes through a constant retarding field; and does not accord with any 
of the qualitative results of those theories that invoke magnetic modulation, 
in one form or another. While not entirely conclusive, an observation 
that the higher energy particles recover before those of lower energy would. 
be strong evidence in favour of a magnetic modulation theory, while. 
equal recovery of particles of all energies would favour the electric 
modulation theory. 

Unfortunately, these flights chanced to occur during a temporary, but. 
very active, period of renewed solar activity which produced intense. 
geophysical activity, and the consequent failure to observe any significant 
difference from the 1957-58 data may have been a temporary consequence, 
or may mean that no recovery had yet occurred. If the apparent 
relationship established later in this paper between sea level neutron 
monitor counts and a-particle fluxes is valid at all times in the solar 
cycle, then no significant difference could have been expected. 


§ 2. Exposure DETAILS 


A stack of 100 Ilford G5, 600 stripped emulsions, measuring 
20x 15cm, was exposed on July 29, 1959. The balloon carrying 
the stack was launched from Duns, Berwickshire (55°47’N, 2°20’ W). 
and the equipment fell at Omagh, Co. Tyrone, Northern Ireland 
(54°36’N, 7°18’ W). The flight curve is shown in fig. 1. While the 
balloon was at ceiling it drifted almost due west, so that the flight was. 
assumed to be at a constant geomagnetic latitude of 57:-7°N. The 
cut-off rigidity predicted by Quenby and Webber (1959) for this location 
is 1-65Bv, corresponding to a cut-off energy of 310Mev per nucleon. 
for w-particles. This is not as low as might be desired, but is as far north 
in the British Isles as is consistent, under normal meteorological conditions, 
with making a balloon flight of useful duration and still recovering the. 
equipment. 

§ 3. EXPERIMENTAL PROCEDURE 

Twelve emulsions near the centre of the stack were each scanned for. 
10cm along a line 2cm below the top edge. The scan-line had to be 
placed so far below the top edge because of the serious curvature of the 
emulsions in the top centimetre. Tracks crossing the scan line were- 
recorded if: . 


(i) They had a projected length of greater than lem. 


(ii) They had an ionization greater than about three times the minimum. 
value. . a 


(iii) They had a zenith angle, , less than or equal to 60° 
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In recent years, it has been customary to use smaller values of 9 in 
similar experiments, in order to reduce the uncertainties of the corrections 
to the top of the atmosphere and to minimize the geomagnetic effects 
that occur at large zenith angles. How serious these geomagnetic effects 
are has not been considered in detail. For this reason, the data obtained 
in this experiment have been separated into two groups of particles, 
those with $< 30°, and those with 30°<9<60°. These two groups have 
been examined for any significant differences. 
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The time-altitude curve of the balloon used to expose the emulsions. 


The efficiency of detection of «-particles achieved by the scanners in 
this experiment was checked by applying the usual tests. 

(i) Re-scanning. The scanning of the emulsions was done at two 
different times, and so the data have been divided into two sets (a) 
-and (b). There were three scanners, and their efficiencies, as determined 
from ré-scanning, were 98%, 90% and 88%. In set (a), each plate was 
rescanned, and the overall scanning efficiency was 99°. Only one plate 
was re-scanned in set (b), and so the average scanning efficiency was 
only 94%. 

These results indicate that there was a small random scanning loss, 
but the determination of efficiencies by cross-scanning does not exclude 
the possible existence of serious systematic losses. _ Some workers appear 
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to have assumed recently that the demonstration of a high efficiency in 
re-scanning is sufficient. to prove a high overall detection efficiency. It 
must be pointed out that such an assumption implicitly assumes that 
all particles have an equal intrinsic probability of being detected. 
Experimentally, this assumption is probably not true. It appears 
reasonable to expect that particles lying near the glass interface, and 
those with a projected length near the minimum accepted, may be more 
difficult to detect than those more favourably situated. For this reason 
additional tests have been applied to the data. 

(ii) Depth Distribution. The depth of each track in the emulsion at 
the scan line was measured, and a depth distribution plotted for sets 
(a) and (b) separately. There was evidence for the missing of tracks in 
the region of the emulsion nearest the glass in set (a), and so tracks found 
in the lower 10% of the emulsion were omitted and the scanning area 
appropriately reduced for the flux calculations. There was no indication 
of the missing of tracks at any depth in set (6). 

(iii) Length Distribution. The projected length of each a-particle 
track was measured and the distribution plotted for sets (a) and (0). 
In set (6), it appeared that some tracks whose length was only slightly 
greater than lcm had been lost, and so the minimum accepted length 
for the flux calculations was set at 1-lem. There was no evidence for 
such a loss in set (a). 

The tracks found by this scanning were aligned along the axis of a 
microscope and inspected in order to distinguish between singly and 
doubly charged particles. In those cases where this inspection did not 
result in an unambiguous identification, a measurement was made of 
the multiple scattering and grain density. 

Each track identified as having been produced by an «a-particle was 
traced back to the top of the stack to verify that it had not been produced 
by the break-up of a heavier nucleus. Those which were not incident 
on the stack were discarded. The remainder were then followed down 
into the stack until they either interacted or left. All «-particles on which 
there was more than 2cm of track length situated more than 1em from 
an outside edge were then examined further. A measurement of multiple 
scattering was made on each of these tracks employing techniques similar 
to those used previously, Fowler and Waddington (1956). <A basic cell 
size of 500, was used, and two and three times overlap calculated. 
Noise was eliminated by assuming that it varied as the square root of 
the cell size. The momentum, and thus the energy of each particle was 
then calculated using a scattering constant of 30-7 (Freier et al. 1959). 


§ 4. Tue Mean Free Pata 


The «-particles found in this experiment were followed through the 
stack until they interacted or left. While this tracing was primarily 
to ensure that each particle was incident on the stack, and that there 
was sufficient track length available for scattering measurements, it 
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was also possible to obtain a value for the interaction mean free path, X. 
In all, in the 2493cm of track length which was available for the 
observation of interactions, 103-5}; events were found, of which three 
Were Qin—Xout interactions. The resulting interaction mean free path 
was 24-0+2-4cm. This may be compared with the rather smaller value 
of A=18-7+0-6cm, which is the weighted mean of those values found 
previously by Waddington (1956a), Hanni (1956), Appa Rao et al. 
(1956), Shapiro et al. (1956) and Lohrmann and Teucher (1959). The 
reason for this apparent discrepancy is not clear. 

This mean value is calculated from the assumption that A is independent 
of energy between a few hundred Mev and at least 10 Bev per nucleon, 
since the median energies of the samples used by these investigators 
range between ~1-5 and 12Bev per nucleon. 

In order to investigate further this assumed independence with energy, 
the data obtained in this experiment have been separated, by the scatter- 
ing measurements, into those particles with H<1Bev per nucleon, 
and those with H>I1Bev per nucleon. For those with H<1 Bev per 
nucleon, median energy 0-6 Bev per nucleon, 38-5 stars were observed 
in a length of 659 cm, where this figure includes a correction for the length 
necessary for the measurement of the energy of each particle. 

The resulting value of A was 17-142:8cm. For those particles with 
H>1Bev per nucleon, median energy 2:65Bev per nucleon, 50 stars 
were observed in a corrected length of 1310cm, giving A= 26-2 +3-7cm. 
These values may be compared with those of 22:0+4-2cm and 
19-95+3-8cm found for samples also separated at 1Bev per nucleon, 
- with median energies of 0-70 and about 2-0 Bev per nucleon, respectively, 
by Waddington (1955). Combining these values gives A= 18-6 + 2:3 cm for 
E<1sev per nucleon and A=23-2+2-6cem for H>1Bev per nucleon, 
which are not significantly different. 


§ 5. FLUX 


The flux of «-particles crossing the scan line, J(x,%), for 0<60° was 
found, after correction for the ascent and descent of the balloon, to be 
105-0+6-6 «-particles/m?stersec. For %<30° it was 109:0+9-0 and 
for 30°<9<60°, 106-04 10-6 «-particles/m*stersec. These values were 
corrected for absorption in the overlying emulsion, using a mean free 
path of 20cm, and for absorption and production in the overlying 
atmosphere, using an absorption mean free path of 50g/cm*. The fluxes 
at the top of the atmosphere, J(0,%), were then: 


167 +10 w-particles/m? ster sec for 9 < 60° ; 
160 + 14 w-particles/m? ster sec for 3 < 30°; 


and 
183 + 18 «-particles/m? ster sec for 30° <9 < 60°. 


+ One track could not be traced from one emulsion to the next, and it was 
assumed that this was probably due to a oe undetected interaction. This 
was given a weight of one-half. 


P.M. 21 
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It is apparent that there is no significant difference between the fluxes 
for ><30° and for 30°<9<60°. Hence the value at the top of the 
atmosphere has been taken as: 

J = 167 +12 «-particles/m? ster sec 
where a possible error of 10% in the absorption mean free path in the 
atmosphere has also been taken into account. No allowance has been 
made for any error in the assumed thickness of the emulsions since 
measurements made on at least four stacks in the last two years have 
shown that the manufacturers have recently kept the average thickness 
very close to the nominal value. 


Fig. 2 
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The a-particle flux values as a function of Ottawa neutron monitor counts. 
Curves (A) and (B) show relations (1) and (2). 


It will be seen later, § 6.3, that this flux value is for «-particles having 
an energy above between 200 and 300Mev per nucleon. In table 1 this 
value is compared with other «-particle flux values measured at similar 
cut-off energies since 1955, Also given in this table are the hourly mean 
neutron counts recorded by the Ottawa neutron monitor on the day of 
each experiment. Figure 2 shows these flux values as a function of the 
neutron counts. A linear relation has been fitted to these data by 
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assuming that the errors in the hourly neutron counts, NV, are negligible 
compared with those in Jy. This relation has the form 
(J go — 189) + 7-2 = (2°73 + 0-41)(N — 285). mre ti) 
However, due to the preferential removal of lower energy particles during 
flux decreases and the greater effectiveness of high energy particles in 
producing counts in neutron monitors, the sensitivity of V to changes in 
J,. should increase as J,,. decreases. Hence it might be expected that 
the relationship between these two quantities is not linear but at least 
quadratic. The following quadratic was derived from the data: 
Jo — 189 = — (13-4 + 9-9) + (NW — 285)(2-48 + 0-41) 
+ (N —285)2(0-043+0-0024). . . . . . (2) 
These two relations are shown in fig. 2. 

Relation (2) is obviously not vaid for N $255, since then dJ,./dN 
becomes negative. However, for larger values of NV relation (2) appears 
to be at least as consistent with the data as relation (1), and has an 
added physical significance. 

The neutron counts are predominantly caused by the primary protons 
of the cosmic radiation, and thus the existence of a simple relation 
between the «a-particle fluxes and neutron counts can be interpreted 
as meaning that the ratio of the «-particle flux to the proton flux has 
remained constant over the wide range of observed values. An observa- 
tion of an a-particle flux associated with a neutron count which signi- 
ficantly diverged from these relations would suggest that, at the time 
of that particular experiment, this ratio had altered. Due to the rather 
wide experimental spread of the data the possibility that small variations 
may have occurred during these experiments cannot be excluded. 

The fact that the flux obtained here is somewhat higher, although 
barely significantly, than those obtained in 1957 and 1958 is supported 
by the observations made on the ending particles in these emulsions, 
Waddington (1959, 1960b). In this stack, exposed according to Quenby 
and Webber (1959) at a cut-off rigidity of ~1-65Bv, the rate of 
production of ending particles, p, was 486+35enders/em%h. This 
value may be compared with the very similar value of 


p=485 + 21 enders/em* h 
obtained in 1957-58 from stacks exposed at the lower cut-off rigidity 
of ~1-0Bv. 

The flux of primary protons has also been measured in this stack by 
examining the interactions produced by singly charged particles, 
Waddington (1960c)f. In the emulsions the flux of singly charged 
particles was determined as 1270 + 180 particles/m?stersec. When cor- 
rected to the top of the atmosphere using the curves given previously 
the flux was 810 +117 protons/m?ster sec, correction (d), or 

860 + 125 protons/m? ster sec, 


a a 
} Footnote added in proof.—Here a factor of 1-74 + 0-04 has been used to correct 
for the loss of small stars, instead of the value of 1-56 + 0-13 used previously. 
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correction (d). Accepting this latter value the proton—alpha ratio is 
5:15 +0-77, which is comparable with the value of 6-23+0-63 found by 
McDonald (1959) on a flight flown on February 16, 1958. 


§ 6. ENERGY SPECTRUM 


The energy values found by scattering, and corrected for the energy 
loss in the material between the top of the atmosphere and the point of 
measurement, are shown in fig. 3, separated into those with 9< 30° and 
those with 30°<%3<60°. These energy values have been used to 
construct integral and differential energy spectra. 


Fig. 3 


9) 
Se 
O 
< 
cc 
er 
Le 
O 
cc 
uJ 
fan} 
= 
=) 
ZnO) Gy << Soy 
5 ~ 1G 
az 
O Oe) ARO ile3) 2:0 


BeV PER NUCLEON 


The distribution in energy at the top of the atmosphere, divided into those 
particles with 3 <30° (a) and those with 30°<9<60° (6). Values which 
are lower limits are shown shaded. 


6.1. Integral Energy Spectrum 


Figure 4 shows the integral energy spectrum for all 9. The flux of 
a-particles, having an energy greater than some value 7’, J(2>T), is 
plotted as a function of the total energy, 7’=m,c?+ £, expressed in 
Bev per nucleon. Also shown are the integral energy spectra obtained 
by similar experimental techniques in 1957 (Freier et al. 1959) and 
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in 1950-54 (Fowler and Waddington 1956). The spectrum obtained 
in this experiment can be seen to lie between these two, which are 
respectively regarded as being characteristic of solar maximum and 
solar minimum. It should be noted that the errors shown on each point 
are not independent, since each point includes all those particles of 
higher energy. The apparent difference between this spectrum and that 
obtained in 1957 is only as significant as that between the two fluxes. 


Fig. 4 
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The integral energy spectrum at the top of the atmosphere obtained in this 
experiment. Also shown are the spectra obtained in 1957 and 1950-54. 


Assuming that these spectra may, over a restricted range of energies 
be represented by a simple power spectrum of the form J(>7)=C. T-», 
where C and y are constants, the values of these constants were deter- 
mined by putting this expression into differential form dJ=K .7-#dT 
where B=y+1 and K=—Q(.y, and considering the values of J (>7) 
at 7’'=1-43 and 2-13Bev per nucleon. The result, for all 9, was 
B=2-124+0-12, y=1:12+0-12 and O=222+31. These values may be 
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compared with those found in 1957 of y=1-:17+0-14 and C=185+30; 
and in 1950-54 of y=1-48+0-12 and C=360+ 40. 

Separating the particles into those with 9<30° and with 30°<9< 60°, 
B, over the same range of energies, is 1:94+0-14 and 2:4240-24 
respectively. The ratio of these is 0:80+0-10. An examination of the 
detailed energy spectrum shows that this apparent difference is due to 
additional low energy particles at large zenith angles, or a lack of 
such particles at small angles. For $<30° the fluxes above and 
below 1 Bev per nucleon are 111+11 and 49+8 a-particles/m? ster sec 
respectively, whereas for 30°<9<60° they are 103+13 and 80+12 
a-particles/m? ster sec. 

Thus, while the total fluxes at the top of the atmosphere of the two 
zenith angle groups agree to within one standard deviation, the exponents 
of the energy spectra differ by two standard deviations. The probability 
of observing such a difference by chance is about 1 in 20, which is not 
negligible. If this difference is real it could only be due to the change 
in cut-off energy with azimuth angle, which should result in a higher 
cut-off energy for particles coming from the east and a lower one for those 
coming from the west. Assuming that the stack rotated during the flight 
the net effect would be a function of the difference between these two 
directions. A consideration of the theory of these geomagnetic effects 
(Schwartz 1959), suggests that the net effect would be a decrease in 
the number of low energy particles at large zenith angles, which is the 
opposite to what is actually observed. Therefore it was concluded that 
not only was there no real difference between particles at different zenith 
angles but there was no physical justification for expecting a difference 
observable under the conditions of this experiment. 


6.2. Differential Energy Spectrum 

The differential energy spectrum is shown in fig. 5, where dJ, expressed 
in w-particles/m? ster sec 100 Mev, is plotted as a function of H, the kinetic 
energy in Bev per nucleon. The individual values are listed in table 2 
These values have been compared with values calculated from previously 
published data obtained in 1957 and 1958 by Freier et al. (1958)}, Aizu 
et al. (1959) and the joint Rochester-Oxford-St. Louis group (Engler, 
private communication). 

While the limited statistical weight of these results makes it difficult 
to draw any firm conclusions from this comparison, it appears that the 
spectrum obtained in this experiment is not too dissimilar from that 
observed previously. There is some indication of an increase in the 
number of particles above 1 Bev per nucleon, which could be evidence in 
favour of a magnetic modulation mechanism, but it is hardly significant. 
The divergence at low energies must be due to the operation of a cut-off 
mechanism. , 


EE ee a ee 1 een 
+ But not including that data obtained during the Forbush. decrease of 
August 31, 1957. 


528 G. RB. Stevenson and CO. J. Waddington on the 


Fig. 5 
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The differential energy spectrum at the top of the atmosphere in this experiment, 
shown by lozenge-shaped points, compared with the mean values found 
previously. The differential energy spectrum above 500 Mev per 
nucleon derived from the integral spectrum is shown as a dotted line. 


Table 2. Differential energy flux values 


dT x0 
dE a-particles/m? ster sec 100 Mev 
100 Mev 

Present experiment 1957-58 
per nucleon All ey eee ie ee eS aes! Previous 

9<30° 30°<9° <60° exp 
1-0-1-5 ee ae oa 2741-1 
1+5-2-0 a i -~ 6141-5 
2-0-2:5 2-0 + 2-0 2-0+2-0 ie 7T-S+1-2 
2-5-3-0 6:8 + 3-0 8-6+4-4 3-64+3-6 7841-2 
3-0-4-0 6-7 42-1 54+ 2-4 8-84+3-9 6-8 +09 
4-0-5-0 T4422 6-5 42-7 8-8+3-9 10-4 + 1-2 
5-0-6-0 10-8 + 2-7 9-7 +3-2 12-344-7 8-7+1-2 
6-0-7-0 9-5 42-5 54424 15-8 45:3 6-8+1-0 
7-0-8-0 11-5+2:8 8-6+3-1 15-8 45-3 67411 
8-0-10-0 5-141:3 4341-5 6-1+42°3 4-4+0-8 

10-0-12-0 5741-4 5-9+1:8 5-B+2-2 = 

12-0-15-0 6-441-2 6-S+15 53418 as 
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6.3. The Cut-off Energy 


In previous experiments to measure geomagnetic cut-off energies in 
emulsions (Waddington (1956b), Fowler and Waddington (1956); Freier 
et al. (1959)), the cut-off energy has been sufficiently great to justify 
the assumption that the integral energy spectrum rises steadily to the 
cut-off value. In this experiment y clearly varies appreciably before 
the cut-off energy is reached and information on the cut-off can only be 
obtained from a comparison of the differential spectrum with that 
obtained previously. Such a procedure is obviously much more sensitive 
to statistical fluctuations. 

Fig. 6 


Area of one track 


NUMBER OF TRACKS 


200 400 400 500 600 
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Histogram of energies at the top of the atmosphere lower than 600 Mev per 
nucleon. See text for explanation of curves. 


Figure 6 shows a histogram of those energy values below 600 Mev per 
nucleon, normalized to a curve obtained by plotting each value as a 
triangle of constant area and base width of two standard deviations. 
Also shown is a curve which, while not precisely normalized, represents 
a differential spectrum in this energy range of the form H*'°. The 
effects of cut-off energies at 300 and 250mev per nucleon are shown by 
curves c and 6 respectively, while that of the cut-off imposed by the 
overlying atmosphere is shown by curve a. It can be seen that while a 
sharp cut-off at 250mev per nucleon is in best agreement with the 
experimental distribution, it is impossible to unambiguously distinguish 
between this cut-off and one at 300 Mev or one lower than the air cut-off 
value. This result is not inconsistent with the value of ~310Mev per 
‘nucleon predicted by Quenby and Webber (1959). 
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ABSTRACT 


Precipitation from an alloy of gold +5 atomic % cobalt has been observed 
by electron microscopy using replication and transmission techniques. The 
precipitate particles form a rectangular grid of lines when observed on {100} 
faces and mainly parallel lines when observed on {110} faces. The simplest 
explanation is that dislocations on {110} planes are decorated by cobalt 
precipitation, These dislocation walls are presumed to have arisen from 
polygonization during ageing. Geometrical arguments are presented. 


In a study of the magnetic nature of cobalt precipitation in a gold matrix, 
we studied the course of the precipitation process by means of electron | 
microscopy using replication and transmission techniques. Single crystals 
of the gold + 5% cobalt alloy were solution treated at 950°c, quenched, and 
then aged at temperatures up to 450°c. Before replication, the specimens 
were etched, with gold being preferentially removed leaving surface 
bumps of cobalt rich precipitate. For {100} faced samples, the cobalt 
particles appear to precipitate along lines which intersect at 90°. These 
lines, within 20°, were along (110) directions. Figure 1|| shows this 
effect in a micrograph obtained by replication using silicon monoxide. 
In {110} faced specimens, the precipitate particles appear to come out 
along parallel lines; occasionally there are deviations with lines shifting 
by 20 to 35 degrees. Figure 2 is a good example of a replica electron 
micrograph obtained from {110} faced samples. 

Using transmission electron microscopy and diffraction on thin foil 
samples (rolled to a few thousands of an inch, homogenized, quenched, 
aged, electropolished) of the same composition, we determined the orienta- 
tion of the-rows of precipitate particles with respect to the gold matrix 
(fig. 3). The foil plane was usually {110}, sometimes {100}. This is in 


+ This research was supported in part by the United States Air Force through 
the Air Force Office of Scientific Research (ARDC). Communicated by the 
Authors. 

t Now at Crucible Steel Company of America, Pittsburgh, Pennsylvania. 

§ On Sabbatical Leave from Temple University, Philadelphia, Pennsylvania. 

|| All figures are shown as plates. 
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agreement with the rolling texture from face-centred cubic metals. For 
the {110} foils, the precipitate particles were always along a (111) direc- 
tion. As well as could be determined for the {100} foils, the rows of 
precipitate particles were in (110) directions. In some {110} cases the 
particles formed an almost perpendicular grid of lines, one of which was 
(111). It should be noted that the situation for foils differs in that they 
were severely cold-worked before the solution anneal. 

The regularities shown in figs. 1, 2 and 3 must be related to some sub- 
structure of the alloy and a plausible mechanism would be that precipita- 
tion occurs along dislocations. Such dislocation decoration has been 
noted by Castaing (1949) and by Wilsdorf and Kuhlmann-Wilsdorf (1954) 
working with aluminium—copper alloys, by Dash (1956) working with 
silicon + copper, by Mitchell (1957) with working silver halides, and by 
Amelinckx (1956) working with NaCl. 

Inasmuch as there is a real question as to whether the precipitates are 
decorating dislocations or whether the precipitates are producing dis- 
locations (Elbaum 1959, p. 237) we shall show the reasonableness of the 
former hypothesis. The production of dislocations may be due to vacancy 
precipitation, thermal or external mechanical shear strains, or to shear 
strains resulting from abrupt changes in composition. The density of 
dislocations in pure, slowly cooled gold foil was found to be 5 x 108/cm? 
(Sileox and Hirsch 1959). It is to be expected that the quenched speci- 
mens would show higher dislocation concentrations, especially since alloys 
should show higher vacancy concentrations (Sosin and Rachal 1959). 
Even if the only dislocations produced are those from vacancy migration 
to form discs and their subsequent collapse, it is expected that this would 
occur before cobalt precipitation. In gold—nickel alloys, gold diffuses 
faster than nickel (Reynolds et al., 1957) and, in a 5 atomic % nickel alloy, 
this means that vacancy diffusion is faster than that of nickel. Since the 
Goldschmidt radius of cobalt is about 5°% larger than that of nickel, we 
expect a similar situation in Au+5% Co. Also, it is rather difficult to 
conceive of a mechanism whereby precipitation would occur with some 
geometrical regularity without some guidance. 

A lower limit to the number of nucleation centres present, as measured 
by counting the number of precipitate particles after brief anneals, was 
about 5 x 10°/em*. Sileox and Hirsch (1959) obtained a value of 101°/em? 
for the dislocation density of quenched gold foil. However, since their 
pure gold was annealed at temperatures up to only 250°c and our gold 
alloys were annealed at 450°c before investigation, this difference is not 
unreasonable. Furthermore, the methods of observation were different, 
the compositions were different, and our number was obtained for solid 
single crystals and theirs for polycrystalline foils. 

There are several ways in which plane arrays of dislocations can be 
produced. One would be to have dislocations (produced by deformation) 
migrate to form simple sub-boundaries; these would tend to be planes 
perpendicular to the active slip planes (Burgers, 1940). For gold the 
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slip plane is {111} and planes perpendicular to {111} include {110}, {211}, 
{321}. Quenching strains have generated dislocations in silicon and 
germanium and it is expected that an even greater production would be 
obtained in metals (Elbaum 1959, p. 238). For a second mechanism we 
look to the collapse of vacancy discs to form dislocation loops on close 
packed planes. These have been observed in gold by Silcox and Hirsch 
(1959). Arrays of such loops would be expected to form on the {111} 
planes (Cottrell 1949) and again sub-boundaries could be formed 
(Amelinckx and DeKeyser, 1959). However, relatively high temperatures 
would be required for the needed dislocation mobility and it is unlikely 
that sufficient time is available during quenching. 

Our micrographs of replicas show precipitate lines which, according to 
the first mechanism given above, are probably traces of planes of dis- 
locations on the plane of observation. The square array of precipitate 
lines can be produced by having two planes of dislocations of type {100}, 
{110}, {111}, {211}. We shall assume that dislocations lie only on two 
planes of a given {hkl}. The requirement then, when we look at the 
micrographs for the {110} faced sample, is to obtain parallel traces from a 
pair of planes producing perpendicular traces on {100} faces. This cannot 
be done with {211} planes. Our calculations also show that arrays on 
{110} faces might include lines at 54-1° and 70-5° for {110} and {111} traces 
and 90° for {100} traces. Small line segments at angles near 55° and 70° 
were actually observed. Square arrays of precipitate lines on {100} faces 
have been seen in deformed silver halide (Mitchell 1957) and alkali 
halides (Amelinckx and Dekeyser 1959) with the traces in (110) direc- 
tions. Mitchell states that, in his silver halides, slip has occurred on 
(101) and (011) planes and that dislocations concentrate on these planes. 
In the case of fig. 1, a possible explanation is that slip has occurred 
along (111) planes, almost horizontal traces, and that there are poly- 
gonized walls along orthogonal (110) planes (Friedel, 1960). In this case, 
the traces will also be orthogonal. However, it is not clear why the 
‘precipitate decoration in the slip planes and in the polygonized walls 
should then be so similar. 

In the transmission electron micrographs rows of particles in (111) 
directions were observed for {110} foils; this requires that the plane 
producing the trace satisfy h+k+/=0. This condition is met by planes 
of the type (110), (211), (321) but not by {100} or {lll}. The (110) 
directions for {100} foils give us the condition h+k=0; this is met by 
(001), (110), (111), (112) planes. Thus it seems most probable that the 
dislocation lines are found in {110} or {211} planes. This would mean that 
the dislocation walls are probably produced as a result of distortion 
during quenching and subsequent polygonization. 

We believe the same geometry is effective in the case of quenched 
single crystals of the alloys. There is good agreement between figs. 2 
and 3; this indicates that the reliability of the replicas is good. However, 
since the directions obtained in these replicas were good only to 20°, we 
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cannot be sure that the same mechanism is effective. The angular 
separation between traces of (110) and (111) planes on {110} planes is only 
18°. It may well be that both mechanisms (polygonization and collapse 
of vacancy discs) are operative and are responsible for the complexities 
of some precipitate patterns. We have used{ 110} above rather than {211}; 
only when selected {211} planes are used do we obtain the results as given 
above. 

It may be of interest here to state some results obtained through mag- 
netic measurements. Torque and rotational hysteresis data indicate that 
the particles were crystallographically aligned, that the initial precipitate 
was face-centred cubic cobalt coherent with the gold matrix, and that 
with increasing time of ageing the f.c.c. cobalt transformed to the hexa- 
gonal variety and the particles become elongated along the c axis. 

Our conclusions are that precipitates are present where we would expect 
to find dislocation walls, that the precipitate particles probably nucleated 
on dislocations, and that the dislocation walls are probably the result of 
polygonization. Furthermore, we believe that most of the dislocations 
originated as a result of quenching strains and that the polygonization 
occurred early during the ageing process. 
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ABSTRACT 


Precipitation of carbon from «-iron during irradiation and thermal 
ageing has been studied using thin film electron transmission microscopy. 
During irradiation at 100°c precipitates formed with a density of 2 x 1014 em-3 
and saturated at 400A diameter after 72 hours in a flux of 


1-5 x 1011 neutron em-? sec-1> 1 mev. 


The precipitates were in the form of plates parallel to {100} and occurred 
individually in the matrix and in rows on dislocation lines. Thermal 
ageing at 100°c produced clusters of plates with a density of 3x 1013 em=-3 
and rows of plates on dislocation lines. At higher ageing temperatures 
the plate-like precipitates were replaced by larger dendritic particles at - 
200°c and needle-shape particles at 250°c with a density of 3 x 1012 cm-*. 
The defects produced during irradiation, which are responsible for irradiation 
hardening, were not detected. Tensile experiments showed that the harden- 
ing was most pronounced when precipitation did not occur. 


§ 1. INTRODUCTION 


A great deal of empirical information has been gathered concerning the 
hardening effects produced by neutron irradiation of various structural 
steels and this work has been reviewed recently by Harries (1960). 
However, relatively little is known about the actual mechanism of 
hardening, which results in an increase in yield stress and a rise in the 
temperature of the transition from brittle to ductile fracture. It has 
been shown from tensile hardening experiments (Churchman e¢ al..1957, 
Hull and Mogford 1958) that in a low carbon steel the hardening can be 
accounted for by an increase in the resistance to dislocation movement 
rather than by an increase in the stress to unpin a dislocation from its 
atmosphere, and it is suggested that this is due to the interaction of the 
moving dislocation with radiation induced defects, e.g. clusters of point 
defects. Irradiation hardening is also produced in pure iron (Kunz 
and Holden 1954, Konobeevsky e¢ al. 1955). In the work of Kunz and 
Holden it was found that the hardening annealed out with an activation 
energy of 3-1 ev suggesting a process involving self-diffusion. 


+ Communicated by the Authors. 
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The development of thin film transmission electron microscopy has 
enabled direct observation of collapsed vacancy aggregates as dislocation 
rings such as those found in face-centred cubic metals quenched from 
just below their melting point (Hirsch et al. 1958). Closed dislocation 
loops ~ 3004 diameter have been observed in irradiated copper (Silcox 
and Hirsch 1959, Makin et al. 1961, Barnes and Mazey 1960). However, itis 
difficult +o correlate the observed irradiation hardening with the 
formation of dislocation loops. 

The present paper describes work originally undertaken to examine 
the nature of the irradiation defects which resist the movement of 
dislocations in ferrous materials. After irradiating iron for only a few 
hours, observation of thin films in the electron microscope revealed a 
large number of defects 1004 to 300A diameter many of which gave a 
contrast similar to that produced by dislocation rings. However, further 
experiments indicated that these defects were a product of irradiation 
and thermal ageing at the pile temperature, and it is believed that they 
are produced by precipitation of interstitial atoms, e.g. carbon, from 
solid solution.. This has led to a comparison of precipitation of carbon 
in iron with and without the influence of neutron irradiation. The effect 
of precipitation and irradiation on the tensile properties of the iron has 
also been determined. 


§ 2, EXPERIMENTAL DETAILS 


The iron used in this work was supplied by the National Physical 
Laboratory. The analysis is given in table 1. 


Table 1. Analysis of material used (all values wt. °%) 


C Si sol. | Si insol. S ie Mn 
0-004. 0-007 | <0-001 | 0:0056 | <0-001 | <0-005 
Cr Ni Cu Mo N O 
<0:001 | <0-001 0-001 | 0-003 0-001 0-0005 


Al sol. | Al insol. 
0-003 0-001 


The 5 in. diameter hot rolled bar was cold swaged to 0-040 in. diameter 
rod with intermediate anneals for 1 hour at 900°c in vacuo followed by 
slow cooling, after every 35°% reduction in area. Foil for transmission 
electron microscopy was prepared from the rod by cold rolling to 
0-001 in. This was then thinned electrolytically using a solution of 
chromic acid in glacial acetic acid and the final foils examined in the 
Siemens Elmiskop Ib electron microscope at 100 kv. It was not possible 
to determine the final thickness of the foil but in calculating the number 
of precipitates and concentration of defects a foil thickness of 30004 
was assumed (McLean 1960). The electrolytic thinning was carried out 
after the thermal and irradiation treatments. 
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Tensile tests were made on the 0-040 in. diameter rod in a hard beam 
autographically recording tensile machine (Adams 1959) with a strain 
rate of 5:7 x 10-4 sec-!. 

The final heat treatment given to the foils and tensile specimens 
before ageing or irradiation was either (i) furnace cool from 900°C or 
(ii) air cool in the evacuated silica tube from 900°c. In the following 
sections these will be referred to as ‘slow cool’ and ‘fast cool’ 
respectively. 

The specimens were irradiated with neutrons in BEPO graphite 
moderated reactor in a flux of 1 to 2x 10" neutron em? sec-!> 1 Mev 
to doses between 2-7x 10% and 7-8x 10!* neutron cm-?>1 Mev. The 
temperature of irradiation was between 60°c and 100°c. Thermal 


ageing treatments to simulate the irradiation conditions were carried 
out in vacuo. 


§ 3. EXPERIMENTAL RESULTS 
3.1. Precipitation During Irradiation 


Foils examined in the electron microscope before irradiation contained 
a low density (<10*cm~*) of dislocations which appeared quite free 
from precipitation. In the fast-cooled material small black spots were 
observed in the matrix in addition to the dislocations. The density of 
the spots was 8-5x1042cm~? and they were about 100A diameter. 
Some of the defects had contrast similar to dislocation loops. 


Table 2. Effect of irradiation on diameter and number of defects 


Time of irradiation Dose, Diameter of Number of 
hours neutron cm~* precipitates precipitates 
> 1 Mev A ema 

5 27x 1018 100 1-7 %:10"4 

12 6-5 x 1915 130 ah Le 

24 i320 1" 170 2-0 x 1044 

72 ao x 10% 380 18310 

144 778 1018 380 273. x10 


When foils were annealed at 900°c for $ hour followed by a fast cool 
and then irradiated, precipitation could be induced after very low doses. 
The first definite signs of precipitation, in the form of small black spots 
~1004 diameter, were formed after 5 hours, with a density of nearly 
2x10%cm-%. The damage produced after 12 hours is shown in fig. If. 
After 24 hours the spots could be resolved into irregular loops. The 
diameter of the loops increased with the irradiation time up to 72 hours 
but the density of loops was approximately independent of time (table 2). 


i 


+ Figures 1-9 are shown as plates. 
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In fig. 2 dark streaks ~350A long can be seen in addition to loops 
and it is assumed that these are produced when the plane of the loop 
is normal to the plane of the foil, i.e. the loops appear edge-on. The 
appearance of the loops depends on the orientation of the foil and the 
results can be summarized as follows : 


(a) In {100} foils two sets of streaks at right angles parallel to {100), 
no loops (fig. 3). 

(b) In {110} foils one set of streaks parallel to (110) and loops : 
the ratio of streaks to loops is approximately unity (fig. 2). 


(c) In {111} foils only loops are observed. 


In addition to apparently homogeneous precipitation throughout the 
material precipitation occurred on dislocations. After 5 hours’ irradiation 
some dislocations were surrounded by an indistinct envelope ~ 100A 
across and after 12 hours definite precipitation could be seen (fig. 1). 
The appearance of the precipitate varied with the orientation of the foil 
and in some foils after 72 hours irradiation there was definite formation 
of a row of loops along the dislocation which appeared similar to the 
loops formed independent of the dislocation, fig. 4 at A-B. Figure 4 
shows a network of dislocations at C-D in which some segments have 
complicated precipitates whilst others are completely free suggesting 
that there is preferential precipitation on dislocations with certain 
Burgers vectors. There is a zone around the network, about 5 x 10-° em 
wide, free from loops and the precipitation on dislocation A—B is clearly 
more pronounced away from the network. 

When foils were annealed at 900°c for } hour followed by a slow cool 
and then irradiated no precipitation effects could be detected after 
six days. The dislocations were free from precipitation prior to irradiation. 


3.2. Precipitation During Thermal Ageing 


All the experiments were made on foils fast cooled from 900°c. The 
first ageing experiments were conducted at 60°c which was believed to 
be the most likely pile irradiation temperature. After three days the 
foil contained precipitates having an average diameter of about 4004 
with a density of 2x 10'%cm-%. A foil given such a treatment was cut 
in half and one half aged for a further three days at 60°c. The precipitates 
were less regular and had increased in diameter, however the density of 
precipitates was approximately the same (table 3). 

The other half of the foil was irradiated for three days and fig. 5 shows 
that there were two forms of precipitates ; large precipitates, somewhat 
larger than those in the foil aged at 60°c for six days, but with the same 
morphology, and fine precipitates (table 3). The growth of the large 
precipitates to 6254 diameter during irradiation compared with 5004 
for the six-day aged specimen suggests that the pile temperature is higher 
than 60°c and subsequent ageing treatments were carried out at 100°C. 
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In material aged at 100°c for successively longer times the formation 
of defects was somewhat similar to ageing at 60°c. After 1 hour small 
loops about 350A diameter were formed and after three days the 
precipitates had a complex morphology (fig. 6). Table 4 gives details 
of ageing at 100°c. Dislocations were decorated with precipitates and 


Table 3. Ageing and irradiation at 60°c 


Diameter of Number of 
Treatment precipitates precipitates 

A em-* 
Aged 72 hrs at 60°c 400 ya A Wie 
Aged 144 hrs at 60°C 500 355108 
Aged 72 hrs at 60°C Sf 625 oe sel Gt 
followed by 72 hrs irradiation 150 3x 10" 


Table 4. Ageing at 100°c 


Time, hours | Diameter of precipitates, A | Number of precipitates, em~% 


1 350 2x 108 
72 800 3 x 108 
120 925 3 x 10% 


again their appearance depended on orientation. When the foil was 
close to {111} the precipitates in the matrix often appeared as a cluster 
of plates with six-fold symmetry (fig. 6). In {100} foils the precipitates, 
either as individual particles or on dislocations, appeared as streaks in 
two directions at right angles (fig. 7) along (100) in the foil. 


3.3. Annealing of Foils Containing Precipitates 


Foils which contained precipitates formed by ageing for three days at 
100°c or by three days’ irradiation were further heated for 1 hour periods 
at 140°c, 175°c, 200°c and 250°C respectively before thinning. In the 
irradiated foils no change occurred after 1 hour at 140°c, but after 1 hour 
at 175°c the precipitates had become larger and more complicated and 
the concentration reduced significantly. At 200°c dendritic particles 
were formed and at 250°c the precipitates assumed the form of 
Widmanstatten-like needles. These changes are illustrated in figs. 8 (a—c) 
and the changes in shape and concentration are listed in table 5. In the 
aged foils no change was observed below 200°C at which temperature 
the precipitates appeared more dense and became elongated, and at 250°C 
the structure was identical to the irradiated foil. In foils parallel to {110} 
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the precipitate needles formed at 250°C were at 60° to each other and 
were along (111). 


Table 5. Effect of heating foils containing precipitates 


Diameter of} Number of 


Treatment precipitates | precipitates Shape 
A eni-= 
Three-day irradiation 350 oom Les Dises 
+1 hour 140°c 350 aroLtre Discs 
+1 hour 175°c 550 9x10" Complicated aggregates 
+1 hour 200°c 5 x 10” Dense dendrites about 
15004 long and 6504 
across 
+1 hour 250°o 1 x 10?” Needles or rods about 
| 50004 long and 3504 
across 
Three-day thermal 
ageing at 100°C 800 a1 Complicated aggregates 
+1 hour 175°c 850 2% 10 Complicated aggregates 
+1 hour 200°c 2x10" Dense dendrites about 
16004 long and 6004 
across 
+1 hour 250°c 3x 10" Needles or rods about 
5000 A long and. 2004 


across 


A number of foils were irradiated and thinned and then heated on the 
hot stage in the electron microscope. Figure 9 (a) is a typical field prior 
to heating which contains loops with and without a black spot in the 
centre, clusters of loops and decorated dislocations. Heating for a few 
minutes between 150°c and 180°c produced no change in the structure 
but at 180°C the single loops were observed to shrink very rapidly (e.g. 
loops at A, fig. 9) and disappear within 1 or 2 minutes. The decoration 
on the dislocation lines also disappeared in a similar way. The more 
complicated loops at B in fig. 9 shrank slowly, leaving behind a black 
spot (often associated with a small piece of dislocation line) which was 
still visible after heating to 235°c. No massive precipitation was observed 
at this temperature similar to that formed in the specimen heated outside 
the microscope before thinning. 


3.4, Tensile Experiments 


A limited number of tensile tests were made to determine whether any 
relation exists between the irradiation hardening and the formation of 
irradiation induced precipitation. The increase of yield stress Aocy 
for tests at —78°o are recorded as a function of irradiation or thermal 
ageing time in fig. 10. For specimens annealed at 900°c and slow cooled 
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Acy increased continuously with the irradiation time but for specimens 
annealed at 900°c and fast cooled, there was little change in Acy between 
5 and 24 hours’ irradiation. In specimens aged at 100°c after a fast cool 
from 900°C no appreciable increase in yield stress was detected after 
72 hours. 


Fig. 10 


ie) 


(@)) 


Aa, dynes cm™* x10® 
iN 


O 1OmmeOursOs5407 2502 «GO, 70.80 
Irradiation or ageing time (sec) 


Increase in yield stress, Ac,, as a function of irradiation or ageing time. 
A Specimens fast cooled from 900°c and irradiated, curve (i) ; 
® specimens slow cooled from 900°C and irradiated, curve (il) ; 
x specimens fast cooled from 900°C and aged at 100°C, curve (iii). 


§ 4. Discussion 
4.1. Precipitation During Irradiation 


Irradiation of iron for 72 hours after a fast cool from 900°c, produced 
a high density of loops, 2x 10 cm~?, about 400A diameter. The loops 
gave a contrast effect similar to loops of dislocation line produced in 
aluminium and copper after quenching from just below the melting 
point and in copper after neutron irradiation, i.e. they consisted of a dark 
ring with little or no contrast in the centre and the contrast, could be 
altered by rotating the specimen in the stereo holder of the microscope. 
In copper and aluminium the dislocation loops are formed by the collapse 
of platelets of vacancies and the present observations were explained 
initially in a similar way ; thus, if it is assumed that each primary 
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knock-on in iron produces 100 Frenkel pairs, an irradiation of 4x 1016 
fast neutrons cm-2 would produce a concentration of approximately 
10-5 point defects. The concentration of defects that would be required 
to produce the observed number of loops after such an irradiation is 
approximately 5x 10-5. However, there are a number of observations 
that indicate that the loops in iron are not collapsed vacancy aggregates. 
Firstly, the number of loops saturates after only 5 hours’ irradiation and 
the diameter of loops saturates after 72 hours, whereas, no saturation 
would be expected if the loops were formed by point defects created 
continuously during irradiation. Secondly, loops were not observed 
in slow cooled material after irradiation. Thirdly, somewhat similar 
loops were observed in foils aged at the pile temperature in the absence 
of irradiation. 

We propose, therefore, that the loops are produced by precipitation 
of interstitial carbon atoms in a thin layer probably one atom thick. 
The appearance of the loops in foils of different orientation is consistent 
with them lying parallel to {100} planes} in the iron lattice; hence 
in {110} foils one set of loops is seen edge-on parallel to (110) and in 
{100} foils two sets of edge-on loops occur at right angles, parallel to 
{100}, in {111} foils no loops appear edge-on. In fig. 2 three streaks 
appear parallel to <100), e.g. at A, and they are probably produced 
by loops lying on {100} planes inclined at 45° to the plane of the foil 
which have intersected the foil surface, so that only a small portion of the 
precipitate lies within the foil The most favourable interstitial site in 
the body-centred cubic cell is 0, 0, } and if all such sites in a (100) plane 
are occupied by carbon atoms the iron atoms parallel to this plane will 
be displaced about 1A normal to the plane. The arrangement of atoms 
around such a layer will be somewhat similar to a closed loop of edge 
dislocation in the (100) plane. 

Since the solubility of carbon in iron at 100°c is very small, ~ 10-7 wt. %, 
practically all the carbon will precipitate from solid solution. After 
72 hours’ irradiation, precipitation appears to be complete and assuming 
that the precipitates are monatomic layers of carbon, the density and 
size of the loops are equivalent to approximately 0-001 wt. % carbon. 
This is close to the solubility of carbon expected in this material fast 
cooled from 900°c. Thus it is considered that precipitation under 
irradiation conditions consists of the growth of a monatomic plate of 
carbon atoms. 

In the presence of grain boundaries and individual dislocations, the 
distribution of the loops was considerably modified. The free zone 
around dislocations, low angle boundaries and grain boundaries presumably 
results from the preferential diffusion of carbon to these sites. In many 
cases the precipitation on the dislocations can be resolved into definite 
i 


+ This is a correction to a preliminary report of this work (Hull and Mogford 
1960) which stated that the loops were parallel to {1 10} in the iron lattice. 
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particles which have the same contrast appearance as the individual 
precipitates, e.g. dislocation A-B in fig. 4, and it is possible that these 
precipitates also consist of monatomic layers of carbon. The decoration 
on dislocations is similar to the effects produced during thermal ageing 
which will be discussed below. 


4.2. Precipitation During Ageing at 100° 

When specimens were aged at the same temperature at which they 
were irradiated, precipitation occurred on dislocations and throughout 
the lattice. At the dislocations the first sign of precipitation was a dark 
zone previously observed by Hale and McLean (1960, private communica- 
tion) and after 24 hours this could be resolved into individual particles. 
In foils with {001} orientation the particles appeared as short streaks 
on two sets of planes at right angles to each other along (100) directions 
indicating that they consisted of small plates parallel to {100} planes. 
The precipitates formed away from dislocations were more complicated 
than in the irradiated specimens and had a much lower density, 
3 x 1018 particles cm~3. They appeared to consist of a cluster of plates 
and again in {100} foils aggregates with plates at right angles could be 
seen. It is possible that the clusters of plates formed by successive 
growth of individual plates and this explains the formation of clusters 
containing 2, 3, 4, 5 and 6 plates in fig. 6. 

The similarity in the appearance of the particles formed on irradiation, 
the particles on dislocation lines in both aged and irradiated foils and the 
individual particles in the clusters formed in aged foils suggests that the 
plate form is common to all. In addition the size of the unit plate is 
practically constant, i.e. about 400-5004 diameter in all cases. If it is 
assumed that the platelets are monatomic layers of carbon atoms, 
comparison of the aged and irradiated foils, when all the carbon is 
precipitated, i.e. 72 hours, requires that the average number of plates 
in the clusters in the aged foils should be the ratio of the number of 
precipitates in the irradiated and aged foils respectively, 

2x 101M 
Le. 5018 = 
which agrees with the observations. 


4.3. Nature of the Nucler 


One of the most important differences between ageing and irradiation 
is the difference in the number of particles. This is most clearly illustrated 
by fig. 5 which shows that a concentration of 3 x 101 particles cm”? is 
formed on ageing at 60°c and that irradiation produces quite independently, 
a concentration of 2 x 10" particles cm~*. The concentration of precipi- 
tates produced during irradiation was constant after about 5 hours. The 
most probable source of nuclei is displacement spikes. The arrangement of 
point defects in this spike can be represented as a small region with 
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a large number of vacancies ~ 10-100 surrounded by a much wider 
region containing interstitials. This central region may provide a ready 
sink for carbon atoms. The number of displacement spikes, calculated 
as the product of the total neutron flux, density of atoms in the metal, 
and the neutron cross-sectional area for elastic scattering, produced in 
5 hours is about 2x 10!4 which agrees closely with the observed number 
of particles. Presumably further nuclei produced are unstable with 
respect to the previously formed nuclei and carbon precipitation continues 
on the existing precipitates. There is good evidence that small precipitates 
can redissolve during irradiation from the width of the free zone around 
dislocations and sub-boundaries. After 12 hours the zone is about — 
10-5 cm wide and after 72 hours it is about 5x 10-° cm. 

The width of the free zone is close to that expected from the diffusion. 
rate of carbon in iron. The diffusion distance x = (Dt)"? can be calculated 
from ¢ the diffusion time in seconds, and D which is given by : 

D=D, exp (—AH/RT). 

Using the values AH =20 400 cal/g mol and D,)=3-16 x 10°? cm? sec + 
obtained by Thomas and Leak (1954) at 100°c, D=4x 10-™ cm? see. 
Thus after 12 hours x=4100A and after 72 hours x=10 000A. However, 
diffusion may not be the rate controlling factor in the deposition of 
carbon atoms on platelets. When the concentration of precipitates is 
2x 1014cm~3, i.e. after irradiation, the average distance apart is about 
17004 so that the growing plate will require a diffusion distance of 
850A for complete precipitation. Precipitation is only complete after 
72 hours, v=10 OOO A. 

The nature of the nuclei in the aged foils is less certain. There appears 
to be a low density of defects after cooling and after 1 hour at 100°c 
there are 2x 104% cm? loops about 300A diameter. It is possible that 
the nuclei are quenched in from the higher temperatures either as 
aggregates of carbon atoms or clusters of vacancies. As in the irradiated 
specimens the distance between particles, i.e. 32004, is much smaller 
than the diffusion distance estimated from the diffusion coefficient and 
activation energy, although the exact time when precipitation was 
complete in aged specimens was not determined. 


4.4. Hffect of Heating Foils and Ageing at Higher Temperatures 


Although the experiments in which specimens were aged at higher 
temperatures are far from complete the results in table 5 clearly show 
that a marked change in the character of the precipitation occurs between 
175°o and 250°c. In the irradiated specimens the loop defects, density 
3x 10! cm™, were replaced at 175°c by complicated precipitates, and 
with dense dendritic precipitates after ageing for 1 hour at 200°c. After 
an hour at 250°c the precipitates had grown into rods 50004 long and 
350 A diameter lying along (111) directions. It is not possible to establish 
whether the rod-shaped precipitates, which are probably Fe,© (Leslie 
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et al. 1959), form as a completely new precipitate or by a rearrangement 
of atoms of the existing precipitates, although the concentration is 
reduced by about 10?. No appreciable change is observed in the aged 
foils held at 175°c for 1 hour but at 200°c and 250°c the morphology 
and density of the precipitates are identical to those in the irradiated 
foils. We can therefore recognize two or three stages in precipitation, 
(1) ageing at 100°c results in the formation of platelets, individually or 
in clusters, depending on the number of nuclei available, throughout 
the lattice and in rows along dislocation lines, (2) between 175°c—200°c 
these dissolve and complex dendritic precipitates are formed, and (3) at 
250°c a Widmanstatten array of rods is formed. 


4.5. Irradiation Hardening 


The present experiments have provided little information on the nature 
of the defects producing an increase in cy onirradiation. It is immediately 
clear from fig. 10 that precipitation is not producing an appreciable 
increase in hardness during irradiation ; for specimens in curve (i) 
precipitation is expected, but there is no increase in cy between 5 and 
24 hours’ irradiation, whereas in curve (ii) for specimens in which no 
precipitation could be detected after 14 days, cy increases continuously 
with irradiation time. This implies for curve (i) that either the irradiation 
hardening is cancelled out by softening induced by ageing, or that the 
irradiation induced point defects, which produced hardening in curve 
(ii), are not effective in the presence of precipitate particles. Curve (iii) 
suggests that no softening occurs when specimens are aged, so it is 
possible that the point defects become associated with the carbon platelets 
and do not produce hardening. 

Since nothing can be seen in slow cooled foils after an irradiation of 
14 days, even though there is an appreciable increase in oy, it follows 
that the irradiation induced defects are extremely small. Further work 
is necessary to determine whether it is possible with higher doses to 
observe irradiation induced defects in material free from. carbon or 
nitrogen. It is also important to make a detailed correlation, in the 
early stages of precipitation, between the tensile properties and the 
appearance of precipitation on dislocation. 


§ 5. CONCLUSIONS 


1. The defects produced during neutron irradiation of iron at 100°c, 
which increase the resistance to the movement of dislocations and hence 
increase to yield stress, have not been detected in the electron microscope 
for the irradiation doses used. 

2. Irradiation creates nuclei for the precipitation of carbon from solid 
solution and they are probably produced at the disorder caused by dis- 
placement spikes. It is possible that irradiation will provide nuclei for 
precipitation in other supersaturated Syston and hence affect the rate 


of transformation. 
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3. At 100°c carbon precipitates in the form of platelets which grow 
in size to about 40048. When a large number of nuclei are present, as 
in the irradiated material, the platelets form individually and when there 
are fewer nuclei the platelets form in clusters. Dislocations provide 
nucleation sites and the platelets form in rows along the dislocations. 

4. In the electron microscope the platelets, which lie along {001} 
planes in the iron lattice, produce a contrast similar to closed dislocation 
loops, and it is suggested that their structure consists of a monatomic 
layer of carbon in the 0, 0, } interstitial sites. The amount of carbon, 
estimated from the size and density of the plates when precipitation is — 
complete, agrees with the amount expected in solid solution before 
precipitation. 

5. Whilst it has been assumed in the discussion that only interstitial 
carbon atoms are precipitating, the experiments do not differentiate 
between carbon and nitrogen. 
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ABSTRACT 


A micro-beam x-ray technique has been applied to the fracture surfaces 
produced by slow-growing fatigue cracks. The cracks were propagated 
in metal sheets subjected to pulsating tension such that the rate of growth of 
the crack was proportional to its instantaneous length. The fracture 
surfaces showed a highly developed sub-grain structure with large mis- 
orientations >13°. The sub-grain size was found to be independent of the 
range of cyclic stress used to propagate the crack and to be characteristic 
of the metal. The progressive development of the sub-grain structure 
was followed in torsion specimens subjected to large ranges of plastic-strain 
~10-. 

If the process of sub-grain formation is regarded as an essential feature 
of the mechanism of fatigue crack propagation then the order of susceptibility 
of metals to fatigue crack growth and the phenomenon of non-propagating 
surface cracks can be interpreted. 


§ 1. INTRODUCTION 


THE emphasis in many recent physical investigations of the effects of 
cyclic stressing has been upon the origin of the fatigue crack. Ina series of 
experiments Thompson ef al. (1956) have shown, for example, how the 
life of a copper specimen is correlated with the development, at an early 
stage, of persistent slip bands which eventually behave as cracks. The 
metallographic changes accompanying the cyclic stressing of aluminium and 
its alloys have been classified by Forsyth (1952) and several structural 
conditions which are associated with the nucleation of surface cracks have 
been identified. In particular, the coarse slip striations, which are 
characteristic of cyclic deformation and which have been described by other 
workers (Hempel ef al. 1955, Wood 1956), appear to cause a surface profile 
which can be identified with the initial crack as calculations by May (1960) 
suggest. 

Although, therefore, the nucleation of fatigue cracks at an early stage of 
the life of the specimen has been related to the behaviour of cyclic slip bands 
at free surfaces the limitations of the resultant physical viewpoint are 
seriously felt when the growth of the initial cracks is considered. The 
transition from the formation of a surface crack to the commencement of 


+ Communicated by the Author. 
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steady crack growth is particularly important since analyses of crack growth 
have shown that the slow, reproducible growth which occurs at low cyclic 
stress ranges, ceases if either the stress range or crack length is too small. 

This phenomenon of non-propagation of fatigue cracks has been shown 
by Frost and Dugdale (1958) to underlie the behaviour of notched mild 
steel fatigue specimens. Notches with high theoretical stress concentra- 
tion factors, Kt, form non-propagating cracks at the notch root at alterna- 
ting stress ranges given by the intrinsic fatigue limit divided by Ay. Fracture 
of the specimen is, however, only produced by raising the alternating stress 
to a propagation value which is characteristic of the material and the length 
of the surface crack formed at the notch root. 

In the present paper, experimental observations resulting from the 
application of an x-ray micro-beam technique to the fracture surfaces 
formed by slow fatigue crack growth, are presented. The cracks were grown 
in centrally notched sheets under pulsating tension at ranges of applied 
nominal alternating stress o, and mean stress, such that the crack growth 
obeyed the relations found by Frost and Dugdale (1958) : 


dl 
see ome 
ep. 
a menemrtie rs eM yy 
in- =k, 


0 

where /=current half-length of crack, N=number of cycles, /)=initial 
half-length. The coefficient of proportionality k is given by o?/N, where 
N, depends upon the metal and is independent of mean stress for the metals 
and mean stresses employed in the present observations. 

The observations (on aluminium, mild steel and zinc) showed that the 
fracture surfaces produced under these conditions possessed a highly 
developed sub-grain structure characteristic of the metal and involving 
large misorientations. The structure was independent of the alternating 
stress range used to produce the fracture and of the initial state of working 
of the material. In order to follow the gradual development of the sub- 
grain structure prior to the state characteristic of fracture it was considered 
that torsion specimens subjected to conditions of large plastic strain range, 
€p] (giving lives up to 5 x 104 cycles) would be most suitable. Under these 
conditions Coffin and Tavernelli (1959) have shown that the relation 


ep! V?=const. holds, a relation which is again relatively free from the 
scatter characteristic of plain fatigue testing. 


§ 2. ExprRIMENTAL Metuop 


The x-ray micro-beam technique introduced by Hirsch and Kellar 
(1952) depends for its application upon the dislocations in the deformed 
metal being so distributed that they form arrays or boundaries, thus 
defining sub-grains which are relatively free from dislocations. Informa- 
tion about the distribution of the dislocations can be obtained if the volume 
irradiated is not too large compared with the sub-grain dimensions and if 
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the difference in dislocation density between sub-grain and boundary is 
sufficiently large. The application of the method to metals deformed by 
rolling by Gay et al. (1954) required the use of very fine micro-beams but, as 
will be shown below, the structures formed by cyclic stressing involve larger 
and better resolved sub-grains, factors which enable the technique to be 
applied in a relaxed form using beams down to 70 u diameter only. 

When the deformed metal is irradiated with a sufficiently fine beam of 
x-rays the Debye-Scherrer ring is observed to be made up of individual 
spots. If the spots are well-enough defined to be counted, the mean sub- 
grain size can be deduced since the number of reflections on a Debve— 
Scherrer ring is given by (Hirsch 1952): 


W=3= peos0 (d+) 


where A=area of cross section of beam, t=effective penetration of 
beam, v=volume of sub-grain, p=multiplicity factor, 0= Bragg angle, 
d@=divergence of beam, A=angle over which the sub-grain reflects 
due to distortion, small size and wavelength spread. 

Further information about the sub-grain in addition to its size can be 
obtained from the radial and tangential spread of individual spots on the 
Debye-Scherrer ring and from the contrast between a spot and the general 
level of background in the Debye-Scherrer ring. When the volume irradia- 
ted is such that only a few original grains are included, the Debye—Scherrer 
ring may be incomplete and consist of separate arcs. The total range of 
misorientation within the grain, due to deformation, can then be estimated 
since sin B/2 =cos 6 sin 4/2 (Hirsch 1952) where y# is the angle subtended by 
the arc at the centre of the ring and f is the maximum angle between normals 
of the reflecting planes of sub-grains within the original grain. 


§ 3. EXPERIMENTAL OBSERVATIONS 
3.1. Aluminium 


In their initial rolled condition the }in. thick aluminium sheets of 
commercial purity (99-8%), used for crack growth, had a wide variation of 
grain size ranging from 20-40. at the outside surfaces to ~ 300 in the 
centre. Using a 250 diameter beam of divergence 2 x 10-3 radians the 
Debye-Scherrer rings from the rolled metal were found to be continuous. 
In the central region of the sheet where the grain size was large the rings 
were incomplete and made up of arcs which again had a continuous 
intensity distribution. A fracture surface was formed in such material by 
growing a fatigue crack at 3 + } ton/sq. in., the crack growing so as to extend 
approximately 1em in 107 cycles at 2000c/m. Using the same camera 
geometry the fracture surface was examined and showed prominent 
spottiness of the Debye-Scherrer rings. Figure 1} illustrates the change in 
detail of the Debye—Scherrer ring for the camera geometry employed. By 


{ All figures are shown as plates. 
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irradiating the mid-section of the sheet and obtaining arcs instead of 
complete rings it was also evident that the building-up of sub-grains by 
cycling involved an increased spread of the ares. Thus the range of 
misorientation within a rolled grain increased as the sub-grains formed. 
The effect of the localized cyclic stressing at the crack tip had, therefore, 
been to drive the dislocations into boundaries defining sub-grains of greater 
size and misorientation than had been produced by the initial rolling 
deformation. 

The fracture surface structure was found to be localized since removal of a 
50 u layer from the fracture surface revealed material producing continuous 
Debye-Scherrer rings. This fact precluded any micro-structural study of 
the sub-grains. It was argued that the fracture surface structure might 
be simulated in a plain direct-stress specimen if the range of stress was high 
enough, since, although the fatigue crack was grown slowly, it would have 
advanced for distances of the order of the sub-grain size in a matter of 
minutes. Plain direct-stress specimens of the commercial purity Al were 
tested for some 30000 cycles at stress levels which were known to give 
endurances of the order of 100000 cycles. Using both annealed and 
unidirectionally strained initial structures the effect of cycling was found 
to be the production of well-developed sub-grains, a structure which was 
found to be common to the interior as well as the exterior of the plain 
specimens. Figure 2 illustrates, for the case of a single grain, how the 
characteristic tensile deformation, as shown in a Transmission Laue 
diagram, is built up into a sub-grain structure by cyclic deformation. 

Closer examination of the plain specimens showed that the fracture 
surface structure had not been completely simulated. Although numeric- 
ally similar, the shapes of the individual diffraction spots on the Debye— 
Scherrer rings were quite different. For the same camera geometry the 
fracture surface spots were more sharply defined, smaller and free from 
radial and tangential distortion. These features suggest that the interior 
of the sub-grains in the fracture surface are more free from residual dis- 
locations than the direct -stress specimens and that the sub-grain boundaries. 
are more sharply defined. It was found impossible to simulate the fracture 
surface structure in respect of the spot shape by any high range of direct. 
stress on a plain specimen. It was concluded that the initiation and rapid 
growth of a surface crack was occurring before the dislocations could be 
driven into boundaries to the extent which occurs at the crack tip in the 
propagation-type specimen. In order, therefore, to follow the develop- 
ment of the sub-grain structure as far as possible, torsion specimens were: 
employed. Aluminium of purity 99-99% was also used for these specimens. 
so that metallographic methods could be used to show up the sub-grain 
structure. For specimens containing a few large grains in the gauge length 
the same grain could be irradiated after increasing numbers of cycles. 
Figure 3 shows the development of structure within a grain at a plastic 
strain range of ~10~*. After a few hundred cycles the deformation of the- 
crystal has reached the stage where the spreading Laue spots have begun. 
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to show a sub-structure and the misorientation range is such that parts of 
Debye-Scherrer rings emerge. By 2000 cycles the build-up of sub-grains is 
illustrated by the increase in extent of the spread of the spotty Debye- 
Scherrer arcs. At 4000 cycles the surface of the specimen shows many small 
cracks, particularly at grain junctions, but the misorientation range of the 
sub-grain structure can be observed to increase for a further 2000 cycles 
since, under torsion conditions, the test is not terminated by the propagation 
of a crack. 

The misorientation at a sub-grain boundary is large enough to allow 
polarized-light techniques to be used to show the range of sub-grain sizes 
and their development. Figure 4(a) shows the characteristic short-slip 
lines (Seeger 1956) on a torsion specimen after a few hundred cycles at a 
plastic-strain range ~ 10~*, fig. 4 (b), shows the same field under polarized 
light. The sub-grains are ‘obaet ved to be forming in the regions between 
the slip-line traces. With increasing number of cycles many cracks appear 
in torsion specimens in the regions of well-developed sub-grain structure. 
In polycrystalline specimens the sub-grain structure is concentrated in the 
region of grain-boundary junctions due to stress-concentration effects and, 
therefore, the cracks tend to be associated with grain boundaries. In 
single-crystal specimens, however, the cracks have a surface trace with 
irregularities on the same scale as the sub-grain size. 

It is noteworthy that the cyclic sub-grain structure is also readily 
formed at —196°c. Thin slices of single-crystal specimens showed pro- 
minent double-Bragg reflection effects (Atkinson 1959) when examined in 
transmission. The cyclic structure is also resistant to recrystallization, 
an anneal of four hours at 330°C, which readily caused complete recrystal- 
lization in an aluminium specimen deformed to 15% strain in tension, failed 
to cause any new grains to appear in the cyclic sub-grain structure. A 
partial recrystallization developed after a further four hours at 420°c. 

The characteristics of the sub-grain structures found in the various 
specimens are tabulated in tables 1 and 2. 


2. Mild Steel and Iron 


_ Fracture surfaces in annealed 0-05°%C mild steel of uniform grain 

size ~20y, resulting from the growth of cracks at 5+3tons/sq. in. 
and 12+ 3 tons/sq.in. were examined with beam diameters of 75, and 
divergences between 1 and 5x 10-%radians. It was found possible to 
estimate the particle size and to measure the total angular misorientation 
within a grain. In contrast, specimens of the same mild steel pulled in 
tension to strains of up to 10% displayed quite continuous Debye—Scherrer 
rings for the same camera geometry. 

As in the case of aluminium, plain direct-stress specimens of mild steel 
were subjected to high stress levels in short-time tests in an attempt to 
simulate the fracture surface structure. Although the Debye-Scherrer 
rings showed spottiness compared with the unidirectionally strained speci- 
men measurements were not possible. Specimens of an iron were next 
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Table 1. Reflection 


Average Total Mean angle 
Specimen size of | misorientation] between | Broadening 
pile sub-grain range sub-grains dd 
phi (1) B a=B/3 
Fracture surface : 4 
AI99-8%3+44 4-6 > 40° >13 ~1x10- 
ton/sq. in. 
| Fracture surface 
Al199-8% 341 4-6 > 40° > 13° ed We Le 
ton/sq. in. 
| Torsion specimen, 
A199-9°%, 6000 4-8 > 40° 15, ~5~x 10-4 


cycles at ep) =10~ 


| Torsion specimen, 
Al 99-9°% 80 000 
cycles at 
€pl= 3 x 10-3 


4-8 > 40° Saal ey ~5~x10-4 


Direct stress, 


A199-9%, 30 000 25°-36° g°_12° 


eycles at 
ép1=8 x 10-3 


Table 2. Transmission 


Direct stress, 
A199-9°% 30 000 
cycles at 

Epl = 8 x 10-3 


5-8 30°—40° 12° Lx ig= 


Table 3. Reflection 


Average Total Mean angle 
eas size of misorientation between 
P sub-grain range sub-grains 

vr/8 (12) B a= /3 

Fracture surface mild steel Some eee 
5 +24 tons/sq. in. Sn eae Baas 

| Fracture surface mild steel é aye Ee 
12+8 tons/sq. in. ae rie a oes 
Direct stress, iron 30 000 S F ghee Pave 

cycles at epi ~ 10-3 ie 20°28 thes 


y 


Pa eee 
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used in direct stress and these proved to have a sub-structure upon which 
measurement was possible. The measurements are summarized in table 3. 

Sheets containing growing fatigue cracks were examined metallo- 
graphically with the plane of section normal to the mean direction of crack 
growth. The crack front was always found to be slightly bowed forward so 
that sections could be made with the leading part of the crack in the middle 
of the sheet. From such sections, taken at intervals of a few microns 
right through the position of the front, it could be seen that the trace of the 
crack at completely broken sections owed its irregularity to the fact that it 
was formed from the joining up of many independent cracks formed ahead 
of the completely broken section at any stage of propagation. Figure 5 (a) 
shows detail of one end of a bowed front at a section where the outside 
faces of the sheet are unbroken down to a depth of 20% of the sheet thick- 
ness. Figure 5 (5) shows a final section of the leading cracks in the middle 
of a sheet; removal of a further 3 at this section eliminated all traces of 
cracks. The irregularities in the small independent cracks are on the scale 
of the observed sub-grain structure. In the case of the mild-steel speci- 
mens it was notable that these elementary cracks were not found to bear 
any relation to the known form of elementary brittle cracks in the mild 
steel or to the particular inclusions which were known to play a part in the 
ductile fracture of the material. 


3.3. Other Metals 


With the micro-beam geometry developed for mild steel attempts were 
made to detect sub-grain structures on fracture surfaces of other metals. 
These were unsuccessful in the cases of copper, stainless steel and two alum- 
inium alloys, one of duralumin type, L.71, and the other a Zn, Mg alloy 
DTD.687A. In the case of a Zn sheet of 99°% purity, well-developed sub- 
grains were detected and measurements would have been possible with the 
75 beam. In the case of copper torsion specimens a local sub-structure 
associated with twin boundaries was observed metallographically after 
x-ray observations had indicated the existence of sub-grains. 


§ 4, DIscuUSSION 


The structure detected by the micro-beam method can be related to some 
long-standing observations on the effect of cyclic stressing. ‘The sharpening 
of Debye-Scherrer rings from iron deformed by cycling was observed by 
Hempel and Moller (1954) and was considered to be a result of the reduction 
of strains between crystals. The formation of sub-grains in annealed 
super-pure aluminium and in iron was observed metallographically by 
Forsyth (1952) who considered the effect to be most pronounced at high 
ranges of alternating strain and interpreted the structure as due to poly- 
gonization consequent upon the production of vacancies in the fatigue 
test. The micro-beam observations show that the cyclic sub-structure 
develops in previously deformed crystals ‘and is, therefore, not simply a 
breaking-down structure but a structure which is formed by the dislocations 
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reaching a preferred distribution under cycles of large alternating strain. 
This distribution develops at —196°c so that the dislocations are not 
dependent upon diffusion mechanisms to form the sub-grain arrays, but 
must reach their positions by cross-slip processes. 

The cyclic structure is shown to reach an extreme condition, i.e. large 
misorientation between particles, «, and sharply defined diffraction spots, 
in material subject to the stress concentration that exists at the tip of a 
slow growing crack. The sub-grain formation is the same for the specimens 
from fracture faces irrespective of the range of alternating stress used to 
produce the fracture. The sub-grain structure appears, therefore, to be 
characteristic of the fracture condition at the crack tip and not of any 
intervening state of plastic deformation. We would expect, therefore, that 
a dislocation model which complies with the data determined by the 
micro-beam technique would contain features suggestive of the mode of 
fracture at the slowly growing crack. 

Following Hirsch (1952) the sub-grain structure can be described in 
terms of the density of excess dislocations of one sign in the boundaries, D,, 
and the residual density of dislocations left within the particles D;. In this 
way the cyclic structure may be compared with the structure found in 
rolled super-pure aluminium. Table 4 gives the comparative data. 


Table 4 


Density of dislocations 
Mean angle 
between 
sub-grains = im 
a 


Total 
misorienta- 
Specimen : tion 
range 


boundary | sub-grain 
D ; 


Super-pure Al 
rolled to 57% 
original thick- . Ox 1U8 
ness Hirsch 
(1952) 


A199:8% 
fracture surface SR 2-3 x 1019; ~1x107 
3+4 ton/sq. in. 


In the case of rolled aluminium, Hirsch considered that the resultant 
sub-grain structure was best described by a polygonized arrangement of 
dislocations such that «=f .v"3/7’ and D,=8/bT where T was the original 
grain size and b the Burgers vector. The cyclic sub-structure can, how- 
ever, be observed with polarized-light techniques known to be insensitive 
to small misorientations, so that the data for this structure are interpreted 
in table 4 on the basis that the sub-grains are not polygonized but mis- 
oriented about a mean position such that w=8/3 and D,=/8b.vt8 
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(Hirsch 1952). The value D, is estimated in the same way for the two 
deformation structures on the basis that there is a random distribution of 
dislocations within the sub-grains (Hirsch 1952). 

The values of D, and D, in table 4 show that the main feature of the 
cyclic sub-structure is the accumulation of high dislocation densities in 
sub-grain boundaries brought about by the removal of dislocations from 
the sub-grain interior. The cyclic sub-grain structure would not be 
expected to be a strong source of hardening (Ball 1957) of the metal, but 
the misorientation of regions of crystal through such large angles, «, would 
cause stresses to develop in adjacent parts of a grain and the densely packed 
sub-boundary regions might be expected to be sites of micro-cracks. Such 
micro-cracks would not be expected to propagate immediately after forma- 
tion since the material within the sub-grains would be still quite ductile. 
The first micro-cracks would serve to cause further local stress concentra- 
tion and aid the formation of cyclic sub-structure at the region of the main 
crack tip. This is in accord with two features connected with slowly 
growing fatigue cracks, that no large micro-hardness changes can be 
detected around the cracks and that the crack front consists of many 
small elementary cracks unconnected with the metallurgical features 
known to be effective in brittle and ductile fracture. 


Table 5 


Thornton and Hirsch 
(1958) Relative rate of crack 
Seeger (1956) propagation Frost (1959) 
if 


Metal 


Stainless steel, 18/8 13 ergs/cm? 0-7 
Mild steel — EO 
Copper 40 4-0 
Aluminium 200 12-0 
Zine Large 150-0 


If the dislocations are considered to perform a symmetric random walk 
before becoming trapped in sub-boundary arrays the density of disloca- 
tions in the arrays would increase as N14 multiplied by the number of 
dislocations which were moving cyclically. Under cyclic conditions of 
constant plastic strain range, <p1, the number of moving dislocations would 
be proportional to €p, and cracks would appear when ep). N1? reached a 
large enough value as is experimentally observed in some metals (Coffin 
and Tavernelli 1959). The ease with which the sub-grain structure could 
form would depend upon the ability of the dislocations to cross-slip and 
thus upon the stacking-fault energy, y. It is, therefore, of interest to 
tabulate the stacking-fault energy and the relative rate of crack propaga- 
_ tion (at unit alternating stress and same crack length) for various metals. 
| The hypothesis can therefore be made that, under the cyclic conditions 
_ prevailing at a slowly growing crack tip, the dislocations tend to form 
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regions of high and low density. The resultant misorientations within 
a grain are detected as well-defined sub-grain boundaries in metals of 
high stacking-fault energy; this ease of forming sub-grains favours the 
formation of the micro-cracks which are an essential stage in the pro- 
gression of the fatigue crack. On this basis some further features of 
crack propagation in sheets can be understood. Thus, if the micro- 
cracks formed in the densely-packed boundaries are unable to propagate 
but serve only to advance the cyclic, sub-grain structure formation, the 
progress of the main crack front will depend on the formation of a sufficient 
number of independent micro-cracks until the applied maximum tensile 
stress pulls off the remaining necks of material between the micro-cracks. 
This pulling-off effect is minimized on the crack surfaces examined in 
the present work since the maximum tensile stress of the cycle was small 
and the crack length was only a small proportion of the sheet width in 
which the crack was growing. It has been observed, however (Frost 
1959), that, as the crack increases in length beyond this condition, the 
fracture surface roughens and, in thin sheets, goes over from a position 
normal to the loading axis to an inclined position. The onset of the 
transition correlates with the maximum tensile stress of the cycle and 
would be interpreted as being due to the gradually increasing effectiveness 
of the pulling-off stage which contributes to the mechanism of advance 
of the main crack front. 

The further main feature of the model of growth is that, since the 
first micro-cracks are dependent upon the formation of high dislocation 
density boundaries defining particles in the range of size of several 
microns, it is clear that an adequate volume of material at the crack 
tip must be subject to alternating cyclic strain. If this is not so a 
sufficient number of sub-grains with their associated stress effects due 
to their large misorientations will not be formed. The range of alternating 
stress, o, and the crack length, J, determine the volume of the region 
subject to alternating strain and these quantities appear in the crack 
growth parameter o*/ in eqns. (1). The process of crack growth as opposed 
to surface crack initiation can, therefore, be defined as occurring when: 


(i) a large-enough volume of material is subject to cyclic plastic strain 
sufficient to cause the development of cyclic sub-grain structure 
to the extent that micro-cracks form at the sub-boundaries; and 


(ii) a sufficient normal stress is present to pull off the necks of material 
which remain between sub-boundary micro-cracks. 


The many processes whereby surface cracks can be formed are regarded, 
therefore, as controlling the fatigue process only until the surface crack 
is long enough to satisfy condition (i) for the stress system on the specimen. 
The length of the surface crack before crack growth as defined commences 
will be much larger, for example, for copper than for aluminium on 
account of the large difference in ease of development of cyclic sub -grain 
structure due to the difference in stacking-fault energy. The importance 
of chemical attack by the atmosphere will also be greater for materials 
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which have a longer period of their fatigue-life taken up in forming a 
crack long enough to grow by the sub-boundary crack mechanism. Ks 
regards diffusion effects, the cyclic sub-structure will resemble a erain- 
boundary system, with the condition that the small size of the sub -orains 
will extend the temperature range over which grain boundaries are 
usually regarded as being effective in contributing to bulk diffusion. 

The failure of surface cracks to propagate is met with experimentally 
when either the cycle of deformation is wholly compressive or the value 
of the parameter o°/ falls below a critical value. Condition (ii) is not 
fulfilled for the compressive cycle and it is suggested that the critical 
value of o%/ for propagation is due to the need to fulfil condition (i). 
Measurements (Frost 1959) on mild steel of the critical value of oJ lead 
to a value of /~ 50, for values of o equal to the plain fatigue limit. The 
processes which form surface cracks on plain specimens, e.g. the develop- 
ment of surface notches as calculated by May, would have, therefore, to 
cause a 50 surface crack before growth by the sub-grain mechanism 
would take over for a range of applied stress equal to the fatigue limit. 
The scatter associated: with fatigue tests on plain specimens at long 
endurances is associated in the hypothesis with the formation of surface 
cracks rather than with their growth. In this connection it may be 
noted that, in experiments in which steps are taken to initiate the surface 
crack by an additional process, e.g. as in fretting, the scatter in the fatigue 
test tends to be diminished. 


§ 5. CONCLUSIONS 


1. Under large ranges of alternating plastic strain the dislocations in 
aluminium and «-iron form sub-grain boundaries irrespective of the 
initial state of distribution. 

2. The cyclic sub-grain structure tends to a characteristic limiting 
condition of sub-grain size, misorientation and dislocation density within 
the sub-grains with increasing number of cycles. The limit is reached 
when micro-cracks form at the sub-grain boundaries. The most developed 
sub-grain structure, i.e. largest misorientation and lowest density of 
dislocations within the sub-grains, is observed on the fracture surfaces 
formed by slow fatigue cracks. 


3. The commencement of fatigue crack growth is, therefore, associated 
with the fracture process based upon the formation of sub-grain cracks, 
a phenomenon not connected with the metal surface. 

4. Two aspects of fatigue crack growth can be interpreted on this 
model. (a) The order of susceptibility to fatigue-crack growth in metals 
relates to the ease of formation of sub-grains and hence to the stacking - 
fault energy; (b) the phenomenon of non-propagation of surface cracks 
is expected when the volume of material ahead of the surface crack which 
is subject to the range of alternating strain is too small to allow the 
characteristic sub-grain structure to develop. 
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ABSTRACT 


The half-life of the naturally occurring odd—odd isotope ®°V for electron 
capture decay to the first excited state of ®°Ti is shown to exceed 8 x 10% 
years and for negatron decay to the first excited state of *°Cr to exceed 
1-2 x 101° years. 


§ 1. [InTRODUCTION 


Ty 1952 Johnson showed by mass measurements that the odd—odd isotope 
°°V was unstable with respect to beta decay to its even—even neighbours, to 
°Ti by (2-39 + 0-13) Mev and to Cr by (1:18+0-12)mev. Early attempts 
to measure the half-life merely agreed that it was long (Sheline 1954, 
Strome 1954, Selig 1954, Heintze 1955). Ti, however, has a first excited 
state, presumably 2+, at 1-58 Mev (Pieper 1952). °°V has a ground state 
spin of 6 units (Kikuchi ef al. 1952, 1953, Baker et al. 1952) and presumably 
even parity, since the odd neutron and proton are likely both to be in fy. 
configurations. Hence the fastest mode of decay in the electron capture 
branch is expected to be the fourth forbidden (4, no) transition to the 
1-58 Mev state of °°Ti since the ground to ground transition involves two 
more units of spin change. Recent attempts to measure the half-life for this 
mode of decay (by finding the rate of emission of 1-58 Mev gamma rays from 
vanadium) agreed that the half-life is around 4 to 5 x 10'4 years (Glover and 
Watt 1957, Bauminger and Cohen 1958). However, unless the log ft 
value for a fourth forbidden beta transition is considerably smaller than 
would be expected from the trend of values shown by first, second and third 
forbidden transitions, this half-life seems to be one or two orders of magni- 
tude too short. It seemed worth investigating the decay again to see if 
this discrepancy from the expected log ft value for a fourth forbidden 
transition were real. 

The recent discovery (Porter ef al. 1958) of an excited state in °°Cr at 
0-78 Mev which fits in well with the general systematics of the first excited 
2+ states of even—even nuclei permits us to examine another fourth for- 
bidden transition mode from ®°V, that by negatron decay to °°Cr via this 
excited state, by searching for 780 kev gamma rays from °°V. A search for 
both 1-58 Mev and 780 kev gamma rays from vanadium in order to arrive 
at a value for the half-life of °°V is described below. 


+ Communicated by Dr. 8. C. Curran. 
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§ 2. APPARATUS 


A 4lin. x 44in. cylinder of Nal(Tl) was used as the detector and the 
pulses were amplified, and recorded on a 100-channel kicksorter. The 
source, 500g of V,0; was compacted in a thin copper holder which sat 
snugly over the crystal, so that the source covered the crystal top and side. 


§ 3. METHOD 


Since the counting rates from the two gamma rays were expected to be 
extremely low, the experiment was designed so that the results would be as 
significant as possible by particular attention to three points: Reduction 
of detector background by massive iron and lead shielding and by geiger 
anti-coincidence techniques, good source-detector geometry combined with 
large quantities of source, and precautions in the method of analysis against 
long term gain drifts. This last precaution consisted of taking short 
~ 2-3 hours runs of background and background plus source alternately, 
with frequent energy calibration of the system. In each set the two 
spectra were subtracted and the count rates within a gamma ray peak 
width about 1-58 Mev and 780 kev were determined and recorded. A large 
number of separate such runs gradually built up useful integrated count 
rates, and in this way gain drifts could be counteracted. 


§ 4. Resutts anp DiscusstIon 


Normally one would simulate the source in such an experiment by 
suitable inactive materials for use in the background measurements. 
Slightly variable °K contamination proved to be troublesome, however, so 
the background was taken with the crystal alone. The increased gamma 
ray interaction in the source material was estimated by measuring count 
rates in the difference spectra at energies above and below the expected 
positions of the two *°°V gamma-ray energies. The trend shown by these 
counts was taken to indicate the true background and the count rates at 
1-58Mev and 780 kev related to this value. Corrections were applied to 
allow for the slight 4°K contamination of the source. The small extra 
uncertainties introduced by the corrections have been allowed for mathe- 
matically in calculation of the statistical error in the results. 

In fact, several separate experiments involving a total counting time of 
about 500 hours agreed in showing that the emission of gamma rays from 
the vanadium was undetectable and that the count rate in the 1-58 Mev and 
780 kev regions of the spectrum did not rise above background. The 
efficiency of the system was calibrated using potassium (assuming a gamma 
emission rate of 3-36 +0-10gammas sec-!g~! of natural potassium) and 
*Na. ‘The final results, allowing for crystal efficiency, were 


(—2+29) 1-58Mev gammas per 100g natural vanadium per hour, 
(+6+14) 780kev gammas per 100g natural vanadium per hour. 


The statistical limits show that the half-life of 5°V for electron capture 
decay to the first excited state of °Ti must exceed 8 x 10! years, and for 
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negatron decay to the first excited state of ®°Cr must be greater than 
1:-2x 101 years. These results give log ft values >22-1 and > 22:5 
respectively which approach more nearly to what would be expected fora 
fourth forbidden (4, no) beta transition. 


§ 5. CONCLUSION 


As expected from consideration of the log ft value, the half-life of 5°V 
for electron capture decay to the first excited state of °°Ti has been shown 
to be much larger than previous measurements of 4 to 4:8 x 10!4 years and is 
not less than 8 x 101° years. The half life for negatron decay to the first 
excited state of *°Cr is similarly shown to exceed 1-2 x 1016 years. 
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ABSTRACT 
A 47 proportional counter system capable of examining thin sources of 
extended area has been used to determine the half-life of 87Rb, which was 
found to be (5:25+0-10)x10! years. Corrections for absorption of 


electrons and for scattering in the source and in the source supporting foil 
are discussed. 


§ 1. INTRODUCTION 


REcENT advances in rubidium-strontium geochronology have shown 
that the method can give self-consistent ages which are more or less in agree- 
ment with ages derived from uranium—lead and potassium—argon analyses 
when the half-life of *’*Rb is assumed to be near 5 x 101° years (Aldrich and 
Wetherill 1958). Lately, differences in the ages derived by the Rb-Sr 
and K—A methods have been shown to indicate the more recent thermal 
history of an ancient rock system, since heating of the rock by late volcanic 
intrusions will be expected to reduce the amount of argon retained in the 
mineral whilst leaving the *’Rb/®**Sr ratio relatively unchanged (Giletti 
1959). The increasing use of the *’Rb—*’Sr method in dating studies has 
made careful evaluation of the half-life of °’Rb a matter of some importance, 
because the accuracy of a measured age depends on the accuracy to which 
the half-life of the transition is known. Independent physical measure- 
ments of the half-life are therefore necessary to check the value of 5 x 101° 
years usually accepted by geochronologists on the basis of correspondence 
between ages determined by the different methods. 

The 8’Rb-—>8’Sr decay is known to be a third-forbidden (3, yes) beta 
transition with a highly forbidden spectrum shape which shows no maximum 
but appears to rise steadily as the energy decreases (Flynn and Glendenin 
1959, MacGregor and Wiedenbeck 1954, Lewis 1952, Curran et al. 1951). 
Although the energy of the transition is about 275 kev the shape of the spec- 
trum is such that the average energy of the beta particles is low, about 45kev, 
and care must be taken to reduce and correct for losses caused by absorption 
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of the softest electrons in the source or source mount. Values of the half- 
life ranging from 4-3 x 101° to 6-4 x 10'° years have been reported, and are 
summarized by Aldrich and Wetherill (1958). The most recent experi- 
ments by Huster and Rausch (1956) using thin sources in a 47 counter gave 
4-95 x 10!° years, but there seemed still to be some doubt about this value 
and later work by Huster (1957), using absorbers and extrapolating to zero 
thickness, gave a half-life of 4-6 x 101° years. In experiments using scin- 
tillation counting Flynn and Glendenin (1959) derived a half-life of 
(4-7+0-1)x 10! years. The only earlier experiment in which source and 
support absorption should not have raised any problem was that of Lewis 
(1952) who used a RbI scintillation spectrometer and measured a half-life 
of 5-9 x 10!° years. Insensitivity to low energy electrons would, however, 
tend to make this a high value. 

In the present work we aimed at reducing absorption of the softest beta 
particles by using very thin uniform sources (produced by evaporation in 
vacuo from a tantalum strip) deposited on thin supporting films. The use 
of enriched rubidium (> 99% 8?Rb) and large source areas helped to improve 
the counting statistics for thin layers. Further improvements in statistics 
were made by shielding the counter in 8 in. of iron and operating it in anti- 
coincidence with the output from a surrounding ring of geiger counters to 
reduce background counts. Internal background was kept to a minimum 
by lining the counter with carefully cleaned copper sheet. The amount of 
87Rb present in the sources was determined to high accuracy by isotopic 
dilution and mass spectrometer analysis. Counting in 47 geometry was 
congidered essential in order to keep back-scattering corrections within 
manageable limits, but in addition to the 47 count we simultaneously 
recorded separate counts from each side of the source and used these counts 
to estimate the effect of absorption and back-scattering within the suppor- 
ting film. 


§ 2. Basic ASSUMPTIONS AND THEORY 


The detector which was used in our study of 8’Rb was a large 47 propor- 
tional counter tube. A solid source gas counter was chosen in preference 
to a scintillation detector because counting could be extended to much 
lower energies. Of gas counters, the proportional was used in preference 
to the geiger counter, since proportional counters appear to be more tolerant 
of changes in the nature of the cathode surface than do geigers and also, in 
our experience, are much more stable over long counting periods. 

The source was mounted on one side of a foil and three simultaneous 
counts taken, as in fig. 1(a). N, is the count rate, extrapolated to zero 
energy, on the side of the counter facing the source. N, is the count rate 
observed at the back of the supporting foil and J, is the count rate obtained - 
when the outputs of the two sides of the counter are added together. It 
follows that electrons scattered from one side of the counter, through the 
foil and source, to the other side count only once in the V, channel and once 
ineachofN,andN,. The block diagram of the system is shown in fig. 1 (0). 
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Suppose now that N electrons actually leave the source and that the 
source support transmits a fraction a of the electrons entering it. (N is not 
the absolute decay rate but is somewhat less due to source absorption.) Of 
the electrons entering the foil and not transmitted let a fraction b be back- 
scattered, and let a fraction s of the electrons counting in either side of the 
counter be scattered backwards in the direction of the other side (cross- 
scattered electrons). Then 


Nie +b(1—a) a + sa? ay 


a 
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N N 
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Counter system. 


We allow for the effect of source absorption on cross-scattered electrons by 
taking measurements with different thicknesses of source as described later 
and extrapolating to zero source thickness. Double cross-scattering or 
scatter in the gas followed by back-scattering from the support does not 
affect the argument, since these events do not increase the count rate, 
because of the integrating action of the counter. 

We are, of course, assuming that degradation of the energy spectrum of 
the electrons in traversing matter does not seriously interfere with the 
value of a,s orb. As we shall show later, we can arrive at essentially the 
same answer for the half-life of *7Rb by a completely independent path | 
without making any assumptions about absorption or scattering, so that 
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our postulates would appear valid. It may be that the shape of the low 
energy spectrum of §’Rb does not alter much on further degradation by 
passage through matter. 

From the above equations we can derive an expression for NV, the number 
of beta particles escaping from the source, 


2N N,b+N 
= es N= 
yoy | ‘ 1+6 |i 


Hence, since we measure N,, VN, and NV, we can find, V, provided we know db. 
If we neglect back-scatter from the support, i.e. b=0, we get 
2N,(N.— Ne) 

2N,—(N, +2) 

This is always larger than when we take back-scatter into account and leads 
to a lower value for the half-life. This is because back-scattering from the 
support tends to increase NV, relative to V, and if we ignore it, the absorption 
of beta rays in the support appears to be worse than it really is, and we over 
correct. 

We have not attempted to measure b at low energies but instead have 
taken published values for the saturation back-scattering factor of electrons 
into 27 geometry since, for electrons which are not transmitted, the support 
appears to be of saturation thickness. 

The values taken were b=0-33 for aluminium (McNair et al. 1956) and 
6=0-5 for the VYNS plus gold supports. This latter figure is based on the 
assumed value for aluminium, and on the way in which the saturation 
back-scattering factor varies with atomic number (Nijgh et al. 1959), 
taking into account the relative thicknesses of gold and plastic. The 
values assumed for the back-scattering factors do not appear to be 
critical as a 10% change in back-scattering factor produces less than 0:5% 
change in the final half-life. 


§ 3. APPARATUS 


The counter is shown in outline sketch in fig. 1(6). It is a copper-lined 
cylinder, 14 cm in diameter divided along a diametral plane by a brass plate 
carrying a sliding shelf which can easily be removed to insert the source 
supporting foils. The effective counting length is 34cm and the source 
area is some 50 cm? or 60 cm?, depending on the trays used. The foil area is 
generally some 30% greater than the source area. The two halves of the 
counter operate independently. The thin tungsten wire anodes lie parallel 
to the axis of the cylinder in the diametral plane perpendicular to the plane of 
the source and some 3cm away from the source on either side. Initially 
we were rather concerned with the effect on the counting characteristics of 
the D-shaped cross section of each counter and placed two more wires in 
each counter in the corners of the D. These wires could be run above or 
NE 


} Dr. D. B. Smith of the Wantage Radiation Laboratory informs us that 
he has deduced a similar expression to allow for absorption in 47 counting. 
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below cathode potential. However, the characteristics of the counters 
appeared to be independent of the voltage on these wires over a wide range, 
and they were removed to keep the counter as simple as possible. 

Outputs from each anode pass, as in fig. 1 (b), to separate amplifiers which 
have been modified so that no after-pulsing takes place following severe 
over-load pulses. By so modifying the amplifier the counter can be run at 
high gain and the limiting factor becomes counter rather than amplifier 
noise. Careful earthing and screening arrangements effectively removed 
external interference and we can now switch off the counter voltage and 
run the apparatus at its most sensitive setting for several days at a time and 
detect no counts at all from external effects or amplifier noise. 

After amplification the outputs are passed through discriminators to 
scaling units. The two amplifier outputs are also added together in a 
separate channel and passed to a third discriminator and scaler to give the 
sum count \,. In later experiments the discriminators and scalers 
were replaced by three 100-channel kicksorters. These give equivalent 
information but permit better counting statistics to be obtained more 
quickly by counting at all discriminator levels simultaneously. 

A geiger anti-coincidence ring placed round the proportional counter 
operates an over-riding gate on the system, blocking all three pulse channels 
each time a cosmic ray particle triggers the geiger array. The counter 
and anti-coincidence ring are shielded by 8in. of steel. Backgrounds of 
around 36c¢.p.m. are obtained in each counter with about 55c¢.p.m. in the 
sum channel. The additional count rate obtained in the sum channel from 
a typical 5 ug/cm? layer of ®’RbCI is of the order of 20 ¢.p.m. 


§ 4. MANUFACTURE OF SOURCE Mounts AND SOURCES 


In the main experiment the source supporting foils were made from VYNS 
by the flotation technique (Pate and Jaffe 1955) and coated on both sides 
with about 15 g/cm? of gold, by vacuum distillation. The thickness of the 
foils was found to be fairly constant and averaged about 42ug/cm?. We 
tried to reduce the foil thickness by coating with half the amount of gold, 
but found that we could not get consistent count rates from one run to the 
next, possibly due to poor conduction in the gold layer. We could produce 
conducting layers of aluminium down to about 2ug/cm® but the VYNS 
invariably became too brittle to use in the counter. 

The source layer was also deposited by vacuum distillation in order to 
obtain uniform thin layers. Only vacuum distillation can give the uniform 
sources with minimum absorption which are essential to an experiment 
of this kind. This has been amply demonstrated by Merritt et al. (1959) 
and confirmed by preliminary work of our own. In these preliminary 
experiments, layers of RbCl of 50 g/cm? and 80 g/cm? were distilled in 

vacuo onto the foil and counted. The layers were then dissolved in distilled 
water in situ and allowed to dry in as uniform a layer as possible. The 
count rate was then some 30 to 40% less than before. On re-dissolving the 


568 A. McNair and H. W. Wilson on the 


layers in distilled water, to which a little wetting agent had been added, the 
count rate rose again but was still 15 to 20% below the original. 

Pure natural rubidium chloride was used for the sources thicker than 
12 ug/em? and enriched rubidium (99-3% *’Rb) for sources thinner than 
12 ug/em?. 


§ 5. ESTIMATION OF THE SOURCE 


The amount of ®’Rb in the source was estimated by the mass spectro- 
meter isotope dilution technique. The entire source and its supporting 
film were carefully removed from the counter shelf and placed in a crucible 
containing, in solution, an accurately measured quantity of rubidium of 
different but known isotopic ratio. The change in isotopic ratio gave the 
number of atoms of ®’Rb in the source. 


§ 6. PROCEDURE 


The procedure, then, is to make a series of VYNS films, coat them with 
gold, and measure the background of the counter with the films in position. 
In the earlier experiments counts were taken at five bias points in each pulse 
channel from an energy corresponding to about 3 kev down to about 300 ev. 
The integral bias curve was then extrapolated to zero energy. The extra- 
polation generally amounted to between one and two per cent of the count 
rate. In the later kicksorter experiments the extrapolation was between 
3 and 4%, but both counting methods were shown to give the same results 
within statistical error. The same films are then coated with source, and 
the counts repeated. The difference counts in the three channels give N,, 
N, and N,, and the number of electrons emitted from the source can be 
calculated according to the formula already derived. Thereafter the 
amount of §’Rb in the source is measured by the isotope dilution technique, 
and the half-life calculated in the normal way. 


§ 7. RESULTS 
7.1. Source Absorption 


We have corrected for self absorption in the source by making half-life 
measurements on a series of sources of different thicknesses up to about 
704g/em?. ‘The half-lives so obtained are plotted against source thickness 
asin fig. 2. We find that, within experimental error, the apparent half-life 
is a linear function of source thickness up to 70.g/em?. Extrapolation to 
zero source thickness is therefore simple, and yields a half-life for infinitesi- 
mal source thicknesses and hence, presumably, zero source absorption. 


7.2. Source Support Absorption and Back-scattering 
We correct for this in fig. 2 by applying the formula for V already derived. 
However, to test the effectiveness of this correction, we made a series of 
measurements using a supporting foil of different material (thus with 
a different back-scattering factor 6) and different thickness (hence different 
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electron absorption). The most convenient material was commercial 
aluminium leaf of mean thickness 180g/em?. The half-lives measured 
with this support, and corrected by the previous formula, are also shown in 
fig. 2 and prove to be collinear with those obtained using the thin support. 
The best straight line through these points, then yields our value of the half 
life of *’Rb corrected for source and supporting foil absorption and foil 
back-scattering. This turns out to be (5:37 + 0-02) x 10! years. 


Fig. 2 
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Apparent half-lives of 8’Rb corrected for back-scattering and absorption of 
electrons in the source support. 


7.3. Maximum Possible Half-life 


If we calculate half-lives directly from the apparent 47 count rate we 
shall obviously get a maximum possible half-life, since we are ignoring 
support absorption. These half-lives are shown in fig. 3, once more 
plotted against source thickness. Once again the relationship is a linear 
one, and by extrapolating to zero we remove the effects of source self- 
absorption. In this case, however, we get two separate parallel lines and 
the values measured from sources deposited on aluminium supports show 
the expected extra absorption over those mounted on thin VYNS plus gold 
supports. 

We find 

T max = (5°79 + 0-03) x 10! years (aluminium support), 
T max = (5°53 + 0:01) x 10! years (VY NS + gold support). 


P.M. 
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We now have an interesting new possibility for arriving at the half-life 
of 87Rb without making any assumptions other than that the maximum 
possible half-life is a linear function of support thickness up to 180 pg/em?. 
This does not appear to be too sweeping an assumption since we know from 
fig. 2 that the half-life is a linear function of source thickness up to 70 peg/em 
If we make this assumption, and extrapolate linearly to zero support thick - 
ness, we arrive at a figure for the half-life which should be the true value. 
We get 5:45 x 10° years, which is in reasonable agreement with the value 
derived from fig. 2, especially as the aluminium leaf was very non-uniform, 
180 g/cm? being only an average value for the thickness. 


Fig. 3 
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To the measured value of 5-37 x 101° years corrections must be made to 
allow for the effects of dead-time in the registers and for counter geometry 
which may be less than the 47 steradians assumed up till now. The geo- 
metry factor was estimated in a subsidiary beta-gamma coincidence 
experiment using *°°Hg which has a fairly soft beta spectrum (maximum 
energy 208 key, rising steeply at low energies). The corrections for geo- 
metry and dead-time amounted to 2% of the experimentally determined 
half-life. 

The error of + 0-02 x 10!° years which has previously been quoted is 
calculated on the basis of counting statistics only, but taking into account 
all possible sources of uncertainty we believe that an error of + 0-10 x 101° 
years is more realistic. We therefore conclude that the half-life of 87Rb is 


(5-25 + 0-10) x 10" years, 


Half-life of Rubidium-87 571 


§ 8. Discussion 


Our value for the half-life of 8’Rb is rather higher than the most recent 
published figures (Flynn and Glendenin 1959, Huster and Rausch 1956, 
Huster 1957). We believe our value to be genuine but, if it is not, then 
about the only possible explanation would appear to be losses in the 
proportional counter caused by quenching of the discharge by gaseous 
impurities from the counter itself. However, we believe this to be a 
negligible effect because a plot of total counts against time is strictly linear 
for periods up to 24 hours, after which slight deviations from linearity 
become just perceptible. Any one counting period was limited, therefore, 
to 15 or 16 hours and the counter was refilled before the next run. Secondly, 


Fig. 4 
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Plot of number of counts observed in sum channel against number of atoms 
of 87Rb in the source. (i) The experimental values have been adjusted 
to what would have been obtained had all the rubidium in the source 
been 8’Rb. (ii) The size of the rectangle represents the statistical error 
in the measurements. 


the counter was operated as a continuous flow counter for about half the 
measurements with complete agreement in the results. Losses in the 
external pulse registers appear unlikely, because of the correspondence 
between results from a discriminator and scaler and those from a 100-channel 
kicksorter. Nor can our high value for the half-life be caused by inactive 
impurities in the source, since the method of estimating the source measures 
only the ®*Rb and the *’Rb content. We also find complete agreement 
between values obtained from sources of enriched §7RbCl and those from 
sources of the natural material. 

In some of our earlier experiments we found an appreciably higher count 
rate (and hence lower half-life), than in others. This effect was traced to 
thorium impurity in the tantalum metal used for the evaporating boats. On 
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heating the boat, a small quantity of ThB (?#2Pb), rather less than 104 
atoms, was deposited on to our source holders. We got round this difficulty 
by heating the boat to a high temperature before using it to evaporate gold 
or RbCl, and we checked for freedom from short-lived impurities by counting 
background and background + source over at least four separate periods of 
15 hours each. 

An attempt was made to see if there was likely to be any zero point error 
in the experiment, i.e. a sudden increase in count rate for very thin sources, 
below our thinnest of just under 5yg/cm?. To do so, we have plotted, in 
fig. 4, the sum count rate obtained (using the VY NS supports) as a function 
of number of atoms of 87Rb in the souree. In order to allow for source 
thickness, when natural RbCl was used we have scaled up the measured 
sum count to what would have been obtained if the source had been 100% 
87RbCl, and scaled up the number of atoms of ®’Rb accordingly. The 
relationship is linear and extrapolates to pass through the origin as would 
be expected if there were no zero error. 


ACKNOWLEDGMENTS 


We wish to thank the Director of the Establishment for permission to 
publish this paper. 

We also thank Mrs. R. C. Jeffery for making the isotopic ratio measure- 
ments, and Mr. J. P. Franklin for assistance with the experimental work. 

Thanks are also due to the Electromagnetic Separation Group, A.E.R.E., 
Harwell, for the supply of enriched ®’Rb. 


REFERENCES 


Aupricu, L. T., and WETHERILL, G. W., 1958, Ann. Rev. nucl. Sci., 8, 257. 

‘ Curran, 8. C., Drxon, D., and Witson, H. W., 1951, Phys. Rev., 84, 151. 

Fiynn, K. F., and Guenpenty, L. E., 1959, Phys. Rev., 116, 744. 

GiteTti, B. J., 1959, Nature, Lond.; 184, 1793. 

Huster, E., and Rauscu, W., 1956, reported by Aldrich and Wetherill 1958. 

Huster, E., 1957, reported by Aldrich and Wetherill 1958. 

Lewis, G. M., 1952, Phil. Mag., 48, 1070. 

MacGrueor, M. H., and Winprensecr, M. L., 1954, Phys. Rev., 94, 138. 

MoNarr, A., GLover, R. N., and Wiison, H. W., 1956, Phil. Mag., 1, 199. 

Pega nae Taytor, J. G. V., and Campion, P. J., 1959, Canad. J. Chem., 

, 1109. 

Niseu, G. J., Wapstra, A. H., and van Liesnovt, R., 1959, Nuclear Spectroscopy 
Tables (Amsterdam : North-Holland Publishing Co.), p. 41. 

Pats, B. D., and Yarrs, L., 1955, Canad. J. Chem., 33, 15. 


[ees7380 | 


The Solar Diurnal Variation of Cosmic Rays during 1958 and 1959, 
at Makerere, Hermanus and Herstmonceux+ 


By D. M. THomson 
Physics Department, Makerere College, Kampala, Uganda 


[Received November 21, 1960] 


ABSTRACT 


The results of observations of the solar diurnal variation in cosmic radiation 
at Makerere (East Africa), Hermanus (South Africa), and Herstmonceux 
(England) are presented for the years 1958 and 1959. The average amplitude 
and phase of the 24 hr component of the variation were obtained at each 
station and the relative values were compared with the values predicted by 
two types of modulation of the primary spectrum. In the first case a 
modulation of the type 4n(P)/n(P)=a. P-! was considered where P is the 
magnetic rigidity of the primary particle. The best account of the main 
features of the observations was given if the modulation was effective for 
rigidities in excess of a cut-off value which averaged 15 ev, the value of a 
being 0-236, and if the direction of maximum modulation was 79°, to the east 
of the sun—earth line. 

In the second case the primary spectrum was considered to be modulated in 
the manner suggested on theoretical grounds by Elliot (1960). 

Neither model gives complete agreement with observation. 


§ 1. INTRODUCTION 


In recent years a number of workers (e.g. Dorman 1957) have attempted 
to explain the diurnal variation in cosmic rays, which is observed by 
monitors at ground level, in terms of an anisotropy produced by some 
means in the region outside the earth’s magnetic field. In this paper an 
attempt is made to find out to what extent two simple types of modulation 
of the primary cosmic ray spectrum can account for the solar diurnal 
variations in cosmic-ray intensities observed at Makerere, Hermanus 
and Herstmonceux during 1958 and 1959. 

These stations were chosen for comparison because they differ so widely 
in latitude but lie in a comparatively narrow band of longitude. Also 
the equipment used at Makerere is closely similar to that at Herstmonceux. 


§ 2. OBSERVATIONS 
2.1. Meteorological Effects on Cosmic-ray Intensities at Makerere 


A standard neutron monitor and standard cubical meson telescopes have 
been in operation at Makerere College (geographical coordinates 0°20-2'N, 
32° 33-8’ E, conventional geomagnetic coordinates — 2° lat., 101-4° long.) 
since the beginning of 1958. The neutron monitor has been calibrated 
from time to time with a Ra—Be neutron source and has shown a long-term 
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stability in sensitivity within 0:5%. The data from both monitors have 
been corrected in the usual way for variations in barometric pressure. 
The locally determined pressure coefficient for the neutron monitor is 
—0-71% per millibar, and for the meson telescope — 0-:16°% per millibar. 
Some care is needed in deriving the barometric correction to cosmic-ray 
data at an equatorial station. This is because the atmospheric tides are 
the predominant cause of pressure variations in the tropics, and the 
barometric pressure, as a result, has a well-marked diurnal and semi-diurnal 
rhythm. Changes in the mean daily pressure from day to day are usually 
less than the range of the diurnal variations in pressure. The barometric 
coefficients at Makerere were determined by means of a partial correlation 
between the daily mean pressure at Makerere, the daily total recorded by 
the cosmic-ray monitor at Makerere, and the daily total, corrected for 
barometric effects, recorded by the Herstmonceux neutron monitor. The 
use of the Herstmonceux data in the analysis enabled us to take most of 
the primary variations into account. Data for a period of over one year 
were used. 

A spurious component of diurnal variation can remain in the pressure 
corrected data if a barograph is used for measuring the pressure. This is 
because there may be a systematic time-lag in the pressure measurements 
due to pen friction. During 1958 pressure data from a local meteorological 
station were employed at Makerere and these, being derived from readings 
of a mercury barometer, should be free from this error. 

Another meteorological factor whose diurnal variation is greater, as a 
rule, than the variation in daily mean at Makerere, is the local temperature. 
In the analysis of the daily variation of the pressure-corrected neutron 
component at Makerere in 1958 it was found that the amplitude (0-43%) 
seemed to be greater than that found at other equatorial stations such 
as Huancayo (0-22°%), and also rather anomalously larger than the 
amplitude of the variation of the meson component at Makerere. Though 
Simpson (1953) has shown that a standard neutron monitor has only a very 
small local temperature coefficient, it was suspected that the Makerere 
neutron monitor may have been affected by a temperature coefficient 
possibly mainly due to effects on the electronic recording equipment. The 
range of local temperature in the cosmic-ray laboratory could be as high 
as 6°C on an exceptional day and on the average the daily range was 
between 3 and 4°c. Careful tests using a standard neutron source and 
artificially increasing the temperature range in the laboratory showed 
that the temperature effect on counting rate was of the order of | 
(—0-:02+0-01% per °c). This effect was such as to reduce the observed 
diurnal variation below its true value but no correction has been made 
for this. The geiger-counter trays were enclosed in thermostatically 
controlled enclosures, and the efficiencies of the geiger counters for meson 
detection should be close to 100% over their sensitive volume, so long 
as the supply voltage remains in the middle of the counter plateau. There 


should: not, therefore, be any instrumental temperature effect on the 
meson counting rate. 
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Annual mean solar diurnal variations in cosmic-ray intensities 1958 and 1959. 
(a) Herstmonceux neutrons, (b) Makerere neutrons, (c) Hermanus 


neutrons, (d) Makerere mesons. 


Fig. 2 
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Harmonic dial vector-sums of the monthly, mean 24 hr component of the 
solar diurnal variations 1958 and 1959. (a) Herstmonceux neutrons, 


(b) Hermanus neutrons, (c) Makerere neutrons, (d) Makerere mesons. 
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Dorman (1957, 1959) has explained how the diurnal variation of the 
meson component of the cosmic radiation should be affected by diurnal 
temperature changes in the atmosphere. He has shown how corrections 
may be worked out from a knowledge of the heights of the isobaric levels. 
During most of the period of this work the only place where atmospheric 
soundings were undertaken in East Africa was at Nairobi, over 350 miles 
from the cosmic-ray station. As a rule flights were made only once every 
24 hrs at night. On selected days two flights were made but always at 
the same times. It is therefore not possible to know in detail how the 
heights of the isobaric layers vary during the 24 hrs. 

Quenby and Thambyahpillai (1960) have worked out the temperature 
effect at Huancayo. It has a maximum at about 05-30 hrs local time and 
amplitude 0-11%. This has been used in correcting the Makerere data, 
in tables 1 and 2. 


2.2. Corrected Diurnal Variations and Harmonic Analysis 


The annual mean diurnal variations, corrected for barometric effects, 
in the neutron and meson components at Makerere and the neutron 
components at Herstmonceux and Hermanus are shown in fig. 1. Figure 2 
shows vector sums of the monthly averages of the 24 hr waves for the 
same components. In both diagrams the times are local solar times. In 
computing the data for these diagrams some days, including those at the 
beginning of Forbush-type decreases, or showing pronounced magnetic 
activity have been omitted. Sandstrém (1959) has shown that this does 
not in general greatly affect the long-term average amplitude and phase 
of the 24 hr wave. 

Table 1 summarizes the results of harmonic analysis for the 24 hr and 
12 hr wave for the components considered. The errors given after the 
amplitudes of the variations have, for the neutron component, been 
taken as 1-3o where ois the standard deviation calculated from the counting 
statistics (McCracken 1958). 

The Makerere meson data are shown in table 1 both with and without 
the correction for the temperature effect. The errors shown are standard 
deviations calculated from counting statistics only. 

The following features of the results may be noted: 

1. There appears to be a greater variability in the amplitude and phase 
of the 24hr component at the middle latitudes (Herstmonceux and 
Hermanus) than at the equatorial station (Makerere). This is revealed 
by fig. 2. 

2. The semi-diurnal component is larger and more significant at 
Makerere than at the other stations. It does not appear to be significant 
at all equatorial stations; Huancayo shows a much smaller effect 
(0-06%). Katzman and Venkatesan (1960) have argued that the entire 
semi-diurnal variation of neutrons observed at some equatorial stations 
is due to residual barometric effects. This conclusion is disputed here 
for two reasons. Firstly, a significant semi-diurnal component in the 
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cosmic-ray intensity is found at Makerere for both neutron and meson 
components, and the ratio between the amplitudes of the semi-diurnal 
and diurnal variations is similar for the two components despite the 
difference in pressure coefficients. Secondly, the phase of the maximum 
of the residual semi-diurnal variations in the neutrons at Makerere differs 
from that of the minimum of the semi-diurnal pressure wave by almost 
90°. This corresponds to a time interval of 3 hrs, which is more than 
could be accounted for by possible barograph errors. 

3. The amplitude of the 24 hr wave at Herstmonceux was larger during 
1959 than in 1958. The change in amplitude at Makerere was smaller but 
in the opposite sense. This indicates a change in the rigidity or energy 
dependence of the modulating mechanism responsible for the variation. 

4, The time of maximum of the 24 hr wave at Makerere is more than 1 hr 
earlier than the time of maximum at Herstmonceux. The time of 
maximum of the 24 hr wave at an equatorial station produced by a 
modulating field system at a fixed orientation with respect to the 
sun-earth line is sensitive to the rigidity or energy dependence of the 
modulation. At latitudes above 50° the dependence is much less, 
owing to the focusing effect of the earth’s magnetic field for particles of 
rigidities between 10 and 30@v which arrive at these latitudes. This has 


already been pointed out by a number of observers including Dorman 
(1957). 


§ 3. CALCULATION OF THE EFFECT OF PRIMARY MODULATIONS 
3.1. Effective Aperture of Cosmic-ray Monitors 


In attempting to see how far the daily variations in cosmic-ray intensities 
recorded at different stations can be related to a common mechanism it is 
necessary to investigate the effect of the angular aperture over which a 
particular type of monitor accepts the incoming particles, and also to 
consider the effect of the earth’s magnetic field on the trajectories of the 
particles. 

The angular apertures of the monitors were considered first. In the 
case of the neutron monitor it was assumed that the angular aperture of 
a monitor deep in the atmosphere is determined only by the atmospheric 
absorption, and that the secondary particles detected by the monitor 
retain the direction of motion of the incoming primary as they pass 
through the atmosphere. Exponential absorption of the nucleonic 
component was assumed, and the characteristic thickness in which the 
intensity is reduced by a factor l/e was taken as 140g/sq.cm. The 
resultant sensitivity per unit angular range at the top of the atmosphere 
of a neutron monitor at sea level is shown in fig. 3(A). It is clear that 
this is only a first approximation to the correct angular aperture, since 
there must be appreciable scattering of the particles which tne the 
nuclear cascade. One would expect this to be more important for particles 
arising from the interactions of the low energy primaries since at high 
energies relativistic effects limit scattering angles in the laboratory frame 
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of reference. It will appear later that the diurnal variation appears to 
affect mainly those primaries whose magnetic rigidities exceed 15 av, 
For such particles half the nucleons emitted from elastic nucleon-nucleon 
collisions must lie within a cone of half-angle of about 20° about the 
direction of the incident particle. Particles emitted from inelastic 
collisions, and which are sufficiently energetic to propagate the cascade, 
are likely to be emitted in a narrower cone. 


Fig. 3 
a (a) ab 
E(8) £8), 
Arbitrary Arbitrary 
Units 5 Units , 
! 1 r 


60 10 a 30 40 so 60 
_ 8deg 9 deg 


Acceptance apertures of cosmic-ray monitors at sea level. (a) Neutron monitor, 
(6) Meson monitor. For an isotropic distribution of particles at the 
top of the atmosphere the number of counts due. to particles incident 
at angles between @ and 6+d6 to the zenith, is proportional to H(6) dé. 


It is reasonable to expect that it will be the first few interactions of 
obliquely incident primaries near the top of the atmosphere which will 
contribute most to widening the effective aperture of a monitor at sea level. 
Secondary particles from such interactions which travel in directions 
nearer the vertical than the incident primary will contribute more to the 
monitor counting rate than those preserving the original direction, because 
of their reduced absorption. A change of direction late in the nuclear 
cascade will be less effective because of the greater atmospheric absorption 
before the scattering occurs. As a rough estimate one might expect 
angles of up to 40° between the primary and secondary trajectories. 
Simple arguments show that the amplitude of the diurnal variation should 
be multiplied by a factor of the order of sin d/¢ where ¢ is a typical angle 
between primary and secondary trajectories. When ¢=40° this factor 
is ~0-9 and the observed amplitude will be 10% less than if there were 
no atmospheric scattering. 

This estimate shows that the recorded diurnal variation is in any case 
comparatively insensitive to the exact form of the effective aperture of the 
monitor until it becomes comparable with 27. 

In the case of cubical meson telescopes it was assumed that the angular 
aperture was as described by Parsons (1957) and shown in fig. 3 (b). 

In the following calculations these acceptance curves were approximated 
by dividing the solid angle over which particles are accepted into sections 
centred at 0°, 16°, 32°, and 48° to the vertical (zenith), and an average 
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sensitivity was taken for each section. The geomagnetic effect on the 
particle trajectories varies not only with the angle the incident particle 
makes with the zenith, but also with the azimuthal angle. A rather crude 
allowance has been made for this by dividing each section of the acceptance 
cone centred round one of the standard zenith angles into four parts, 
centred on the azimuth directions N, S, E and W and taking the central 
trajectories as typical. 


3.2. Geomagnetic Effects 

The complicated problem of transforming the distribution of particle 
trajectories which affect a monitor deep in the atmosphere was thus 
approximated and simplified to the problem of tracing a few standard 
trajectories to infinity, through the earth’s magnetic field. This was 
done with the aid of the curves presented by Brunberg and Dattner (1953) 
and set out in a review article by Singer (1958). These curves give the 
two coordinates of the asymptotic velocity vector; the north latitude 
dx, and 7, the longitude angle measured towards the east from the 
observer’s meridian plane. 

The yield functions worked out by Webber and Quenby (1959) were 
used in calculating the contributions to the counting rates due to particles 
of different rigidities. 

Hence it was possible to calculate the contribution to the counting rate 
of a particular monitor made by one primary particle of magnetic rigidity P 
incident on the earth’s magnetic field from an asymptotic direction 
determined by ¢y and y,. In order to proceed further and attempt to 
locate the asymptotic direction of the source of the modulation of intensity 
of the primary cosmic rays which produces the daily variation observed 
by the monitors at sea level, it is necessary to make some assumption 
about the rigidity dependence of the modulation. 


3.3. Choice of Primary Modulation Spectrum 
Dorman (1957) following Alfven (1950) suggested that the main 
contribution to the modulation arises from streams of ionized matter 
proceeding outward from the sun and carrying frozen-in magnetic fields. 


If it be assumed that the undisturbed primary cosmic-ray differential 
rigidity spectrum is of the form 


n(P)dP=AP~’dP per unit solid angle . . . . (1) 


and supposed that for particles of particular magnetic rigidity P the 
rigidity increases by AP on crossing the stream, then it can be shown 
that the effects due to deflection of the trajectories, and the change in 
rigidity on crossing the stream, combine to increase the intensity in the 
range of rigidity P to P+dP in such a way that 


An(P) _ (AP 
n(P) = const. 3 a eee ee © 


If further, as in the simplest model, the electric field in the stream increases 


Solar Diurnal Variation of Cosmic Rays during 1958 and 1959 581 


the energy of all the particles by the same amount AK, a simple 
transformation shows that the change in rigidity AP will also be 
almost constant. 

Over most of the relevant part of the cosmic-ray primary spectrum 
the relative change in intensity is therefore approximately given by 


An(P) a 

AP) > P SE ooo). x) et (3) 

‘There is some empirical evidence in favour of this modulation spectrum 
in that it seems to be of the right type to explain long-term changes in 
average cosmic-ray intensities, though the long-term changes do not seem 
to need the low rigidity cut-off found necessary in attempting to account 
for the diurnal effect (McDonald and Webber 1959). 

Additional calculations were made using the modulation spectrum 
worked out by Elliot (1960). Elliot assumed that the diurnal variation 
arises from the action of a solar dipole field. 


3.4. Comparisons with Observations 


Relative total undisturbed average intensities for the various monitors 
were calculated using the modified primary spectrum at solar maximum 
given by McDonald (1959) and geomagnetic cut-offs based on those given 
by Quenby and Webber (1959) for vertical incidence. At the equator 
the cut-off varies markedly with angle of incidence and this was allowed 
for. The primary spectrum was combined with the yield functions and 
angular acceptance factors mentioned above to determine the relative 
counting rates of the monitors. 

Dorman (1957) has shown that, in order to account for the amplitudes 
and phases of the 24 hr component of the diurnal variation of cosmic rays 
seen by monitors at different latitudes, using the P~* modulations, it is 
necessary to apply to the modulation spectrum An(P) a low rigidity 
cut-off higher than that determined by the earth’s magnetic field at middle 
or high latitudes. He attributes this cut-off to the finite thickness of the 
corpuscular streams which produce the modulation. - Particles with less 
than a certain minimum rigidity will be reflected from the streams with 
unchanged rigidity. Invoking Liouville’s theorem for such particles he 
predicted that their intensity at the upper boundary of the atmosphere 
and their isotropy will not be affected by this process. The required 
cut-off rigidity is determined empirically so that the predicted ratio of 
the amplitudes of the 24 hr variation in intensity of cosmic radiation at 
different stations is as near as possible to that found experimentally. 
The same cut-off is then applied at other stations and predictions compared 
with observations. 

In the present work the amplitude and phase of the 24 hr wave in neutron 
intensity at Herstmonceux in 1959 were taken as standard. It was found 
that the other data could best be accounted for by assuming a low-rigidity 
cut-off of 15 av for the modulation. The results are rather insensitive to 
_ the value of this cut-off within 2 or 3ayv. 
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It was supposed that the field system giving rise to the 24 hr component 
applies to the intensity of particles arising from asymptotic directions 
whose longitude measured to the east from local zenith is #,, a modulation 
of the form : 

An(P, t, tp) =(A/P) [cos (15(¢— 12) + py —a)n(P, Pe) ] (4) 
where ¢ is the local solar time in hours and « is the angle between the 
sun—earth line and the effective centre of the modulating field system 
(see fig. 4). When’ expression (4) is multiplied by solid angle and yield 
functions, and integrated over all values of P above cut-off, and over 
the result is of the form 


AN|N = const. cos [15(t— 12) +4, —«] 2 SA roe 
where 7, is an effective average value of ¥,, for the particular monitor and 
the particular type of modulation spectrum chosen. The maximum 


Fig. 4 


To Effective Centre 
of Modulating Field 


Diagram showing the relation of the direction of maximum modulation to the 
sun-earth line and local vertical. 


amplitude of the variation will be achieved when the argument of the 
cosine is zero, i.e. when 


t=tmax= 12+ SF hrs, . . . . . . (6) 


ii, was found for each monitor as follows. 


The counting rate, in arbitrary units, produced by the modulating 
spectrum of primaries was worked out for the primaries arising from each 
ten degree longitude belt in ,. 

The totals in each longitude belt were then multiplied by a factor 
cos (7, +6) where @ was arbitrarily chosen, and the results summed. 
This is a numerical approximation to an integration. The summation 
was done for several different values of 6. From the results the value of 0 
which makes the sum a maximum was found by interpolation. This 
value of 6 equals %,. The constant was found from the maximum 
amplitude of the variation at Herstmonceux,. 
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The value of #,, thus found was substituted in eqn. (6) together with the 
value of tmax for the Herstmonceux neutron monitor, found experimentally, 
in order to find a. 

Using eqn. (6) for the other monitors, with the appropriate value 
of Das and the value of « found as stated above, expected times of maximum 
at the other monitors were calculated. In these calculations no account 
was taken of the possible variation of the modulation with celestial latitude. 
Dorman (1957) gives an estimate of this which shows that there may be 
a significant reduction in amplitude at angles greater than 20° to the 
ecliptic. 

The inclination of the earth’s axis of rotation and the axis of its magnetic 
dipole to the ecliptic was also neglected. If this is taken into account, 
and if the modulation amplitude varies as the cosine of the celestial latitude, 
the calculated values of the mean diurnal amplitudes over a solar year 
should be reduced by some 7%. 

In table 2 the predicted values of amplitude and times of maximum 
predicted by both types of modulation of the primary spectrum are given. 


§ 4. DiscussIon 


It is difficult to assign standard errors to the predicted values of 
amplitude and time of maximum contained in table 2. The estimation 
of these values involves a number of stages and a number of approximations 
whose accuracy is not easy to determine. In particular the earth’s magnetic 
field has been assumed to be that of a simple dipole when working out the 
deflections of cosmic-ray trajectories. It is certain that this is not a good 
approximation, particularly in the neighbourhood of Hermanus (Quenby 
and Webber 1959). Errors in the times of maximum can also arise through 
the grouping of trajectories at an average ys, which in fact must be spread 
over a finite range. Dorman (1957) estimated the accuracy of his 
predictions of the time of maximum of the daily variation by a similar 
method at +0-5 hr. If we assume that the error can be approximated by 
a variation of random phase superimposed on the calculated value, the 
error in the predicted amplitude corresponding to a phase error of 0-5 hr 
would be in the region of 10%. 

Another difficulty, in comparing the Hermanus data with those from 
the other stations, arises from the fact that the average latitude of origin 
éy for particles affecting a monitor is different at Hermanus from what 
it is at Herstmonceux and Makerere. The majority of particles responsible 
for the daily variation at Makerere and Herstmonceux, according to our 
assumptions, arise within a few degrees north and south of the plane of the 
magnetic equator. Particles incident at Hermanus arise, on the average, 
some 20° north of the equatorial plane, and cross the equatorial plane 
on their path to the earth. The average modulating amplitude for these 
particles will therefore be some 6% less than that which is effective for 
Herstmonceux and Makerere, if the amplitude falls off as cos ¢. This 
has not been allowed for in the tabulated values. 
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With these considerations in mind it is possible to make a critical 
examination of table 2. 

Agreement between the observations and the amplitudes calculated 
from the P+ modulation spectrum is better for 1958 than for 1959. 
However, the calculated amplitude at Makerere is smaller, as compared 
with that at Herstmonceux, than is observed. Only relative amplitudes 
are significant in this case. 

In the case of the spectrum proposed by Elliot (1960) absolute values 
of the variation are predicted. The previous discussion of the effect of 
atmospheric scattering and of the inclination of the axis of the magnetic 
dipole of the earth to the ecliptic shows that the tabulated amplitudes 
should perhaps be reduced by between 15 and 20%. This still leaves a 
significant discrepancy with the Herstmonceux data and with the relative 
amplitudes at the three stations. 

It may be that both the modulation spectra which have been considered 
are weighted too strongly in favour of low rigidities. Below about 30 av 
a monitor at 50° Nis sensitive predominantly to primaries arising within 
10° of the geomagnetic equator and scans a latitude belt very similar to 
that scanned by an equatorial monitor. Above about 30 ev this similarity 
ceases, and the relative amplitudes of the variations at the two stations 
become more dependent on the ¢-dependence of the primary modulation, 
which could well reduce the relative amplitude at 50°N as required by 
the data. 

On the other hand the fact that the average times of maximum for 
neutrons differ by over one hour indicates that the mean rigidity cannot 
be greatly in excess of 30 Gv. 

Evidence from underground experiments (Sandor 1960) confirms that 
primaries of 30 Gv or more are subject to diurnal variations to an appreciable 
extent. 

Since the period surveyed in this study, the diurnal variation at Makerere 
has decreased in amplitude and the maximum has moved to an earlier time. 
The average amplitude between April 1959 and April 1960 was 0-35% 
and the time of maximum was 13:1 hrs. This is consistent with a reduction 
in the mean rigidity of the modulated particles, as solar activity declined. 

Any adequate account of the daily variation must of course take account 
of the great variability of the effect from day today. In discussing average 
values over periods of one year, one makes the assumption that the 
disturbances which produce the variability average out, but this is almost 
certainly an oversimplification. 
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ABSTRACT 


The appearance of ‘ zig-zag’ dislocation lines during the decomposition 
of crystals of lead iodide in the electron microscope can be explained by 
the climb of existing dislocations into a flattened helical form. The simple 
zig-zags are often distorted into cusped and looped forms and it is shown how 
these distortions may be accounted for by the elastic interaction between 
neighbouring segments of the helices. Such internal re-arrangement within 
helices is possible only as a result of the marked anisotropy of climb and glide 
processes in these particular crystals. The climb of arrays of dislocations 
leads to interactions between neighbouring helices. Some examples of 
such interaction are described and it is shown how these may lead to the 
formation of isolated loops of dislocation and double helices. An examination 
of the nature of loops formed in this way suggests an interesting possibility 
for the creation of new regenerative sources of glissile dislocations. 


§ 1. LyrRopucTION 


Some preliminary observations on the mode of decomposition of lead 
iodide crystals under irradiation in an electron microscope have already been 
described in a previous paper (Forty 1960). As reported there, it is possible 
to classify the events leading to complete decomposition according to a 
number of stages in the process, and the first of these is one of structural 
re-arrangement within the crystal. At this stage the distribution and 
configuration of arrays of dislocations already existing in the crystals 
changes in detail and, at the same time, new loops of dislocation begin to 
appear. The configurational changes can be accounted for by considering 
that the dislocations are climbing out of their glide planes and this, together 
with the observations of loop formation, may be taken as evidence that 
point defects are being created within the crystal during the electron 
bombardment. The climb of the dislocations appears to proceed along 
certain preferred planes. The orientation of these planes in the crystal 
structure suggests very strongly that the point defects involved in the 
process are molecular vacancies (i.e. PbI, vacancies or associated complexes. 
of lead and iodine ion vacancies) which are probably created thermally as a. 
result of the heating of the crystal by the electron beam. 
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One of the most interesting features of the initial phase is the formation of 
‘zig-zag’ dislocation lines. Simple examples of these have already been 
described and discussed in the previous paper. They are thought to be the 
projected images of segmented or flat helices whose axes lie in the plane of 
the crystal, and which have been developed by the climb of dislocations 
originally lying in this plane. and predominantly in the screw orientation. 
In general, the forms of these helices are more complex than those which 
might be expected for more isotropic crystals, and it is the main purpose of 
this paper to describe the various forms which have been observed and to 
attempt to relate these to the crystallographic features of climb and glide 
processes in the hexagonal structure of lead iodide. 

The basic unit in the lead iodide structure may be described as a 
sandwich of divalent lead ions between two close-packed layers of iodine 
ions. The hexagonal structure is built by stacking these sandwiches in the 
e-direction. In this paper, however, we find it most convenient to describe 
the structure as a hexagonal close-packing of iodine layers with the smaller 
lead ions filling one-half of the octahedral positions between alternate 
layers. Thus, all planes and directions in the lattice may be indexed 
simply in the h.c.p. system, a procedure which will be adopted throughout 
the following discussion. One feature of the structure which deserves 
emphasis at this stage, since it relates the observed climb planes with 
the production of molecular vacancies, is the observation that low-index 
planes having molecular composition are limited to two kinds. These 
are the basal planes, which contain PbI, sandwiches, and the first 
order pyramidal planes {1011}, which are followed closely by slightly 
corrugated Pbl, layers. The climb of the dislocations appears to be 
restricted to follow these planes and this might reflect the ease with which 
molecular vacancies can be created or movedinthem. Moreover, in such a 
structure easy glide occurs in one plane only, namely the basal plane where 
glide presumably proceeds most readily between the neighbouring iodide 
layers which are not separated by the interstitialleadions. The restrictions 
on climb and glide of dislocations thus imposed by the crystal structure 


can account for many of the configurations observed and described in 
the following section. 


§ 2. Tot OBSERVATIONS 


Techniques for preparing the crystals and detecting dislocations in them 
have been described already in the earlier paper (Forty 1960). It is 
convenient, however, to give a further brief account here. Crystals thin 
enough for direct transmission electron microscopy can be grown from a 
solution of the re-crystallized salt in aqueous alcohol by just saturating the 
solution at 100°c and allowing this to cool slowly to room temperature. It 
is found that, whereas the thicker crystals precipitate to the bottom of the 
container, those which are thin enough for this work remain floating as 
very thin platelets on the surface of the solution. These can be collected 
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on a microscope specimen grid for direct observation without further 
preparation. The observations described here have all been made with 
specimens prepared in this way and examined in the standard Philips 
electron microscope operating at 80 kv. | 

{t is found that, although a few of the crystals appear to be free of 
dislocations, many of them contain extensive arrays which appear to lie 
in the plane of the platelet and have Burgers vectors lying in this plane 
(i.e. in the basal plane). These are introduced by deformation during 
mounting. Clearly the geometrical conditions are ideal for observing any 
climb out of the basal plane which may occur to modify the shapes of the 
dislocations. The geometrical conditions are also extremely favourable 
for observing glide processes in these crystals; a feature which can be 
exploited and which will be discussed more fully in a later paper. 

After a short period of exposure to the electron beam the crystals begin 
to undergo some structural rearrangement. At this stage the dislocations 
lose their initially smooth appearance to take up more complex configura- 
tions. Figure 17, for example, shows part of a crystal in which a simple 
network of dislocations has become distorted in this way. It will be noticed 
that the basic feature in this pattern is the zig-zag. Similar shapes have 
been observed by Thomas and Whelan (1959) in their investigation of the 
climb of dislocations in quenched aluminium/copper alloys. These have 
been attributed to the formation of helical dislocations by the condensation 
of vacancies on screw dislocations, and it seems reasonable that a similar 
interpretation might be applied to our present case. As pointed out by 
Thomas and Whelan, a simple circular helix viewed at right angles to its 
axis should be projected as a sine waveform. Variations from sinusoidal 
form can be expected if the helix is not circular or if it is not viewed at right 
angles. The extreme sharpness of the zig-zags in lead iodide in fact lead 
us to the conclusion that, if they represent helices, these must have a 
flattened cross section with long segments lying in the basal plane and 
shorter segments lying in some plane, probably one of the {1011} planes, 
inclined to this. Such a form is drawn in perspective in fig. 4 (a) where the 
different thicknesses of line are intended to represent segments at different. 
levels in the crystal. This conclusion is also necessary on geometrical 
grounds if helices having projected widths equal to those of the zig-zags 
(i.e. about 5000 4) are to be accommodated in a platelet whose thickness is 
only a fraction of this. Evidence that the zig-zags do, in fact, represent 
flattened helices with segments at different levels in the crystal rather than 
zig-zag dislocations confined to a single plane can be drawn from an 
inspection of contrast effects in the image. This is discussed later in 
connection with figs. 5 and 6. 

Figure 2 shows further examples of zig-zags. Again, these have the 
sharpened form which suggests that they represent flattened helices. 
There are, however, some additional features of interest in this photograph. 


+ Figures 1, 2, 3, 5, 7, 8, 10, 11 and 12 are shown as plates. 
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In the areas marked A, for instance, the simple zig-zag form is replaced bya 
cusped shape. Other examples of this cusped form are to be found in the 
photographs illustrating this paper. In particular, a number of good 
examples can be seen in the helices within region Ainfig.3. The appearance 
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Illustrating the distortion of a flattened helix to produce a cusped zig-zag in 
projection. (a) Perspective view of flattened helix. (6) Projection of 
flattened helix into simple zig-zag. (c) Cusped zig-zag produced by 
attraction between neighbouring long segments of the helix. (d) 
More accurate description of cusped zig-zag, showing the projected 
portions of the short segments. (e) Perspective view of cusped helix. 


of cusps on one side of a zig-zag can be expected if the helix is viewed at an 
angle, and many examples of this kind have been published by Thomas and 
Whelan. Cusps on both sides, however, cannot be accounted for in this 
way and a different explanation providing some real distortion of the helix 
must be sought. Figure 4 illustrates one possible mechanism. A flattened 
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helix has its long segments lying in the basal plane of easy glide whilst its 
short segments lie in inclined pyramidal planes in which glide is not possible. 
Thus, although the long segments are free, the helix is effectively pinned 
against translation along its axis by the short segments. Therefore only 
internal rearrangement is possible. Adjacent long segments of the helix 
have opposite Burgers vectors and since they lie in parallel glide planes they 
tend to glide by mutual attraction into a more stable relationship to lie 
closer together. This leads to a distortion of the long segments, giving 
rise to the cusped appearance in the projected image. 

It is difficult to calculate the final shape of the helix resulting from this 
kind of internal interaction or even to estimate the magnitude of the forces 
experienced by the long segments. An approximate description might be 
found by considering a pair of neighbouring segments as a pair of edge 
dislocations lying in parallel glide planes and having opposite Burgers 
vectors. If the dislocation lines were continuous and the structure was 
isotropic then a position of stability against further glide would be reached 
when their separation in projection equalled that between their glide planes 
(see, for example, a discussion of this kind of situation by Read (1953)). 
This conclusion must, of course, be modified by taking into account the 
presence of the short segment terminating the long segments, the fact that 
the long segments are not continuous and are not pure edge dislocation 
but have a mixed character, and also the crystal anisotropy. However, 
refinements of this kind are not likely to have a significant effect on the 
final shape. 

A closer inspection of the cusped zig-zags shows that the tips are not 
perfectly sharp but that the segments are separated here by a distance of a 
few hundred angstrom units, where the dislocation line becomes less visible. 
According to the present discussion this should represent the portion 
occupied by the short segment and gives an order of magnitude for the 
separation between the basal planes containing the long segments. A 
value of a few hundred A is consistent with the geometrical restriction 
imposed by the crystal thickness. It appears, therefore, that a more 
- accurate description of the cusped helix should be represented by fig. 4 (e), 
giving a projected zig-zag image as shown in fig. 4 (d). 

A further interesting feature of fig. 2 is shown in the area marked B. 
In this region the normal zig-zag has been modified to one in which the 
tips are terminated by chains of loops. Similar examples are shown in 
fig. 5 at the places marked A and C. This feature has been reproduced in 
fig. 6 where the different thicknesses of line are intended to represent the 
different depths at which the dislocations lie in the crystal. The variation 
in depth is indicated by the variation in contrast of the dislocation image in 
the photograph; since the zig-zag must be a continuous dislocation 
having the same Burgers vector everywhere, and there is only one set of 
reflecting planes contributing to the dislocation contrast, variation in 
depth is the only factor which can change the contrast in this way. The 
extent of this variation in contrast depends on the particular diffracting 
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conditions, on the positions of the dislocations, and on the thickness of the 
crystal. This explains why the variation is more pronounced in some of the 
photographs shown here than in others. , 

‘A number of mechanisms leading to the formation of chains of loops have 
been examined, but none of these is completely satisfactory. The most 
attractive is illustrated in stages by fig. 6. It is first necessary to consider 
the formation of single loops at the tips of the normal zig-zag. This stage 
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Illustrating the mechanism by which chains of loops may be formed in helices. 
(a) Perspective view of flattened helix in which the short segments lean 
backward. (6) Projection of flattened helix. (c) Interaction between 
neighbouring long segments, leading to loop formation. (d) Growth of 
chains of loops from looped helix after the outward climb of the short 
segments. (e) Growth of chains of loops from cusped helix after the 
outward climb of the short segments. 


in the process is frequently observed and examples are visible in the area 
marked B in fig. 5. Clearly both loops and cusps are possible in helices in 
these crystals and it is thought that the choice is determined by the orienta- 
tion of the climb plane in which the short segments lie. If the short 
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segments lean forward along the helix cusps are formed (as in fig. 4) ; if they 
lean backward (as in fig. 6) loop formation is favoured since neighbouring 
long segments must always overlap to some extent. The length and apparent 
width of the loops in fig. 5 is surprising in view of the strong attraction which 
must exist between neighbouring long segments and which has already been 
used to explain the cusping in some of the helices. It would appear that 
some additional shear stress must operate in a systematic manner in the 
vicinity of the loops in order to hold them further apart. It is considered 
that such a stress might be accounted for by the observation that the turns 
of the helix are not pure edge dislocation but have some screw character. 
Our elementary model of a pair of edge dislocations lying in parallel glide 
planes must therefore be modified by the addition of a screw dislocation 
orthogonal to these and lying in a parallel glide plane between them. An 
inspection of such a model with the appropriate choice of Burgers vectors 
suggests that the equilibrium separation of the edge components can be 
increased by the addition of the screw dislocation. Thus the representation 
of a looped helix given in fig. 6(c) appears to be a permissible form of 
distortion produced by elastic interaction between neighbouring turns. It 
is interesting to note that the crossing over of the segments to develop the 
loops is of the symmetrical form expected for a pair of long dislocations of 
this kind, and which has been discussed previously by Read (1953) and 
Amelinckx (1956). 

The last stage in the process is indicated in fig. 6 (d) which shows how the 
final appearance of a chain of loops in the projected image can be produced 
by the glide of the two segments of the loop towards each other to link at a 
few isolated points. This leads to a stable situation but there is no obvious 
reason why the attraction should exist to initiate this localized glide. <A 
more satisfactory explanation of the final stage follows if the short segments 
of the helix climb further outwards in the allowed {1011} planes to extend 
the loops. The segments of the loops then separate to the equilibrium 
spacing except at a few places where they remain pinned together by some 
isolated attractive force. Pinning might well result from. localized 
precipitation of lead on the dislocations during the later stage of climb. 
Some further observations which will be described in a later paper show that 
precipitation of lead formed by the final decomposition of the crystal often 
occurs preferentially along helical dislocations. This explanation accounts 
for the observation that the chains of loops extend preferentially along close- 
packed directions in the basal plane (see also Forty 1960). The {1011} 
planes in which the short segments of the helices can climb to extend the 
loops intersect the basal plane along these directions. 

Chains of loops can also be formed in a similar manner from cusped 
helices if the short segments of these climb further outward. It is 
important though to note that this form does not involve a crossover in the 
long segments and therefore the alternative linkage of the chains shown in 
fig. 6 (e) must be developed. If the interpretation of contrast variations for 
the dislocation line in fig. 5 in terms of variation in level in the crystal is 
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reliable, then the chains at A are developed from looped helices (giving the 
form of fig. 6(d)) whilst those at C are grown from cusped helices (form 
6(e)). Some indication that loops may be formed from cusped helices may 
also be found in the enclosed area B in fig. 3. 

The dislocation in fig. 5 appears to be broken in the region D. This must 
be a place where the line has climbed completely out of the crystal, pre- 
sumably in an attempt to overcome an obstacle surrounded by the strongly 
diffracting patch. It is difficult to interpret the cusping and apparent 
reversal of the dislocation at the break, although this appears to be a fairly 
common mode of termination since several other examples can be seen in 
fig. 2. No satisfactory explanation is available at present. 

The climb of arrays and networks of dislocations leads to even more 
complicated configurations. These are of special interest when the helices 
produced are sufficiently close to allow interactions to occur. Some 
examples of this kind have been given by Thomas and Whelan in their 
discussion of helical dislocations in the aluminium copper alloys. They 
show how interactions between turns of neighbouring helices in the arrays 
can lead to tangled networks of helices and isolated loops. Similar 
observations have been made on the climb of arrays of dislocations in lead 
jodide and some of these will now be described in the following paragraphs. 

Figure 7 shows a fairly simple example of the climb of an array. This 
appears to be a crossed grid of screw dislocations, each of which has become 
helical. Such a configuration represents the climb of a twist boundary. 
It is interesting to note from this that if a twist boundary undergoes climb it 
loses its planar character as its component helical dislocations occupy a 
finite slab of the crystal. This is not true for a tilt boundary since the 
component edge dislocations here simply climb in the same plane. 

A more complex array of helices is shown in fig. 8. _ This includes a family 
of closely spaced parallel helices crossed by two other independent 
dislocations. In the region marked A neighbouring helices have interacted 
to form a tangled network. It is difficult in this case to analyse the details 
of the interaction although the general form of the individual helices is still 
detectable. The area B, however, provides a particularly interesting 
example of interaction which can be followed in detail. In this region two 
helices have attracted each other and at one point some overlap and 
annihilation has occurred to produce a pair of double helices and a single 
closed loop. This kind of interaction is represented in detail in fig. 9. 
Further interaction between segments of the double helices can lead to the 
formation of more closed loops, and a simple example of this kind of event 
is provided by fig. 10. This shows part of a hairpin dislocation, one side of 
which has become helical whilst the other side has remained undistorted. 
At the bend in the hairpin segments of the helix have touched the other side, 
and, after some annihilation, a closed loop has become detached. 

Interactions of this kind lead to complex distributions of loops such as 
those shown in figs. 11 and 12. In fig. 11 segments of the original helices 
are still recognizable, but in fig. 12 the annihilation has proceeded so far 
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that their identity has been destroyed completely. There is some indica- 
tion of lining-up of loops, however, and this probably occurs along the lines 
of the original dislocations. It might be noted that the loops thus formed 
by climb have Burgers vectors lying in the basal plane and can therefore be 
expected to glide or regenerate other glissile dislocations to contribute to the 
plastic flow of the crystal under suitable applied stress. This is in marked 
contrast to loops which are formed directly by the collapse of condensed 
aggregates ofvacancies. The Burgers vectors of these loops do not lie in the 
basal plane and they are therefore sessile. Loops of this kind have been 
observed in lead iodide (see Forty 1960) and, in fact, some are visible in 
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Representing the interaction at B in fig. 8. 


the photographs illustrating this paper (for example, see fig. 1, areas A and 
B, fig. 5, area E and the numerous small loops in the background of fig. 11). 
It is interesting to observe the marked difference in contrast for the two 
kinds of loops. This is clearly associated with the different orientations of 
their Burgers vectors relative to those strongly diffracting planes of the 
erystal lattice which are contributing to the microscope image in these 
regions. : 


§ 3. CONCLUSION 


These observations on the formation of helices in lead iodide are useful 
in that they draw attention to the marked anisotropy of glide and climb 
processes in crystals of this kind. As we have seen from the foregoing 
discussion, the shapes and interactions of helices are largely dependent on 
the crystallographic details of the processes involved in their formation. 

This provides a marked distinction between the observations reported here 
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and those of Thomas and Whelan on helices formed in aluminium/copper 
alloys, where climb and glide may proceed in a much more isotropic manner 
and where, in fact, the helices are generally rounded and lack crystallo- 
graphic detail. The peculiarities of the features reported in the present 
paper clearly result from the anisotropy of the hexagonal structure. It 
should be noted, however, that the anisotropy of the particular crystals 
with which we are dealing here is: almost certainly enhanced by the 
molecular composition, and it does not necessarily follow that similar 
features might be observed in other hexagonal structures. In the present 
case, the structure and the composition of the crystals restricts the glide 
motion of dislocations to follow the basal (0001) plane and the climb 
motion to this and the {1011} planes. 

One of the most interesting features of helices, both in lead iodide and in 
aluminium/copper, is the extreme stability of these dislocations against 
further glide. ‘No overall movement which might be described as glide 
has been observed, although some internal re-arrangement of helices 
in lead iodide occurs. In the case of the aluminium/copper alloys, this 
stability has been attributed to the possibility that the helices are favoured 
sites for the separation of 0’ precipitates during the ageing of these alloys, 
and that these pin the dislocations against further movement. This 
might account also for the lack of any distortion of the helices from their 
original rounded shape. There is also some evidence that precipitation of 
metallic lead formed during the later stages of decomposition occurs 
preferentially along helical dislocations in lead iodide. It is not necessary, 
however, to seek an explanation of this kind for stability in this case ; indeed, 
the distortions inside the helices indicate that some segments remain free 
to move. Stability is provided by the short segments of the helices which 
appear to lie in pyramidal planes and are therefore unable to glide. Thus 
the dislocations in lead iodide are pinned only at isolated points. Moreover, 
this pinning is effective only against glide since some further movement 
can take place by the climb of the short segments. The distortions of the 
helices from the elementary form are possible only as a result of this greater 
freedom of movement. 

Some explanation is required for the extreme flatness of the helices in 
these crystals. A flattened cross section clearly indicates that climb 
proceeds more readily along directions in the basal plane; that is, the 
climb of the short segments in these directions proceeds faster than that 
of the long segments out of the basal plane. This anisotropy in the rate 
of climb might be associated with a greater flux of vacancies through the 
basal planes or it might be determined simply by the relative capacities 
of the long and short segments for accepting vacancies. The former 
possibility is perhaps the controlling factor as it is thought that the vacancies 
are being created near the centre of decomposition in the crystal, where the 
electron beam is most intense, and flow outwards into the surrounding good. 
crystal to create the supersaturation required to promote climb. This 
_ flow can be expected to proceed more readily in the basal planes. 
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The formation of loops by the interaction of neighbouring or overlapping 
helices during the climb of an array of dislocations suggests an interesting 
means of creating regenerative sources for further plastic flow in these 
crystals. For loops created by such a process have Burgers vectors 
lying in the glide plane and are therefore able to move during further 
stressing. Moreover, as they are formed from helices, many of these loops 
must be jogged out of the basal plane in the way shown in fig. 13. The 
jogs, like the short segments of the helices, lie out of the glide plane and 
might therefore act as pinning points against the translation of the long 
segments. In this way the long segments might well become sources of 
new dislocation loops by the regenerative mechanism suggested by Frank 


Fig. 13 
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Illustrating one possible form of jogged loop created by intersection of flattened 
helices. This particular loop could be that formed by the interaction 
shown in fig. 10. 


and Read (1950) and developed further in a recent paper by Kuhlman- 
Wilsdorf (1958). Evidence for the generation of new loops of dislocation 
from existing arrays has been found in some of the crystals, but it has not 
yet been possible to resolve the details of the operative mechanism. 
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NOTICE 


X-ray Evidence for Segregation of Solute to Stacking Faults in a Copper— 
Aluminium Alloy, R. W. Cann and R. G. Daviss, 1960, Phil. Mag., 
5, 1119. 


It has been brought to our attention that K. Nakajima has observed 
an enhanced small angle scattering from cold worked and heat treated 
copper-aluminium alloy (Nakajima, Sci. Rep. Res. Insts Téhoku Univ., 
A, 1960, 12, 39), and not from copper-nickel alloy as we had deduced from 
a brief note published by him earlier (J. Phys. Soc. Japan, 1959, 14, 1825). 
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REVIEWS OF BOOKS 


Reports on Progress in Physics. Vol. XXIII. By Tur Puystcan Socrury. 
(John Wright & Sons Ltd., 1960.) [Pp. 629.] Price £3 3s. Od. 


TuIs volume contains ten stiff doses of contemporary physics—doses that will 
cure the malady of ignorance about any of the specialties covered. However, 
the volume is no panacea; the ailing physicist who cannot keep himself in- 
formed about annual progress along the whole front of physics cannot look to 
this volume for a cure. No one could reconstruct what happened in physics 
during 1959-1960 by reading the book and no one will be expected to read it 
from page | through page 629, except a reviewer. 

There was a time when the Reports on Progress in Physics informed you about 
the year’s progress in many active parts of the subject ; the spectroscopist 
could keep informed about crystallography, the nuclear physicist about mole- 
cular physics. This age has passed and the Reports have changed; they are 
no longer—as they once were—the ‘Annual Reports’ of the company of 
physicists. They no longer show signs of a consistent, comprehensive plan to 
deal with the year’s physics in a single volume. The twenty-third volume leaps 
from one corner of physics to another, boring in with vigour and virtuosity 
wherever it lights, but the volume could as easily have been ripped into ten 
separate little tracts. The contents of the volumes are as unplanned as the 
content of an issue of The Physical Review; it is a collection of miscellaneous 
review papers disguised as progress reports. 

On every side one hears physicists complaining more and more frequently 
that keeping up with the literature is a job that has become impossible. The 
Reports do little to solve the problem behind this complaint; perhaps in the 
face of the enormous variety and complexity of physics there is no longer a 
way to meet the complaint, but I wonder if the editors of the Reports did not 
abandon their early efforts too soon. The turning point seems to have been in 
the early years of World War II. Did war pressures cause an expedient course 
to be taken, a course that has not been changed back because it has been an 
easier course to follow ? Or has physics really become so much more complex 
than in the 1930’s that the editorial policies of that decade would be impossible 
in the 1960's ? 

The quality of each separate article is good—each is thorough, reasonably 
well written, and competent. There are three articles on subjects of solid state 
physics. One of these, on ‘ Group Theory in Solid State Physics ’, is really a 
short text-book on the subject. It, like each of the other nine articles, is avail- 
able separately and I would guess that it may often be used separately by solid 
state theoreticians. Another article, ‘ Experimental Analysis of the Electronic 
Structure of Metals ’, uses 30 pages to review—with a novel touch, to be sure— 
the traditional independent particle model of electrons in metals. The 
remaining two-thirds of the article develops a physical understanding of the 
four principal experimental techniques that have been used with remarkable 
success to map the Fermi surfaces of several real metals. The third article, 
‘ Band Structure Calculations in Solids ’, covers an active corner of theoretical 
physics, and this article goes over each of the principal methods of calculating 
band structures that theorists have been exploiting. 

Three articles in the field of optics are included: ‘ Optical Properties of Thin 
Films ’, ‘ Photoelectric Photometry ’, and ‘ New Developments in Interference 
Spectroscopy ’. These three articles treat the theory and practice of experi- 
mental techniques in optics. The second article is short and almost completely 
oriented to experiment; the other two are long and carry along developments of 

‘the theory of the experimental techniques and apparatus. 
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One article, ‘ Planetary Nebulae’, treats a problem in theoretical astro- 
physics: how can the characteristics of these nebulae be deduced from the 
observational data? Planetary nebulae are mentioned in almost every 
elementary astronomy text, a distinction not shared in elementary physics 
texts by many of the nine companion topics; however, this article is definitely 
for the specialist, being filled with detailed calculations based on known physical 
principles. 

The last three articles deal with topics in nuclear physics. In ‘ Orbital 
Electron Capture by the Nucleus ’, the authors review the theory of this process, 
show how the transition probabilities are to be calculated and then run through 
experimental results for each type of transition. The article, * Precision 
Measurement in Gamma-Ray Spectroscopy ’, is another contribution on the 
theory and practice of an experimental technique; this article should guide one 
to almost anything he wants to know about precision gamma-ray spectro- 
meters, including details of the leading instruments in various laboratories. 
The final report on ‘ 7He Induced Reactions ’ covers the use of *He as a nuclear 
projectile and the nuclear reactions induced by it, surveys experimental results, 
and ends with a section on ‘ Suggested Experiments ’. 

I am sure that some of these articles may become the “ Bible ’ for workers in 
the specialized fields just as have the famous articles by Bleaney and Stevens 
and by Bowers and Owen which nurtured a scientific generation of workers on 
electron paramagnetic resonance in crystals. However, I see no reason why a 
high quality annual review journal should be camouflaged by the green covers 
of a Progress Report. Nor do I see why the early purposes of this series should 
have been permanently abandoned. The Report of the Council of the Physical 
Society for the year ending February 28, 1934, has this to say: “‘ Some time 
ago the Council sought the opinion of the Fellows of the Society on the desir- 
ability of initiating a series of annual reports on the progress of physical science. 
The unequivocal approval expressed by the Fellows . . . have encouraged the 
Council to proceed with this project.” These words, inaugurating the series, 
and the early volumes of the series show that the intent was to cover the 
important new developments in physics as they occurred and were recognized. 
Lest my account of the contents of this volume still leave any reader believing 
this coverage has been achieved, I will add that nowhere in this volume did I 
find reference to the work of Méssbauer. R. W. 8. 


Modern Physics. By M.S. Smirn. (London: Longmans, Green & Co., 1960.) 
[Pp. 254.] Price 15s. ; 
THIs is a brave and, on the whole, successful attempt to introduce the ideas 
of modern physics to the intelligent but uninstructed reader. It is perhaps 
not surprising that the treatment is least satisfactory for the more highly 
developed branches of classical physics, such as statistical mechanics, which 
are not susceptible to abbreviated presentation. But the author’s accounts of 
relativity and wave mechanics are excellent, and form a good, though rather 
stiff, introduction to more detailed study. On the descriptive side the writing 
is clear; perhaps there is occasionally too much emphasis on historical experl- 
ments and outworn ideas, but this is a matter where differences of opinion 
need not be construed as adverse criticism. The book may be highly recom- 
mended to the sixth-former who really wants to now about physics, and to the 
undergraduate scientist whose speciality is not physics. Physics specialists 
will soon outgrow it, but there is very little in it which they will have to unlearn 
and this is most creditable. A.B. P 


[The Editors do not hold themselves responsible for the views 
expressed by their correspondents.]. 
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An enlarged view of a zone-refined lead specimen exhibiting the surface 
condition produced by cold-work. 


R. B. CAMPBELL and L. MULDAWER Phil. Mag. Ser. 8, Vol. 6, Pl. 66. 


Fig. 1 


Replica from (100) face of alloy annealed at 450°c for 45 min. Not shadowed. 
Approximate size of precipitate particles: 4004. The macroscopic 
etch pits are believed to be a secondary etching effect arising from the 
lining up of the precipitate particles. 


R. B. CAMPBELL and L. MULDAWER Phil. Mag. Ser. 8, Vol. 6, PI. 67. 


io, 2 
Fig. 2 


Replica from (110) face of alloy annealed at 450°c for 30min. Shadowed, 
Approximate size of precipitate particles : 400 A. 


Fig. 3 


Transmission electron micrograph of alloy annealed at 450°c for 200 min. 
The dark streaks are the cobalt precipitate. These dark regions became 
empty when the foil was preferentially cobalt etched. 


D. HULL and I. L. MOGFORD Phil. Mag. Ser. 8, Vol. 6, Pl. 68. 
Fig. 1 


Specimen irradiated for 12 hours, fine precipitation and decorated dislocation 
line. Plane of the foil is close to (111). x30 000. 


Specimen irradiated for 72 hours, irregular loop shaped particles some of which 


appear edge-on lying parallel to <110>. Note denuded zone around 
low angle boundary. Plane of foil close to (110). 30000. 


D. HULL and |. L. MOGFORD Phils Mag.-ser. 8, Vol. 6, Pl._69. 


Fig. 3 


Specimen irradiated for 72 hours, edge on loops at right angles. 
Plane of the foil is close to (100). x30 000. 


Specimen irradiated for 72 hours, network of dislocations 
decorated with precipitates. x30 000. 


D. HULL and I. L. MOGFORD Phil. Mag. Ser. 8, Vol. 6, Pl. 70. 


Specimen aged for 72 hours at 60°c, then irradiated for a further 72 hours, 


coarse irregular precipitate with fine precipitates. Plane of foil is close 
to (111). x30 000. 


Fig. 6 


Specimen aged for 72 hours at 100°C. Clustered ipita i l 
: precipitates with 6-fold 
symmetry. Plane of the foil is close to (LAT), x30 000. ‘ 


D. HULL and |. L. MOGFORD Phil. Mag. Ser. 8, Vol. 6, Pl. 71. 


Fig. 7 


Specimen aged for 72 hours at 100°c, edge-on precipitates at right angles along 
dislocation lines. Plane of the foil is close to (100). x30 000. 


Fig. 8 


D. HULL and |. L. MOGFORD Phil. Mag. Ser. 8, Vol. 6, Pl. 72. 
Fig. 8 (continued) 


(c) 
Specimen irradiated for three days followed b 
° y 1 hour at 175° ° 
and (c) 250°C. Plane of the foil is close to (111). (a) 5 (b) 200°, 


D. HULL and |. L. MOGFORD Phil. Mag. Ser. 8, Vol. 6, Pl. 73. 
Fig. 9 


(0) 
Specimen irradiated for three days (a) before heating on the hot stage in the 
microscope, (b) after heating for 3 minutes at 180°c. Plane of the foil 
is close to (110). x30 000. 


J. HOLDEN Phil. Mag. Ser. 8, Vol. 6, Pl. 74. 


(0) 


Detail of Debye-Scherrer ring structure (Fe radiation; line 16) from: (a) original 


rolled 99-8°%, aluminium sheet; and (b) fracture surface formed by slow 
fatigue crack growth in the sheet. 


J. HOLDEN Phil. Mag. Ser. 8, Vol. 6, Pl. 75. 


Transmission photographs (unfiltered Fe radiation) divergence d0=10~ from: 
(a) annealed 99-8% aluminium crystal ; (5) crystal deformed 20% in 
tension; and (c) crystal deformed in tension as for (b) and then subjected 
to 30000 cycles at e),=7x 107°. 


J. HOLDEN Phil. Mag. Ser. 8, Vol. 6, Pl. 76. 


Fig. 3 


(a) i 
2000 
(0) (e) 
156 4000 
( (f) 


450 6000 


Development of cyclic sub-grain structure in a crystal subjected to increasing 
numbers of cycles of torsion at ¢=5x 10-3. Unfiltered Fe radiation, 


(a) 0 cycles, (b) 156 cycles, (c) 450 cycles, (d) 2000 cycles, (e) 4000 cycles, 
(f) 6000 cycles. 


J HOLDEN Phil. Mag. Ser. 8, Vol. 6, PI. 77. 


(2) (0) 


Development of cyclic sub-grain structure: (a) dark ground; and (b) same 
field polarized light ( x 500). 


(a) (2) 


(a) Section of part of crack front in mild steel sheet specimen showing joining of 
numerous small cracks (x 350); and (b) leading cracks at centre of 
mild steel sheet specimen (x 900). 


A. J, FORTY Phil. Mag. Ser. 8, Vol. 6, PI. 78. 


Fig. 1 


The distortion of an array of dislocations during the initial stage of 
decomposition. x 16 000. 


The climb of dislocations into flattened helices. Showing simple Zig-Zags, 
cusped zig-zags (A) and looped zig-zags (B). x 20 000. 


vo ORD Y: Phil. Mag. Ser. 8, Vol. 6, Pl. 79. 


The climb of a family of dislocations. Note the cusped zig-zags in region A, 
and the creation of loops at the tips of cusped zig-zags in area B. 
x 7000. 


Illustrating the formation of loops in zig-zag dislocations. Note the single 
loops in region B and the multiple loops in regions A and C. x 27 000. 


A. J. FORTY Phil. Mag. Ser. 8, Vol. 4, Pl. 80. 


Biged 


The climb of a crossed grid of screw dislocations. x 7000. 


The climb of an array of dislocations, leading to interactions between neigh- 
bouring helices. Note the tangled network at A and the interaction 
between a pair of helices at B. x 8000. 


PalALOR UY Phil. Mag. Ser. 8, Vol. 6, Pl. 81. 


The formation of an isolated loop from a helical dislocation. x 20000 


Fig. 11 


Distribution of elongated loops formed by the interaction of looped zig-zags. 
The long segments of the helices are still recognizable. x 14 000. 


A. J. FORTY Phil. Mag. Ser. 8, Vol. 6, Pl. 82. 


Distribution of loops formed by the interaction of helices. The helices are no 
longer recognizable but there is some alignment of the loops along the 
original dislocations. x LO 000 


